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This chapter provides a comprehensive overview of the Bohmian formulation of quantum mechan- 
ics. It starts with a historical review of the difficulties found by Louis de Broglie, David Bohm, and 
John S. Bell to convince the scientific community about the validity and utility of Bohmian mechan- 
ics. Then, a formal explanation of Bohmian mechanics for nonrelativistic, single-particle quantum 
systems is presented. The generalization to many-particle systems, where the exchange interaction 
and the spin play an important role, is also presented. After that, the measurement process in 
Bohmian mechanics is discussed. It is emphasized that Bohmian mechanics exactly reproduces the 
mean value and temporal and spatial correlations obtained from the standard, that is the Copen- 
hagen or orthodox, formulation. The ontological characteristics of Bohmian mechanics provide a 
description of measurements as another type of interaction without the need for introducing the 
wave function collapse. Several solved problems are presented at the end of the chapter, giving 
additional mathematical support to some particular issues. A detailed description of computational 
algorithms to obtain Bohmian trajectories from the numerical solution of the Schrodinger or the 
Hamilton-Jacobi equations are presented in an appendix. The motivation of this chapter is twofold: 
first, as a didactic introduction to Bohmian formalism, which is used in the subsequent chapters, 
and second, as a self-contained summary for any newcomer interested in using Bohmian mechanics 
in his or her daily research activity. 
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Bohmian Mechanics: From Nanoscale Systems to Cosmology". The final 
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I. PREFACE 

Most of our collective activities are regulated by other people who decide whether they are well done or not. One 
has to learn some arbitrary symbols to write understandable messages, or to read those from others. Human rules 
over collective activities govern the evolution of our culture. On the contrary, natural systems, from atoms to galaxies, 
evolve independently of the human rules. Men cannot modify physical laws. We can only try to understand them. 
Nature itself judges, through experiments, whether a plausible explanation for some natural phenomena is correct 
or not. Nevertheless, in forefront research where the unknowns start to become understandable, the new knowledge 
is still unstable, somehow immature. It is supported by unclear experimental evidences or the evidences are still 
subjected to different interpretations. Certainly, novel research grows up closely tied to the economical, sociological 
or historical circumstances of the involved researchers. A period of time is needed in order to distil new knowledge, 
separating pure scientific arguments from cultural influences. 

The past and the present status of Bohmian mechanics cannot be understood without these cultural considerations. 
The Bohmian formalism was proposed by Louis de Broglie even before the standard, i.e., Copenhagen, explanation of 
quantum phenomena was established. Bohmian mechanics provides an explanation of quantum phenomena in terms of 
point particles guided by waves. One object cannot be a wave and a particle simultaneously, but two can. Specially, 
if one of the objects is a wave and the other is a particle. Unfortunately, Louis de Broglie, himself, abandoned 
these ideas. Later, in the fifties, David Bohm clarified the meaning and applications of this explanation of quantum 
phenomena. Bohmian mechanics agrees with all quantum experiments done up to now. However, it remains almost 
ignored by most of the scientific community. In our opinion, there are no scientific arguments to support its marginal 
status, but only cultural reasons. One of the motivations for writing this book is helping in the maturing process that 
the scientific community needs about Bohmian mechanics. 
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Certainly, the distilling process of Bohmian mechanics is being quite slow. Anyone interested enough to walk this 
causal road of quantum mechanics can be easily confused by many misleading signposts that have been raised in the 
scientific literature, not only by its detractors, but unfortunately very often also by some of its advocates. Nowadays, 
following opinions from other reputed physicists (we are easily persuaded by those scientists with authority) is far 
from being a proper scientific strategy to get our own opinion about Bohmian mechanics. 

In any case, since the mathematical structure of Bohmian mechanics is quite simple, it can be easily learned by 
anyone with only a basic knowledge of classical and quantum mechanics who makes the necessary effort to build his 
own scientific opinion based on logical deductions, free from cultural influences. The introductory chapter of this 
book, including a thorough list of exercises and easily programmable algorithms, provides a reasonable and objective 
source of information in order to achieve this later goal, even for undergraduate students. 

Curiously, the fact that Bohmian mechanics is ignored and remains mainly unexplored is an attractive feature for 
some adventurous scientists. They know that very often new cutting edge ideas come from outside of the main stream 
and find in Bohmian mechanics a useful tool in their research activity. On the one hand, it provides an explanation 
of quantum mechanics, in terms of trajectories, that results to be very useful in explaining the dynamics of quantum 
systems, being thus also a source of inspiration to look for novel quantum phenomena. On the other hand, since 
it provides an alternative mathematical formulation, Bohmian mechanics offers new computational tools to explore 
physical scenarios that presently are computationally inaccessible, such as many-particle solutions of the Schrodinger 
equation. In addition, Bohmian mechanics sheds light on the limits and extensions of our present understanding of 
quantum mechanics towards other paradigms such as relativity or cosmology, where the internal structure of Bohmian 
mechanics in terms of well-defined trajectories is very attractive. With all these previous motivations in mind, this 
book provides eight chapters with practical examples showing how Bohmian mechanics helps us in our daily research 
activity. 

Obviously, there are other books focused on Bohmian mechanics. However, many of them are devoted to the 
foundations of quantum mechanics emphasizing the difficulties or limitations of the Copenhagen interpretation for 
providing an ontological description of our world. On the contrary, this book is not at all devoted to the foundations 
of quantum mechanics, but only to discuss about the practical application of the ideas of de Broglie and Bohm to 
understand the quantum world. Several examples of such practical applications written by leading experts in different 
fields, with an extensive updated bibliography, are provided here. The book, in general, is addressed to students in 
physics, chemistry, electrical engineering, applied mathematics, nanotechnology, as well as to both theoretical and 
experimental researchers who seek new computational and interpretative tools for their everyday research activity. 
We hope that the newcomers to this causal explanation of quantum mechanics will use Bohmian mechanics in their 
research activity so that Bohmian mechanics will become more and more popular for the broad scientific community. 
If so, we expect that, in the near future, Bohmian mechanics will be taught regularly at the Universities, not as the 
unique and revolutionary way of understanding quantum phenomena, but as an additional and useful interpretation 
of all quantum phenomena in terms of quantum trajectories. In fact, Bohmian mechanics has the ability of removing 
most of the mysteries of the Copenhagen interpretation and, somehow, simplifying quantum mechanics. We will be 
very glad if this book can contribute to shorten the time needed to achieve all these goals. 

Finally, we want to acknowledge many different people that has allowed us to embark into and successfully finish 
this book project. First of all, Alfonso Alarcon and Albert Benseny who became involved in the book project from 
the very beginning, as two additional editors. We also want to thank the rest of the authors of the book for accepting 
our invitation to participate in this project and writing their chapters according to the general spirit of the book. Due 
to page limitations, only eight examples of practical applications of Bohmian mechanics in forefront research activity 
are presented in this book. Therefore, we want to apologize to many other researchers who could have certainly been 
also included in the book. We also want to express our gratitude to Pan Stanford Publishing for accepting our book 
project and for their kind attention during the publishing process. 



II. INTRODUCTION 

The beginning of the twentieth century brought surprising non-classical phenomena. Max Planck's explanation 
of the black body radiation [T], the work of Albert Einstein on the photoelectric effect [5] and the Niels Bohr's 
model to account for the electron orbits around the nuclei [3H5] established what is now known as the old quantum 
theory. To describe and explain these effects, phenomenological models and theories were first developed, without 
any rigorous and global justification. In order to provide a complete explanation for the underlying physics of such 
new nonclassical phenomena, physicists were forced to abandon classical mechanics to develop novel, abstract and 
imaginative formalisms. 

In 1924, Louis de Broglie suggested in his doctoral thesis that matter, apart from its intrinsic particle-like behavior, 
could exhibit also a wave-like one [B^. Three years later he proposed an interpretation of quantum phenomena 
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based on non-classical trajectories guided by a wave field [7]. This was the origin of the pilot-wave formulation 
of quantum mechanics that we will refer as Bohmian mechanics to account for the following work of David Bohm 
[SI [S]- In the Bohmian formulation, an individual quantum system is formed by a point particle and a guiding 
wave. Contcmporanely, Max Born and Werner Heisenberg, in the course of their collaboration in Copenhagen with 
Niels Bohr, provided an original formulation of quantum mechanics without the need of trajectories [101 lllj . This 
was the origin of the so-called Copenhagen interpretation of quantum phenomena and, since it is the most accepted 
formulation, it is the only one explained at most universities. Thus, it is also known as the orthodox formulation of 
quantum mechanics. In the Copenhagen interpretation, an individual quantum system is simultaneously both a wave 
and a particle and exhibits its wave or its particle nature depending on the experimental arrangement. 

The present status of Bohmian mechanics among the scientific community is quite marginal (the quantum chemistry 
community is an encouraging exception). Most researchers do not know about it or believe that is not fully-correct. 
There are others that know that quantum phenomena can be interpreted in terms of trajectories, but they think that 
this formalism cannot be useful in their daily research activity. Finally, there are few researchers, the authors of this 
book among them, who think that Bohmian mechanics is a useful tool to make progress in front-line research fields 
involving quantum phenomena. 

The main (non-scientific) reason why still many researchers believe that there is something wrong with Bohmian 
mechanics can be illustrated with Hans Christian Andersen's tale The Emperor's new clothes. Two swindlers promise 
the Emperor the finest clothes that, as they tell him, are invisible to anyone who is unfit for their position. The 
Emperor cannot see the (non-existing) clothes, but pretends that he can for fear of appearing stupid. The rest of the 
people do the same. Advocates of the Copenhagen interpretation have attempted to produce impossibility proofs in 
order to demonstrate that Bohmian mechanics is incompatible with quantum phenomena [13^. Most researchers, who 
are not aware of the incorrectness of such proofs, might conclude that there is some controversy with the Bohmian 
formulation of quantum mechanics and they prefer not to support it, for fear of appearing discordant. At the end of 
the tale, during the course of a procession, a small child cries out 'Hhe Emperor is nakedP\ In the tale of quantum 
mechanics, David Bohm [51 [5] and John Bell [13] were the first to exclaim to the scientific community ^'Bohmian 
mechanics is a correct interpretation of quantum phenomena that exactly reproduce the predictions of the orthodox 
interpretationP' . 



How Bohmian mechanics helps 

Researchers who have spent the necessary time to analyze the ideas of de Broglie and Bohm with the pertinent 
scientific rigor accept that there is no scientific argument against them. Bohmian mechanics agrees with all quantum 
experiments done up to now. However, most of such researchers do also believe that Bohmian mechanics is not a 
useful tool to do research. In the words of Steven Weinberg, in a private exchange of letters with Sheldon Goldstein 
[15j : "In any case, the basic reason for not paying attention to the Bohm approach is not some sort of ideological 
rigidity, but much simpler — it is just that we are all too busy with our own work to spend time on something that 
doesn't seem likely to help us make progress with our real problems. ". 

The history of science seems to give credit to Weinberg's sentence. In spite of the controversies that have always 
been associated with the Copenhagen interpretation since its birth a century ago, its mathematical and computational 
machinery has enabled physicists, chemists and (quantum) engineers to calculate and predict the outcome of a vast 
number of experiments, while the contribution of Bohmian mechanics during the same period is much less significative. 
In our opinion, the differences are due to the fact that Bohmian mechanics remains mainly unexplored. 

Contrarily to Weinberg's opinion, we believe that Bohmian mechanics can help us make progress with our real 
problems. There are, at least, three clear reasons why one could be interested in studying quantum problems with 
Bohmian mechanics: 

1. Bohmian explaining: Even when the Copenhagen mathematical machinery is used to compute observable 
results, the Bohmian formulation ofently offers better interpretational tools. We can find descriptions of elec- 
tron dynamics such as "an electron crosses a resonant tunneling barrier and interacts with another electron 
inside the well". However, an electron crossing a tunneling region is not rigourously supported within orthodox 
quantum mechanics, but it is within the Bohmian picture. Thus, in contrast to the Copenhagen formulation, 
the Bohmian formulation allows for an easy visualization of quantum phenomena in terms of trajectories that 
has important demystifying or clarifying consequences. In particular, it provides a single-event description of 
the experiment, while Copenhagen quantum mechanics accounts for its statistical or ensemble explanation. We 
will present several examples in chapters 2 and 3 emphasizing all these points. 
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2. Bohmian computing: Although the predictions of the Bohmian formulation reproduce the ones of the or- 
thodox formulation of quantum mechanics, its mathematical formalism is different. In some systems, Bohmian 
equations might provide better computational tools than the ones obtained from the orthodox machinery, 
resulting in a reduction of the computational time, an increase in the number of degrees of freedom directly 
simulated, etc. We will see examples of these computational issues in quantum chemistry in chapters 4 and 5, 
as well as in quantum electron transport in Chap. 6. 

3. Bohmian thinking: From a more fundamental point of view, alternative formulations of quantum mechanics 
can provide alternative routes to look for the limits and possible extensions of the quantum theory. As we 
will discuss later, the work of John Bell on non-locality is a clear example of the unquestionable utility of 
understanding quantum phenomena with Bohmian mechanics. In particular. Chap. 7 presents the route to 
connect Bohmian mechanics with geometrical optics and beyond. The natural extension of Bohmian mechanics 
to the relativistic regime and to quantum field theory are presented in Chap. 8, while Chap. 9 discusses its 
application to cosmology. 

The fact that all measurable results of the orthodox quantum mechanics can be exactly reproduced with Bohmian 
mechanics is the relevant point that completely justifies why Bohmian mechanics can be used for explaining or 
computing different quantum phenomena in physics, chemistry, electrical engineering, applied mathematics, nan- 
otechnology, etc. In the scientific literature, the Bohmian computing technique to find the trajectories (without 
directly computing the wavefunction) is also known as a syntectic technique, while the Bohmian explaining tech- 
nique (where the wavefunction is directly computed first) is referred to the analytic technique ^16] . Furthermore, 
among others, the fact that Bohmian mechanics provides a picture in terms of single events which is missing in 
the orthodox explanation, is an attractive feature for those researchers interested in thinking about the limits or 
extensions of the quantum theory. 

In order to convince the reader about the practical utility of Bohmian mechanics for explaining, computing 
or thinking, we will not present elaborated mathematical developments or philosophical discussions, but provide 
practical examples. Apart from the first chapter, devoted to an overview of Bohmian mechanics, the book is divided 
into eight additional chapters with several examples on the practical application of Bohmian mechanics to different 
research fields, ranging from atomic systems to cosmology. These examples will clearly show that the previous 
quotation by Weinberg does not have to be always true. 

On the name Bohmian mechanics 

Any possible newcomer to Bohmian mechanics can certainly be quite confused and disoriented by the large list 
of names and slightly different explanations of the original ideas of de Broglie and Bohm that are present in the 
scientific literature. Different researchers use different names. Certainly, this is an indication that the theory is still 
not correctly settled down among the scientific community. 

In his original works [HI E] : de Broglie used the term pilot-wave theory [T7] , to emphasize the fact that wave- fields 
guided the motion of point particles. After de Broglie abandoned his theory, Bohm rediscovered it in the seminal 
papers entitled "A suggested interpretation of the quantum theory in terms of hidden variables" The term hidden 

variable^ refering to the positions of the particles, was perhaps pertinent in 1952, in the context of the impossibility 
proofs [13} . Nowadays, these words might seem inappropriate because they suggest something metaphysical on the 
trajectorie^ 

To give credit to both de Broglie and Bohm, some researchers refer to their works as the de Broglie-Bohm f/ieorj[^20J . 
Some reputed researchers argue that de Broglie and Bohm did not provide the same exact presentation of the theory 
[T71 m]. While de Broglie presented a first order development of the quantum trajectories (integrated from the 
velocity), Bohm himself did a second order (integrated from the acceleration) emphasizing the role of the quantum 
potential. The differences between both approaches appear when one considers initial ensembles of trajectories which 



^ Note that the term hidden variables can also refer to other (local and non-local) formulations of quantum mechanics. 

^ Sometimes it is argued that the name hidden variables is because Bohmian trajectories cannot be measured directly. However, what 
is not directly measured in experiments is the (complex) wavefunction amplitude, while the final positions of particles can be directly 
measured, for example, by the imprint they leave on a screen. John S. Bell wrote 14 (page 201): "Absurdly, such theories are known 
as 'hidden variable' theories. Absurdly, for there it is not in the wavefunction that one finds an image of the visible world, and the 
results of experiments, but in the complementary 'hidden' (!) variables. " 

^ In fact, even de Broglie and Bohm were not the original names of the scientists' families. Louis de Broglie's family, which included 
dukes, princes, ambassadors and marshals of France, changed their original Italian name Broglia to de Broglie when they established in 
France in the seventeenth century 1181 . David Bohm's father, Shmuel Diim, was born in the Hungarian town of Munkacs and, was sent 
to America when he was young. Upon landing at Ellis Island, he was told by an immigration official that his name, Diim would mean 
"stupid" in English. The official himself decided to change the name to Bohm |19l . 
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FIG. 1. Number of citations per year for (a) the two 1952 Bohm's papers entitled "A suggested interpretation of the quantum 
theory in terms of hidden variables" [S] [5] and (b) de Broglie's paper "La mecanique ondulatorie et la structure atomique de 
la matiere et du rayonnement" [7]. Data retreived from ISI Web of Knowledge |26] in December 2010. 



are not in quantum equilibriun:j2j Except for this issue, which will not be addressed in this book, both approaches 
are identical. 

Many researchers prefer to use the name Bohmian mechanics |22j. It is perhaps the most popular name. We 
know directly from his alive collaborators, Basile Hiley and David Peat [53], that this name irritated David Bohm 
and he said about its own work "if's Bohmian non-mechanics^\ He argued that the 'quantum potential' is a non- 
local potential that depends on the relative shape of the wavefunction and thus it is completely different from other 
mechanical (such as the gravitational or the electrostatic) potentials. See this particular discussion in the last chapters 
of Bohm and Hiley's book entitled The Undivided Universe: An Ontological Interpretation of Quantum Theory [2i] . 
He preferred the names causal or ontological interpretation of quantum mechanics (20ti24) . The latter names emphasize 
the foundational aspects of its formulation of quantum mechanics. 

Finally, another very common term is quantum hydrodynamics [16 that underlines the fact that Bohmian trajec- 
tories provide a mathematical relationship between the Schrodinger equation and fluid dynamics. In fact, this name 
is more appropriate when one refers to the Madelung's theory [25 , which is considered as a precursor of de Broglie's 
and Bohm's work, see Sec. |III.8[ 

From all these different names, we chose Bohmian mechanics because it is short and clearly specifies what we 
are referring to. It has the inconvenience of not giving credit to the initial work of Louis de Broglie. Although it 
might be argued that Bohm merely reinterpreted the prior work of de Broglie, we think that he was the first person to 
genuinely understand its significance and implications. As we mentioned, Bohm himself disliked this name. However, 
as any work of art, the merits and influence of Bohm's 1952 paper evolve independently of the own author opinion. 
f\ The paper becomes part of the scientific heritage. In any case, we understand Bohmian mechanics as a generic 
name that includes all those works inspired from the original ideas of Bohm and de Broglie. In the Fig. [TJ we plot the 
numbers of citations per year for the 1952 Bohm's seminal papers [SJ E] , certifying the exponentially growing influence 
of these papers, which is not the case for the original work of de Broglie [7]. 



* Quantum equilibrium assumes that the initial positions and velocities of Bohmian trajectories are defined compatible with the initial 
wavefunction. Then the trajectories computed from Bohm's or de Broglie's formulations will become identical. However, if the quantum 
equilibrium hypothesis is not considered, one can select completely arbitrary initial positions from the (first order) de Broglie explanation 
and arbitrary initial velocities and positions from the (second order) Bohm work, see Sec. |IV.6[ 

^ For example, Erwin Schrodinger, talking about quantum theory, wrote: "/ don't like it, and 1 m sorry I ever had anything to do with 
it", but his opinion did not influence the great applicability of his famous equation. 
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On the book contents 

The book contains nine chapters. It can be clearly divided in two different parts. The first part is formed by chapter 
1 alone that provides an accessible introduction to Bohmian mechanics. The second part is composed of the rest of 
chapters with practical examples on the applicability of Bohmian mechanics. 

Chapter 1 is the longest one and it is entitled Overview of Bohmian mechanics. It is written by Xavier 
Oriols and Jordi Mompart, the editors of the book, both from the Univcrsitat Autonoma de Barcelona, Spain. This 
chapter is intended to be an introduction to any newcomer interested in Bohmian mechanics. No previous knowledge 
on Bohmian mechanics is required. Only basic concepts of classical and quantum mechanics is assumed. The chapter 
is divided into four different sections. First, the historical development of Bohmian mechanics is presented. Then, 
Bohmian mechanics for single particle and for many particle systems (with spin and entanglement explanations) are 
discussed in sections two and three, respectively. The topic of Bohmian measurement is addressed in sc;ction four. 
The chapter contains also a list of exercises, mathematical demonstrations an easily implementable algorithms for 
computation of Bohmian trajectories. 

Chapter 2 is entitled Hydrogen photoionization with strong lasers. It is written by Albert Benseny and 
Jordi Mompart from the Universitat Autonoma de Barcelona, Spain; Antonio Picon from ANL, the Argonne National 
Laboratory, Illinois, USA; Luis Plaja from the Universidad de Salamanca, Spain; and Luis Roso from the CLPU, 
the Laser Center for Ultrashort and Ultraintense Pulses, also in Salamanca. They discuss the dynamics of a single 
hydrogen atom interacting with a strong laser. In particular, the Bohmian trajectories of these electrons represent 
an interesting illustrating view of both the above threshold ionization and the harmonic generation spectra problems. 
They do also present a full three dimensional model to discuss the dynamics of Bohmian trajectories when the light 
beam and the hydrogen atom exchange both spin and orbital angular momentum. The chapter does also provide 
a practical example on how Bohmian mechanics are computed, with an analytical (i.e., Bohmian explaining) 
procedure, when full (scalar and vector potentials) electromagnetic fields are considered. 

The title of chapter 3 is Atomtronics: Coherent control of atomic flow via adiabatic passage. It is written 
by Albert Benseny, Joan Baguda, Xavier Oriols, and Jordi Mompart from the Universitat Autonoma de Barcelona, 
Spain; and Gerhard Birkl from the Institut fiir Angewandte Physik, from the Technische Universitat Darmstadt, 
Germany. Here, it is discussed an efficient and robust technique to coherently transport a single neutral atom, a single 
hole, or even a Bosc-Einstein condensate between the two extreme traps of the triple-well potential. The dynamical 
evolution of this system with the direct integration of the Schrodinger equation presents a very counterintuitive 
effect: by slowing down the total time duration of the transport process it is possible to achieve atomic transport 
between the two extrem traps with a very small (almost negligible) probability to populate the middle trap. The 
analytical (i.e., Bohmian explaining) solution of this problem with Bohmian trajectories enlightens the role of the 
particle conservation law in quantum systems showing that the negligible particle presence is due to a sudden particle 
acceleration yielding, in fact, ultra- high atomic velocities. The Bohmian contribution opens the discussion about the 
possible detection of such high kinetic energies or the need for a relativistie formulation to accurately describe such 
a simple quantum system. 

Chapter 4, entitled The role of trajectories in quantum chemistry and chemical physics, is written by 
Angel S. Sanz and Salvador Miret-Artes from the Instituto de Fisica Fundamental, Consejo Superior de Investigaciones 
Cientfficas, Spain. Up to very recently, the Bohmian computing abilities have been explored almost exclusively by 
the quantum chemistry community. In the own words of the authors: "It is interesting to note that Bohmian mechanics 
is better accepted within the chemistry community than within the physics one, which is often quite reluctant. . . To some 
extent, it could be said that the pedagogical advantages that Bell associated with Bohmian Mechanics have been better 
appreciated by chemists than by physicists. ". There are many different Bohmian computational techniques developed 
by the quantum chemistry community, mainly, during the last ten years. Since it was impossible to dedicate a chapter 
to each, here, the authors provide a clear summary of the problems and Bohmian solutions. 

Chapter 5, whose title is Adaptive quantum Monte Carlo approach states for high-dimensional systems, 
is written by Eric R. Bittner, Donald J. Kouri, Sean Derrickson, from University of Houston; and Jeremy B. Maddox, 
from the Western Kentucky University, in USA. They provide one particular example on the success of Bohmian 
mechanics in the chemistry community. In this chapter, the authors explain their Bohmian computing development 
to calculate the ah initio qiiantum mechanical structure, energetics and thermodynamics of multi-atoms systems. They 
use a variational approach that finds the quantum ground sate (or even excited states at finite temperature) using 
a statistical modeling approach for determining the best estimate of a quantum potential for a multi-dimensional 
system. 

Chapter 6 is entitled Nanoelectronics: Quantum electron transport. It is written by Alfonso Alarcon, 
Guillem Albareda, Fabio Lorenzo Traversa and Xavier Oriols from the Universitat Autonoma de Barcelona, Spain. 
The authors explain the ability of their own many-particle Bohmian computing algorithm to understand and model 
the behavior of nanoscale electron devices. In particular, it is shown that the adaptation of Bohmian mechanics to 
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electron transport in open systems (with interchange of particles and energies) leads to a quantum Monte Carlo 
algorithm, where randomness appears because of the uncertainties in the number of electrons, their energies and the 
initial positions of (Bohmian) trajectories. A general, versatile and time-dependent 3D electron transport simulator for 
nanoelectronic devices, named BITLLES (Bohmian Interacting Transport for nonequiLibrium eLEctronic Structures), 
is presented showing its ability for a full prediction (DC, AC, fluctuations) of the electrical characteristics of any 
nanoelectronic device. 

Chapter 7, entitled Beyond the eikonal approximation in classical optics and quantum physics, is 

written by Adriano Orefice, Raffaele Giovanelli and Domenico Ditto from the Universita degli Studi di Milano, Italy. 
It is devoted to discuss how Bohmian thinking can also help in optics, exploring the fact that the time-independent 
Schrodinger equation is strictly analogous to the Helmholtz equation appearing in classical wave theory. Starting 
from this equation they obtain indeed, without any omission or approximation, a Hamiltonian set of ray-tracing 
equations providing (in stationary media) the exact description in term of rays of a family of wave phenomena (such 
as diffraction and interference) much wider than that allowed by standard geometrical optics, which is contained as a 
simple limiting case. They show in particular that classical ray trajectories are ruled by a wave potential presenting 
the same mathematical structure and physical role of Bohm's quantum potential, and that the same equations of 
motion obtained for classical rays hold, in suitable dimensionless form, for quantum particle dynamics, leading to 
analogous trajectories and reducing to classical dynamics in the absence of such a potential. 

Chapter 8, entitled Relativistic quantum mechanics and quantum field theory, is written by Hrvoje Nikolic 
from the Rudjer Boskovic Institute, Croatia. This chapter presents a clear example on how a Bohmian thinking on 
superluminal velocities and nonlocal interactions helps in extending the quantum theory towards relativity and quan- 
tum field theory. A relativistic-covariant formulation of relativistic quantum mechanics of fixed number of particles 
(with or without spin) is presented, based on many-time wavefunctions and on an interpretation of probabilities in 
the spacetime. These results are used to formulate the Bohmian formulation of relativistic quantum mechanics in a 
manifestly relativistic covariant form and are also generalized to quantum field theory. The corresponding Bohmian 
formulation of quantum field theory describes an infinite number of particle trajectories. Even though the particle 
trajectories are continuous, the appearance of creation and destruction of a finite number of particles results from 
quantum theory of measurements describing entanglement with particle detectors. 

Finally, chapter 9, whose title is Sub-quantum accelerating universe, is written by Pedro F. Gonzalez- 
Dfaz and Alberto Rozas-Fernandez from the Instituto de Ffsica Fundamental, Consejo Superior de Investigaciones 
Cientfficas, Spain. Contrarily to the general belief, quantum mechanics does not only govern microscopic systems, 
but it has influence also on the cosmological domain. However, the extension of the Copenhagen version of quantum 
mechanics to cosmology is not free from conceptual difficulties: the probabilistic interpretation of the wavefunction of 
the whole universe is somehow misleading because we cannot make statistical 'measurements' of different realizations 
of our universe. This chapter deals with two new cosmological models describing the accelerating universe in the 
spatially flat case. Also contained in this chapter is a discussion on the quantum cosmic models that result from the 
existence of a nonzero entropy of entanglement. In such a realm, they obtain new cosmic solutions for any arbitrary 
number of spatial dimensions, studying the stability of these solutions, as well as the emergence of gravitational waves 
in the realm of the most general models. 

Only eight chapters of practical applications of Bohmian mechanics in forefront research activity, ranging from 
atomic systems to cosmology, are presented in this book. Many more chapters could have certainly been also included. 
In any case, we believe that these examples are enough to convince the reader about the practical utility of explaining, 
computing and thinking about quantum phenomena within Bohmian mechanics. 



III. HISTORICAL DEVELOPMENT OF BOHMIAN MECHANICS 

In general, the history of quantum mechanics is explained in textbooks as a chronicle where each step follows 
naturally from the one preceding it. However, it was exactly the opposite. The development of quantum mechanics 
was a zigzagging route full of personal disputes. It was a painful history, where scientists were forced to abandon 
well-established classical concepts and to explore new and imaginative routes. Most of the new routes went nowhere. 
Others were simply abandoned. Some of the explored routes were successful in providing new mathematical formalisms 
capable of predicting experiments at the atomic scale. Even such successful routes were painful enough, so relevant 
scientists, such as Albert Einstein or Erwin Schrodinger, decided not to support them. In this section we will briefly 
explain the history of one of these routes: Bohmian mechanics. It was first proposed by Louis de Broglie [27], who 
abandoned it soon afterward, and rediscovered by David Bohm Oil] many years later, and it has been almost ignored 
by the scientific community since then. Here, we will discuss the historical development of Bohmian mechanics to 
understand its present status. Also, we will introduce the basic mathematical aspects of the theory, while the formal 



11 



and rigorous structure will be presented in the subsequent sections. 

III.l. Particles and waves 

The history of physics is a difficult topic because it does not only deals with its proper history but also deals with 
physics. In particular, quantum theory revolves around the notions of particles and waves. In classical physics, the 
concept of a particle is very useful for the description of many natural phenomena. A particle is directly related to 
a trajeciorj^l^ [t] that defines its position as a continuous function of time, usually found as a solution of a set of 
differential equations. For example, the planets can be considered particles orbiting around the sun, whose orbits are 
determined by the classical Newton gravitational laws. 

In classical mechanics, it is natural to think that the total number of particles (e.g., planets in the solar system) 
is conserved, and the particle trajectories must be continuous in time: if a particle goes from one place to another, 
then, it has to go through all the trajectory positions between these two places. This condition can be summarized 
with a local conservation law: 

^ + vMt) = o (1) 

where p{r,t) is the density of particles and j{r,t) is the particle current density. For an ensemble of point particles 

at positions fi[t] with velocites Vi[t], it follows that /9(r, t) = J2H^~ ^iM) ^^'^ j(^7^) — J2^i[^]^i^ ^ ^iM): '^ith S{r) 
being the Dirac delta function. 

However, the total number of planets in the solar system could be conserved in another (quite different) way. A 
phenomenon where a planet disappearing (instantaneously) from its orbit and appearing (instantaneously) at another 
point far away from its original location would certainly conserve the number of planets but would violate Eq. ([I]) . 
We must then think of Eq. ([T]) as a local particle number conservation law. 

Fields, and particularly waves, also appear in many explanations of physical phenomena. The concept of a field 
was initially introduced to deal with the interaction of distant particles. For example, there is an interaction between 
the electrons in an emitting radio antenna at the top of a mountain and those in the receiving antenna at home. Such 
interaction can be explained through the use of an electromagnetic field. Electrons in the transmitter generate an 
electromagnetic field, a radio wave, that propagates through the atmosphere and arrives at our antenna, affecting its 
electrons. Finally, a loudspeaker transforms the electron motion into music at home. 

The simplest example of a wave is the so-called plane wave: 

^(r,i) = e'M-^-*") (2) 

where the angular frequency w and the wave vector k refer respectively to its temporal and spatial behavior. In 
particular, the angular frequency w specifies when the temporal behavior of such wave is repeated. The value of 
ip{'fl,ti) at position fl and time ti is identical to '4'ii^i,t2) when t2 = ti + 2Trn/u! for n integer. The angular frequency 
to can be related to the linear frequency as lo = lixv. Analogously, the wave vector A; determines the spatial repetition 
of the wave, that is, the wavelength A. The value of il}(f\^t\) at position f\ and time t\ is identical to '^{fi^t\) when 
k ■ f2 = k ■ f\ \- 2TTn with n integer. Unlike a trajectory, a wave is defined at all possible positions and times. Waves 
can be a scalar or a vectorial function and take real or complex values. For example, Eq. Q is a scalar complex wave 
of unit amplitude. The shape of waves is determined by a set of differential equations. In our broadcasting example. 
Maxwell equations define the electromagnetic field of the emitted radio wave that is given by two vectorial functions, 
one for the electric field and one for the magnetic field. 

Whenever the differential equations that govern the fields are linear, one can apply the superposition principle to 
explain what happens when two or more fields (waves) traverse simultaneously the same region. The modulus of 
the total field at each position is related to the amplitudes of the individual waves. In some cases, the modulus of 
the sum of the amplitudes is much smaller than the sum of the modulus of the amplitudes; this is called destructive 
interference. In other cases, it is roughly equal to the sum of the modulus of the amplitudes; this is called constructive 
interference. 



In order to avoid confusion, let us emphasize that in the orthodox formulation of quantum mechanics, the concept of a particle is not 
directly related to the concept of a trajectory. For example, the electron is a particle but it has no precise trajectory, as we will discuss 
later. 
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III. 2. Origins of the quantum theory 

At the end of the nineteenth century, Sir Joseph John Thomson discovered the electron, and in 1911, Ernest 
Rutherford, a New Zealander student working in Thomson's laboratory, provided experimental evidence that inside 
atoms, electrons orbited around a nucleus in a similar manner as planets do around the sun. Rutherford's model of 
the atom was clearly in contradiction with well-established theories, since classical electromagnet ism predicted that 
orbiting electrons should radiate, gradually lose energy, and spiral inward. Something was missing in the theory, 
since it seemed that the electron behavior inside an atom could not be explained in terms of trajectories. Therefore, 
alternative ideas needed to be explored to understand atom stability. 

In addition, at that time, classical electromagnetism was unable to explain the radiated spectrum of a black body, 
which is an idealized object that emits a temperature-dependent spectrum of light (like a big fire with different 
colors, depending on the flame temperature). The predicted continuous intensity spectrum of this radiation became 
unlimitcdly large in the limit of large frequencies, resulting in an unrealistic emission of infinite power, which was 
called the ultraviolet catastrophe. However, the measured radiation of a black body did not behave in this way, 
indicating that a wave description of the electromagnetic field was also incomplete. 

In summary, at the beginning of the twentieth century, it was clear that natural phenomena such as atom stability 
or black-body radiation, were not well explained in terms of a particle or a wave description alone. It seemed necessary 
to merge both concepts. 

In 1900, Max Planck suggested pQ that black bodies emit and absorb electromagnetic radiation in discrete energies 
hv, where v is the frequency of the emitted radiation and h is the (now-called) Planck constant. Five years later, 
Einstein used this discovery in his explanation of the photoelectric effect 2] , suggesting that light itself was composed 
of light quanta or photona^oi energy hiy. Even though this theory solved the black-body radiation problem, the fact 
that the absorption and emission of light by atoms are discontinuous was still in confiict with the classical description 
of the light-matter interaction. 

In 1913, Niels Bohr [3-5] wrote a revolutionary paper on the hydrogen atom, where he solved the (erroneously 
predicted in classical terms) instability by postulating that electrons can only orbit around atoms in some particular 
nonradiating orbits. Thus, atom radiation occurs when electrons jump from one orbit to another of lower energy. 
His imaginative postulates were in full agreement with the experiments on spectral lines. Later, in 1924, de Broglie 
proposed, in his PhD dissertation, that all particles (such as electrons) exhibit wave-like phenomena like interference 
or diffraction |27j . In particular, one way to arrive at Bohr's hypothesis is to think that the electron orbiting around 
the proton is a stationary wave that imposes simple restrictions on the shape of such a stationary wave. Physics at 
the atomic scale started to be understandable by mixing the concepts of particles and waves. All these advances were 
later known as the old quantum theory. The word quantum referred to the minimum unit of any physical entity (e.g., 
the energy) involved in an interaction. 

III. 3. "Wave or particle?" vs. "wave and particle" 

In the mid-1920s, theoreticians found themselves in a difficult situation when attempting to advance Bohr's ideas. 
A group of atomic theoreticians centered on Bohr, Max Born, Wolfgang Pauli, and Werner Heisenberg suspected that 
the problem went back to trying to understand electron trajectories within atoms. In under two years, a series of 
unexpected discoveries brought about a scientific revolution |12| . 

Heisenberg wrote his first paper on quantum mechanics in 1925 |11| and two years later stated his uncertainty 
principle It was him, with the help of Born and Pascual Jordan, who developed the first version of quantum 

mechanics based on a matrix formulation fTHl HH [30l |3T] . 

In 1926, Schrodinger published An Undulatory Theory of the Mechanics of Atoms and Molecules '32], where, inspired 
by de Broglie's work tSlEllSJi described material points (such as electrons or protons) in terms of a wave solution 
of the following (wave) equation: 

^^dM^ ^ _|lvV(r,t) + V^(r,t)^(f,i) (3) 

where V{r,t) is the potential energy felt by the electron, and the wave (field) ip{f,t) was called the wave function. 
Schrodinger, at first, interpreted his wave function as a description of the electron charge density q\Tlj{f,t)\ with q 



In fact, the word photon was not coined until 1926, by Gilbert Lewis 1281 . 
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the electron charge. Later, Born refined the interpretation of Schrodinger and defined t)p as the probabifity 
density of finding the electron in a particular position r at time t |12j . 

Schrodinger's wave version of quantum mechanics and Heisenberg's matrix mechanics were apparently incompatible, 
but they were eventually shown to be equivalent by Wolfgang Ernst Pauli and Carl Eckart, independently ^TH I33j . 

In order to explain the physics behind quantum systems, the concepts of waves and particles should be merged in 
some way. Two different routes appeared: 

1. Wave or particle? The concept of a trajectory was, consciously or unconsciously, abandoned by most of 
the young scientists (Heisenberg, Pauli, Dirac, Jordan,...). They started a new route, the wave or particle? 
route, where depending on the experimental situation, one has to choose between a wave or a particle behavior. 
Electrons are associated basically to probability (amplitude) waves. The particle nature of the electron appears 
when we measure the position of the electron. In Bohr's words, an object cannot be both a wave and a particle 
at the same time; it must be either one or the other, depending upon the situation. This approach, mainly 
supported by Bohr, is one of the pillars of the Copenhagen, or orthodox, interpretation of quantum mechanics. 

2. Wave and particle: de Broglie, on the other hand, presented an explanation of quantum phenomena where 
the wave and particle concepts merge at the atomic scale, by assuming that a pilot-wave solution of Eq. ^ 
guides the electron trajectory. This is what we call the Bohmian route. One object cannot be a wave and a 
particle at the same time, but two can. 

The differences between the two routes can be easily seen in the interpretation of the double-slit experiment. A 
beam of electrons with low intensity (so that electrons are injected one by one) impinges upon an opaque surface 
with two slits removed on it. A detector screen, on the other side of the surface, detects the position of electrons. 
Even though the detector screen responds to particles, the pattern of detected particles shows the interference fringes 
characteristic of waves. The system exhibits, thus, the behavior of both waves (interference patterns) and particles 
(dots on the screen). 

According to the wave or particle? route, first the electron presents a wavelike nature alone when the wave 
function (whose modulus gives the probability density of finding a particle when a position measurement is done) 
travels through both slits. Suddenly, the wave function collapses into a delta function at a (random) particular 
position on the screen. The particle-like nature of the electron appears, while its wavelike nature disappears. Since 
the screen positions where collapses occur follow the probability distribution dictated by the modulus of the wave 
function, a wave interference pattern appears on the detector screen. 

According to the wave and particle route, the wave function (whose modulus means the particle probability density 
of being at a certain position, regardless of the measurement process) travels through both slits. At the same time, 
a well-defined trajectory is associated with the electron. Such a trajectory passes through only one of the slits. The 
final position of the particle on the detector screen and the slit through which the particle passes is determined by 
the initial position of the particle. Such an initial position is not controllable by the experimentalist, so there is an 
appearance of randomness in the pattern of detection. The wave function guides the particles in such a way that they 
avoid those regions in which the interference is destructive and are attracted to the regions in which the interference 
is constructive, giving rise to the interference pattern on the detector screen. Let us quote the enlightening summary 
of Bell [Hj: 

Is it not clear from the smallness of the scintillation on the screen that we have to do with a particle? And 
is it not clear, from the diffraction and interference patterns, that the motion of the particle is directed by 
a wave? De Broglie showed in detail how the motion of a particle, passing through just one of two holes in 
screen, could be influenced by waves propagating through both holes. And so influenced that the particle 
does not go where the waves cancel out, but is attracted to where they cooperate. This idea seems to me 
so natural and simple, to resolve the wave-particle dilemma in such a clear and ordinary way, that it is a 
great mystery to me that it was so generally ignored. 

Now, with almost a century of perspective and the knowledge that both routes give exactly the same experimental 
predictions, it seems that such great scientists took the strangest route. In fact. Bell argued that their choice was 
unprofessional [13]. Imagine that a student asks his or her professor, "What is an electron?" The answer of a 
(Copenhagen) professor could be, "The electron is not a wave nor a particle. But, do not worry! You do not have to 
know what an electron is to (compute observable results) pass the exam.'j^ If the student insists, the professor might 
reply, "Shut up and calculate."]^ 



For example, in the book Quantum Theory, 1341 written by Bohm before he formulated Bohmian mechanics in 1952, he wrote, when 
talking about the wave-particle duality; "We find a strong analogy here to the fable of the seven blind men who ran into an elephant: 
One man felt the trunk and said that 'an elephant is a rope'; another felt the leg and said that 'an elephant is obviously a tree,' and so 
on." 

This quote is sometimes attributed to Dirac, Richard Feynman, or David Mermin 1351 1361 . It recognizes that the important content of 
the orthodox formulation of quantum theory is the ability to apply mathematical models to real experiments. 
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Another example of the vagueness of the orthodox formulation can be illustrated by the question that Einstein 
posed to Abraham Pais: "Do you really think the moon is not there if you are not looking at it?" The answer of a 
Copenhagen professor, such as Bohr, would be, "I do not need to answer such a question, because you cannot ask me 
such question experimentally." 

On the other hand, an alternative (Bohmian) professor would answer, "Electrons are particles whose trajectories 
are guided by a pilot field which is the wave function solution of the Schrodinger equation. There is some uncertainty 
in the initial conditions of the trajectories, so that experiments have also some uncertainties." With such a simple 
explanation, the student would understand perfectly the role of the wave and the particle in the description of quantum 
phenomena. Furthermore, in the Bohmian interpretation, the position of an electron (or the moon) while we are not 
looking at it, is always defined, even though it is a hidden variable for experimentalists. 

One of the reasons that led the proponents of quantum mechanics to choose the wave or particle? route is that 
the predictions about the positions of electrons are uncertain because the wave function is spread out over a volume. 
This effect is known as the uncertainty principle: it is not possible to measure, simultaneously, the exact position 
and velocity (momentum) of a particle. Therefore, scientists preferred to look for an explanation of quantum effects 
without the concept of a trajectory that seemed unmeasurable. However, the argumentation about the impossibility 
of dealing with trajectories is, somehow, unfair and unjustified, since it relies on the "principle" that physical theories 
should not contain entities that cannot be observedF^I 



In addition, everyone with experience on Fourier transforms of conjugate variables recognizes the quantum uncer- 
tainty principle as a trivial effect present in any wave theory where the momentum of a particle depends on the slope 
of its associated wave function. Then, a very localized particle would have a very sharp wave function. In this case 
such a wave function would have a great slope that implies a large range of possible momenta. On the contrary, if 
the wave function is built from a quite small range of momenta, then it will have a large spatial dispersion. 



Perhaps the most relevant event for the development of the quantum theory was the fifth Solvay Conference, which 
took place from October 24-29, 1927, in Brussels [18]. As on previous occasions, the participants stayed at Hotel 
Britannique, invited by Ernest Solvay, a Belgian chemist and industrialist with philanthropic purposes due to the 
exploitation of his numerous patents. There, de Broglie presented his recently developed pilot wave theory and how 
it could account for quantum interference phenomena with electrons [H] . He did not receive an enthusiastic reaction 
from the illustrious audience gathered for the occasion. In the following months, it seems that he had some difficulties 
in interpreting quantum measurement with his theory and decided to avoid it. In fact, one (nonscientific) reason that 
forced de Broglie to give up on his theory was that he worked isolated, having little contact with the main research 
centers in Berlin, Copenhagen, Cambridge, or Munich. By contrast, most of the Copenhagen contributors worked 
with fluid and constant collaborations among them. 

Finally, let us mention that the elements of the pilot wave theory (electrons guided by waves) were already in place 
in de Broglie's thesis in 1924 [27], before either matrix or wave mechanics existed. In fact, Schrodinger used the de 
Broglie phases to develop his famous equation (see Eq. ([s])). In addition, it is important to remark that de Broglie 
himself developed a single-particle and a many-particle description of his pilot waves, visualizing also the nonlocality 
of the latter [18^. Perhaps, his remarkable contribution and influence have not been fairly recognized by scientists and 
historians because he abandoned his own ideas rapidly without properly defending them |18l 137] . 



Not even Einstein gave explicit support to the pilot wave theory |12j . It remains almost unknown that in 1927, 
the same year that de Broglie published his pilot wave theory |?J, Einstein worked out an alternative version of the 
pilot wave with trajectories determined by many-particle wave functions. However, before the paper appeared in 
print, Einstein phoned the editor to withdraw it. The paper remains unpublished, but its contents are known from a 
manuscript [33] . 

It seems that Einstein, who was unsatisfied with the Copenhagen approach, did not like the pilot wave approach 
either because both interpretations have this notion of action at a distance: particles that are far away from each 



From a philosophical point of view, this is known as "positivism" or "empiricism" and represents a nonphysical limitation on the possible 
kinds of theories that we could choose to explain quantum phenomena. For example, the wave function cannot be measured directly 
in a single experiment but only its modulus from an ensemble of experiments. However, there is no doubt that the (complex) wave 
function comes to be a very useful concept to understand quantum phenomena. Identically, in the de Broglie and Bohm interpretation, 
the trajectories cannot be measured directly, but they can also be a very interesting tool for understanding quantum phenomena. 




III. 4. Louis de Broglie and the fifth Solvay Conference 



III. 5. Albert Einstein and locality 
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other can profoundly and instantaneously affect each other. As the father of the theory of relativity, he believed that 
action at a distance cannot travel faster than the speed of light. Let Bohm explain the difficulties of Einstein with 
both the Bohmian and the orthodox interpretations |40) : 

In the fifties, I sent [my Quantum Theory book] around to various quantum physicists - including Niels 
Bohr, Albert Einstein, and Wolfgang Pauli. Bohr didn't answer, but Pauli liked it. Albert Einstein sent 
me a message that he'd like to talk with me. When we met he said the book had done about as well 
as you could do with quantum mechanics. But he was still not convinced it was a satisfactory theory. 
Einstein's objection was not merely that it was statistical. He felt it was a kind of abstraction; quantum 
mechanics got correct results but left out much that would have made it intelligible. I came up with the 
causal interpretation (that the electron is a particle, but it also has a field around it. The particle is never 
separated from that field, and the field affects the movement of the particle in certain ways). Einstein 
didn't like it, though, because the interpretation had this notion of action at a distance: Things that are 
far away from each other profoundly affect each other. He believed only in local action. 

Einstein, together with Boris Podolsky and Nathan Rosen, presented objections to quantum theory in the famous 
EPR article in 1935, entitled "Can Quantum-Mechanical Description of Physical Reality Be Considered Complete?" 
|41j . There, they argued that on the basis of the absence of action at a distance, quantum theory must be incomplete. 
In other words, quantum theory is either nonlocal or incomplete. Einstein believed that locality was a fundamental 
principle of physics, so he adhered to the view that quantum theory was incomplete. Einstein died in 1955, convinced 
that a correct reformulation of quantum theory would preserve local causality. 



III. 6. David Bohm and why the "impossibility proofs" were wrong? 

Perhaps the first utility of Bohm's work was the demonstration that von Neumann theorem had limited validity. 
In 1932, von Neumann put quantum theory on a firm theoretical basis [13^. Some of the earlier works lacked mathe- 
matical rigor, and he put the entire theory into the setting of operator algebra. In particular, von Neumann studied 
the following question: "If the present mathematical formulation of the quantum theory and its usual probability 
interpretation are assumed to lead to absolutely correct results for every experiment that can ever be done, can 
quantum-mechanical probabilities be understood in terms of any conceivable distribution over hidden parameters?" 
von Neumann answered this question negatively. His conclusions, however, relied on the fact that he implicitly re- 
stricted his proof to an excessively narrow class of hidden variables, excluding Bohm's hidden variables model. In 
other words, Bohmian mechanics is a counterexample that disproves von Neumann's conclusions, in the sense that it 
is possible to obtain the predictions of orthodox quantum mechanics with a hidden variables theory [T31 [20] ■ 

Bohm's formulation of quantum mechanic^^ appeared after the orthodox formalism was fully established. Bohm 
was, perhaps, the first person to genuinely understand the significance and fundamental implications of the description 
of quantum phenomena with trajectories guided by waves. Ironically, in 1951, Bohm wrote a book, Quantum Theory 
|34j . where he provided "proof that quantum theory is inconsistent with hidden variables" (see page 622 in [21]). In 
fact, he wrote in a footnote in that section, "We do not wish to imply here that anyone has ever produced a concrete 
and successful example of such a [hidden variables] theory, but only state that such theory is, as far as we know, 
conceivable." Furthermore, the book does also contain an unusually long chapter devoted to the quantum theory of the 
process of measurement, where Bohm discusses how the measurement itself can be described from the time evolution 
of a wave function, rather than invoking the wave function collapse. It seems that Bohm became dissatisfied with 
the orthodox approach that he had written in his book and began to develop his own causal formulation of quantum 
theory, which he published in 1952 ^ [5j . 

The original papers by Bohm [3 H] provide a formal justification of the guidance equation developed 25 years 
before by de Broglie. Instead of reproducing his exact mathematical development in terms of the quantum Hamilton- 
Jacobi equation, here, we discuss a very simple explanation of the guidance equation. By a simple mathematical 
manipulation of the (wave) equation, Eq. ([3|, we can find the local (particle) continuity equation, Eq. ([T]), discussed 
at the beginning of this sectionFj From the standard definition of the current density j{f,t), as a product of the 



Apart from these works, the history of science has recognized many other relevant contributions by Bohm I19| . As a postgraduate at 
Berkeley, he developed a theory of plasmas, discovering the electron phenomenon now known as Bohm diffusion. In 1959, with his 
student Yakir Aharonov, he discovered the Aharonov-Bohm effect, showing how a magnetic field could affect a region of space in which 
the field had been shielded, although its vector potential did not vanish there. This showed for the first time that the magnetic vector 
potential, hitherto a mathematical co nvenien ce, could have real physical (quantum) effects. 
See the formal demonstration in Sec. IV. 3. 2 or in Ref. 1421. 
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particle density p(f,t) = \%lj{r,t)\'^ and the velocity v{r,t), we can exactly obtain the guidance equation that was 
predicted by de Broglie and Bohm for the particle velocity: 

dt ' p{r,t) ^ ' 

There was another important point explained by Bohm in 1952. If one considers an ensemble of trajectories whose 
initial positions at time t = are distributed according to p{r,0) = \ip{r,0)\'^ , such an ensemble of trajectories will 
reproduce p{r, t) = t)P at any other time if the trajectories follow the guidance equation, Eq. Q. Therefore, we 
are able to exactly reproduce the time evolution of the wave function solution of Eq. ([s]) from an ensemble of trajec- 
tories guided by waves. Let Bohm himself explain this revolutionary point in the first page of his original reference [5]: 

The usual [Copenhagen] interpretation of the quantum theory is based on an assumption [. . . ] that the 
physical state of an individual system is completely specified by a wave function that determines only the 
probabilities of actual results that can be obtained in a statistical ensemble of similar experiments. [. . . ] 
In contrast, this alternative interpretation permits us to conceive of each individual system as being in 
a precisely definable state, whose changes with time are determined by definite laws, analogous to (but 
not identical with) the classical equations of motion. Quantum mechanical probabilities are regarded (like 
their counterpart in classical mechanics) as only a practical necessity and not as a manifestation of an 
inherent lack of complete determination in the properties of matter at the quantum level. 

Bohm's original papers do also provide a different path to find the trajectories by introducing the polar form of the 
wave function ip^fyt) = R{r,t)e^^^^'*^^^ into the (nonrelativistic) Schrodinger equation, Eq. (jsj. Let us emphasize 
that R{f, t) and S{r, t) are real functions. Then, after a quite simple manipulation, one obtains from the real part a 
quantum Hamilton-Jacobi equation; 

+ 2^ (V5(r, t)) ' + y (f, t) + Qir, t) ^ (5) 
where we have defined the quantum potential Q{r,t) as: 

which is the only difference with respect to the classical Hamilton-Jacobi equation if the velocity is defined as 
V = Such an alternative explanation provides an additional justification for the guidance equation, Eq. 

It allows a second-order interpretation of Bohmian trajectories in terms of acceleration, forces, and energies. 
In particular, in this second-order point of view, the new quantum potential Q{f,t) is responsible for the deviations 
of the Bohmian trajectories from the classical behavior that can be expected from the classical potential V{'r,t). As 
we will discuss later, the quantum potential of a system of N particles will be the responsible for the (nonclassical) 
nonlocal causality of Bohmian mechanics. 

Bohm completed the work of de Broglie in two fundamental aspects. First, as explained before, he demonstrated 
that Bohmian mechanics leads to exactly the same predictions as the ones obtained by orthodox quantum mechanics. 
Second, he provided a theory of measurement. He developed an explanation of the measurement problem without 
invoking the wave function collapse. The theory of Bohmian measurement will be discussed in Sec. |VI| Some authors 
argue that if we could change history, allowing Bohm to help de Broglie defend his pilot wave theory in the Solvay 
Conference, Bohmian mechanics would now be certainly taught at universities |43| . The enlightening work of Bohm, 
however, appeared 25 years too late, once the Copenhagen interpretation of quantum phenomena was already too 
well established. 

After the Solvay Conference, Bohr, Heisenberg, and their colleagues spread the new interpretation around the 
world and convinced the vast majority of the physics community that both quantum mechanics and the Copenhagen 
interpretation worked with extraordinary precision. A lot of young physicists were attracted to European institutes to 
study with the "fathers" of this new theory, and during the second quarter of the twentieth century, as good disciples, 
they spread the Copenhagen interpretation over the entire globe. In his 1976 Nobel lecture, Murray Gell-Mann 



13 See Sec. 



IV.3.4 



See problem |b] 



for the detailed formal demonstration. 
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referred to this question: "Niels Bohr brainwashed an entire generation of physicists into believing that the problem 
[of the interpretation of quantum mechanics] had been solved fifty years ago" [44] . 

Finally, we summarize the importance of Bohm's work with another magistral quote from Bell that appears in a 1982 
paper entitled "On the Impossible Pilot Wave" and collected in his famous book entitled Speakable and Unspeakable 
in Quantum Mechanics |14j : 

But in 1952 I saw the impossible done. It was in papers by David Bohm. Bohm showed explicitly 
how parameters could indeed be introduced, into nonrelativistic wave mechanics, with the help of which 
the indeterministic description could be transformed into a deterministic one. More importantly, in my 
opinion, the subjectivity of the "orthodox" version, the necessary reference to the "observer," could be 
eliminated. Moreover, the essential idea was on that had been advanced already by de Broglie in 1927, 
in his "pilot wave" picture. But why then had Born not told me of this "pilot wave"? If only to point 
out what was wrong with it? Why did von Neumann not consider it? More extraordinarily, why did 
people go on producing "impossibility" proofs after 1952, and as recently as 1978? When even Pauli, 
Rosenfeld, and Heisenberg, could produce no more devastating criticism of Bohm's version than to brand 
it as "metaphysical" and "ideological" ? Why is the pilot wave picture ignored in textbooks? Should it not 
be taught, not as the only way, but as an antidote to the prevailing complacency? To show that vagueness, 
subjectivity, and indeterminism, are not forced on us by experimental facts, but by deliberate theoretical 
choice? 



III. 7. John Bell and nonlocality 

Along this introduction, we have mentioned many different times Bell's positive opinion on Bohmian mechanics. 
Furthermore, his opinion about the Copenhagen interpretation was that the orthodox theory is "unprofessionally vague 
and ambiguous" |141 145H47I in so far as its fundamental dynamics is expressed in terms of "words which, however, 
legitime and necessary in application, have no place in a formulation with any pretension to physical precision" |47| . 

Bell spent most of his professional career at the European Organization for Nuclear Research (CERN), working 
almost exclusively on theoretical particle physics and on accelerator design, but found time to pursue a major avocation 
investigating the foundations of quantum theory. As seen in many of his quotes used in this introduction, his didactic 
ability to defend Bohmian mechanics against many unjustified attacks has been of extraordinary importance for 
maintaining the work of de Broglie and Bohm alive among the scientific community. Fortunately, Bell himself had his 
own reward from this unbreakable support of Bohmian mechanics. His outstanding work on locality and causality was 
directly inspired by his deep knowledge of Bohmian mechanics. Bell's theorem has been called "the most profound 
discovery of science" [15]. 

Bell's most relevant contribution to physics is probably the demonstration that quantum mechanics is nonlocal, 
contrarily to what Einstein expected. In 1964, inspired by the EPR paper [41j and Bohm's work on nonlocal hidden 
variables. Bell elaborated a theorem establishing clear mathematical inequalities, now known as Bell inequalities, for 
experimental results that would be fulfilled by local theories but would be violated by nonlocal ones [45 . Over the past 
30 years, a great number of Bell test experiments have been conducted. These experiments have confirmed that Bell's 
inequalities are violated (see, for example, Ref. (49j). Therefore, we have to conclude that quantum experimental 
results cannot be explained with local hidden variable theories. According to Bell, we must accept the real existence, 
in nature, of faster-than-light causation. 

The experimental violation of Bell's inequalities gave direct support, not only to the Copenhagen interpretation, 
but also to Bohm's formulation of quantum theory, since both are nonlocal theories. 

In the Bohmian case, as we have discussed, the quantum potential of A'^ entangled particles is defined in a SN 
configuration space so that an action on the first particle can have an instantaneous (i.e., faster-than-light) causal 
effect on the last particle. In Bohmian mechanics, one can understand that the quantum potential is responsible for 
the instantaneous nonlocal changes on the trajectories of quantum particles. Let Bell explain this point in his own 
words [H]: 

That the guiding wave, in the general case, propagates not in ordinary three-dimensional (3D) space but in 
a multidimensional configuration space is the origin of the notorious "nonlocality" of quantum mechanics. 
It is a merit of the de Broglie-Bohm version [of quantum mechanics] to bring this out so explicitly that it 
cannot be ignored. 

Unfortunately, as already happened with the von Neumann theorem jl3j, there is an historical misunderstanding 
about the consequences of Bell's theorems on Bohmian mechanics. Let us mention just one example that appeared 
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in a prestigious journal in 2000 [50] : 

In the mid-1960s John S. Bell showed that if hidden variables existed, experimentally observed probabili- 
ties would have to fall below certain limits, dubbed Bell's inequalities. Experiments were carried out by 
a number of groups, which found that the inequalities were violated. Their collective data came down 
decisively against the possibility of hidden variables. 

The author of this sentence has omitted the adjective "local" when he mentions hidden variables. Therefore, a 
confident reader, who has no time to read Bell's and Bohm's works, will understand that Bohmian mechanics is 
refuted by Bell's theorem. However, it is exactly the contrary. Bell's inequalities give direct support to Bohmian 
mechanics. Unfortunately, this misunderstanding appeared, and continues to appear, in many scientific articles, 
propagating into textbooks, websites, etc., provoking further comments and replies in the scientific literature. This 
discussion can give the impression that there still exists some controversy about the validity of Bohmian mechanics 
or that there is "something unclear" about it which is not the case. 

III. 8. Quantum hydrodynamics 

It is sometimes claimed that ideas similar to those developed by de Broglie were put forward by Madelung in 1926 
pS] , What Madelung proposed, however, was to regard an electron with mass m and wave function ilj{f,t) not as a 
particle with a determined trajectory but as a continuous fluid with mass density m |'!/'(r, |161 154] . In Madelung's 
hydrodynamic interpretation of Eq. ([s]), the fluid velocity coincides mathematically with de Broglie's guiding Eq. 
Q, but the physical interpretation is quite different. Let us return again to the double-slit experiment to understand 
Madelung's point of view. 

Orthodox quantum mechanics does not predict what happens to a single electron crossing a double slit, but it 
predicts what is the statistical probability of detecting electrons when we consider an infinite ensemble of such 
experiments. In contrast, the proposal from de Broglie and Bohm intends to predict what happens, in principle, to 
a single electron. At the end of the day, however, because it is not possible to determine the initial position of such 
a single electron with better uncertainty than that obtained from the initial wavepacket spatial dispersion, Bohmian 
mechanics also provides statistical results. In this regard, Madelung was not interested in dealing with a single-electron 
trajectory but with the ensemble. When we deal with a very large (infinite) number of particles (trajectories), the 
particle and current densities that appear in Eq. M are no longer discrete functions (sum of deltas) anymore but they 
have to be interpreted as continuous functions. Aerodynamics and hydrodynamics are examples of such continuous 
material systems. The concept of an individual trajectory becomes irrelevant in these disciplines but Eq. ([ij is still 
present. In summary, from a physical point of view, Madelung's proposal is completely orthodox. However, from 
a mathematical point of view, its formalism (and computational abilities) is very similar to de Broglie's proposal. 
Therefore, a reader who does not feel comfortable with the ideas of Bohm and de Broglie can assume Madelung's 
physical point of view and use many of the concepts explained in this book. 

III. 9. Is Bohmian mechanics a useful theory? 

Our historical introduction to Bohmian mechanics ends here. The reader has certainly noticed that we have some 
preference for the Bohmian interpretation over the Copenhagen, orthodox, explanation of quantum phenomena. How- 
ever, obviously, we have a profound and sincere respect for the Copenhagen interpretation of quantum phenomena 
and its extraordinary computing capabilities. Our opinion about what we have to expect from a physical theory that 
describes natural phenomena is quite pragmatic and close to the great fathers of the quantum theory. As Bohr said 
(see page 228 in [55]): 

In my opinion, there could be no other way to deem a logically consistent mathematical formalism than 
by demonstrating the departure of its consequence from experience or by proving that its predictions did 
not exhaust the possibilities of observations. 



See, for example, the experience of J. T. Gushing |51| or some recent works, "demonstrating" that Bohmian mechanics was wrong [52] . 
and the comment to the work |53| . 
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Curiously, similar arguments were used by Bohm to defend his causal theory against its many detractors (see page 
18 in 56 ): 



In conclusion, the author would like to state that we would admit only two valid reasons for discarding a 
theory that explains a wide range of phenomena. One is that the theory is not internally consistent, and 
the second is that it disagrees with experiments. 

In summary, from all the discussions done up to here, the only statement that wc will use in the rest of the book is that 
one cannot discard an interpretation of quantum phenomena in terms of Bohmian trajectories because it is internally 
consistent and it agrees with experiments. Once such a statement is accepted, other practical questions about Bohmian 
mechanics appear naturally, for example, is Bohmian mechanics a useful computational tool for predicting the results 
of quantum experiments? Although it provides identical predictions as the orthodox formalism, are there advantages 
in understanding/visualizing/explaining quantum experiments in terms of quantum trajectories? Or, even, is it a 
recommended formalism for thinking about limits and extensions of (nonrelativistic) quantum theory? This book, 
through several examples in the following chapters, will try to convince the reader that the answer to these three 
questions is affirmative. 



IV. BOHMIAN MECHANICS FOR A SINGLE PARTICLE 



After the previous historical introduction, we start with the formal presentation of Bohmian mechanics for a single 
particle. Usually, new scientific knowledge is built from small variations of old ideas. This explains why the initial 



development of quantum theory (see Sec. III. 3) was so traumatic. Quantum explanations did not evolve from a 



small variation of classical mechanics but from radically new ideas. The languages of both theories were completely 
different. Classical theory provides an explanation of a physical experiment on particles in terms of a well-defined 
trajectory, while the quantum orthodox theory needs a wave function. Here, we will see that both a wave description 
of an ensemble of classical particles and a description of quantum systems with trajectories are possible. Starting by 
setting a common language for both classical and quantum theories will certainly improve our understanding on the 
similarities and differences between themF^I 



IV. 1. Preliminary discussions 



For a reader without any previous knowledge on quantum mechanics, the following sections would be not a presenta- 
tion on Bohmian mechanics, but on quantum mechanics in general. Apart from their behavior during the measurement 
process discussed in section 1.4, the Copenhagen and the Bohmian wave functions have exactly the same time evolu- 
tion. In addition, the Bohmian wave function is complemented by a Bohmian trajectory which appears in a natural 
way when developing a trajectory-based explanation of quantum mechanics as we will discuss in section 1.2.3. 

The initial formulation of quantum mechanics with Bohmian trajectories developed by de Broglie and Bohm was 
perfomed for dynamical systems with associated velocities much slower than the speed of light. Although it is also 
possible to extend Bohmian mechanics to relativistic systems (see, for example, chapter 8 and Refs. [57Vl59j ). we will 
only deal in this chapter with nonrelativistic systems. 

In order to simplify as much as possible our mathematical notation, we will first study a single spinless particle 
living in a one-dimensional (ID) space. The spatial degree of freedom of the particle will be represented by x. The 
theoretical generalization of all the arguments mentioned in this section into a single particle in a 3D space is quite 
simple. However, the practical solution of a single-particle Schrodinger equation in a 3D space has considerable 
computational difficulties, as we will discuss in Sec. [Vjand appendix A.l. In Sec. V.4 we will consider the role of the 
spin. 



IV. 2. Creating a wave equation for classical mechanics 



In this section wc will derive a wave equation similar, although not identical, to the Schrodinger equation for an 
ensemble of classical trajectories. Certainly, this approach allows us to compare the quantum and classical theories 
by using a very similar language. 



It is very rare to find such descriptions in the standard literature but, in our opinion, it is very important to be able to compare classical 
and quantum mechanics on an equal footing. 
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IV.2.1. Newton's second law 



Our starting point to derive such a classical wave equation will be Newton's second law Let us consider a 

particle with mass m in a classical system that moves under the action of a potential V{x), where x is the position 
coordinate. Here, we assume the potential to be time independent so as to simplify the mathematical treatment. We 
define the particle trajectory as x\t] and its velocity as v\t] = dx[t]/dt = x[t]. Since we are considering a classical 
system, the trajectory of the particle will be obtained from its acceleration, a[t\ = <Px[t]/dt^ = from Newton's 
second law: 



dV{x) 
dx 



Since Eq. ([t]) is a second-order differential equation, we need to specify both the initial position a;[<o] 
initial velocity v[to\ = vq of the particle. 



(7) 

xo and the 



IV. 2. 2. Hamilton's principle 



Apart from Eq. (It]), there are other alternative ways to describe a classical system. For example, according to 
Hamilton's principle^GO] , the trajectory Xp [^P^ solution of Eq. (7^ between two different times, to and i/, provides a 
stationary value for the action function, S{x\t\;x(),to;xf,tf), where xq — '' " 
possible, we will omit the dependence of the action function on xq, Iq, Xf. 



x[tQ] and Xf — x[tf]. Hereafter, whenever 
and tf] 



6S{x[t] 
Sx[t] 







x[t]=Xp[t] 



where the action function is defined as: 



S{x[t]) 

The function L(x, v) is the Lagrangian function: 



L{x[t], x[t])dt 



L{x,v) = K{v) - V{x) 



(8) 



(9) 



(10) 



with K{v) = rnu^/2 the kinetic energy of the particle. The Lagrangian equation can also be defined as L{x,v) = 
pv — H{x, v), where H{x,p) = K{v) + V{x) is the Hamiltonian functiorj^ and p — mv is the (linear) momentum. 

Let us discuss the meaning of a stationary (or extremal) value of the integral in Eq. ([9|. We denote the physical 
trajectory solution of Eq. ([7| as Xp[t] and choose a shghtly different trajectory x^^lt] — Xp[t] + 5x[t] with the same 
initial and final conditions, that is, the same xq, io: xj, and tf. Therefore, we have (5a;[io] — and = 0. See 

Fig. [2] for a schematic representation of both trajectories. Equation ([8| means that the value of the action function 
'S'np — S{x[t\ + 6x[t];xo,ta;x f ^tf) will always be greateip^ than S^] that is, S'np > Sp. Thus, the trajectory that 
provides a stationary value for the action function is the physical trajectory, Xp[t]^ while a;np[t] ~ Xp[t\ + 5x[t] is a 
nonphysical solution incompatible with Eq. ([?]). 

The numerical evaluation of Eq. ([9]) to find Sp requires the prior knowledge of the trajectory Xp[t]. If we do not 
know the trajectory, we would have to evaluate <5'(a;[t]; xq, Iq; Xf,tf) for all possible trajectories (starting from xq and 
tg and ending ai Xf and tf) and choose the one that provides a stationary value. Certainly, it seems that Eq. Q has 
little practical utility. However, it provides an interesting starting point to develop our classical wave equation. 



IV. 2. 3. Lagrange's equation 

From Hamilton's principle it is possible to derive the Lagrange's equation |:60^ that gives us the differential equation 
that the physical trajectory satisfies. 



^'^ The suffix "p" means physical in order to distinguish from nonphysical trajectories, "np," but it will be omitted, when unnecessary. 
These definitions of Lagrangian and Hamiltonian functions are valid for the simple system described here. In any case, a different 
definition does not change the main results developed here. 

Strictly speaking, the physical trajectory can also correspond to a maximum of the action function. Such a maximum value is also an 
(stationary) extremal value of Eq. (Isl. 
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FIG. 2. Schematic representation of pliysical (solid lines) and nonphysical (dotted lines) trajectories in the {x, t) plane. The 
trajectory Xp[t] has the initial time to, the final time tf, the initial position a;p[to] ~ xq, and the final position a;p[t/] = Xf. The 
trajectory Xp-Sx[t] is a physical trajectory with identical initial and final conditions as Xp[t\, except for a different final position 
Xf + Sx. The trajectory a:;p_5t[f] is a physical trajectory with identical initial and final conditions as Xp[t] but reaching the final 
point at a larger time t + St. 



Let us define a trajectory a;np[i] = ^pM + Sx[t\ by adding a small arbitrary displacement Sx[t\ to the physical 
trajectory. In particular we fix (5a;[to] = and Sx[tf] = 0, so both trajectories have the same initial and final 



conditions. A Taylor expansion of the Lagrange function of Eq. (10) around the physical trajectory Xp[t\ reads 



L{Xnp[t],Xnp[t]) = L{Xp[t],Xp[t]) + 



dL 

dx 



6x[t] 



x=x-p[t\ 
±=ip[*] 



(11) 



The action function of Eq. (9| with the Lagrangian of Eq. (11) gives two contributions. The first term of the right 



hand side (r.h.s.) of Eq. ([TT) is Sp. After performing an integration by parts on the second part and using (5a;[io] = 
and = 0, we obtain: 



S{xnp[t\;xo,to;xf,tf)^Sp+ I 5x[t] y. 

Jt„ 



dL 

dx 



i=ip[t] 



d 

dt 



dL 

dx 



X = Xp[t] 

x=x„\t] 



dt 



(12) 



Hamilton's principle tells us that the action function must take a stationary value for the physical trajectory. This is 
equivalent to impose that small variations around the physical trajectory do not modify the value of the action, that 
is, for any (small) variation Sx[t], Snp = Sp. Thus: 



dL{x, x) 
dx 



X = X-p[t\ 

i=ip[t] 



d 

dt 



dL{x, x) 



dx 



X = Xp[t] 

x=x„\t] 



= 



(13) 



Equation ([13| is the Lagrange equation and gives, for a classical system, a differential equation that classical 
trajectories must fulfill. In order to enlighten the meaning of the Lagrange equation, Eq. (13), we can see that by 
substituting dL/dx = —dV{x)/dx and dL/dx = mv in Eq. (Ilo]), we recover: 



dV{x,t) 
dx 



X=Xj,[t] 
X=Xp[t] 



d 
dt 



dx [t] 

I ; 

dt 



x=Xplt] 
x=x-p[t] 



(14) 



which is the original Newton's second law of Eq. ([7|. In fact, what we have done is to check that Newton's second 
law is included in the fundamental Hamilton's principle. 
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IV. 2. 4- Equation for an (infinite) ensemble of trajectories 



The formalism based on the action function allows us to find not only a single physical trajectory but also the 
equation for an (infinite) ensemble of physical trajectories with slightly different initial or final conditions. For 
example, let us define: 



stl 



tf+St 



L{xp^st [t\,Xp^st [t])dt 



(15) 



to 



Xp_st\i\ being a trajectory with identical initial and final positions as Xp[t], but taking a longer time tf+St (see 
trajectory a;p_5t[t] in Fig. [2]). The new physical trajectory Xp_5t[t] can be written as Xp-st[t] = Xp[t] + Sx[t] for 
t G [to,tf]. For larger times, Xp-st[tf + 5t] = Xp[tf + St] + Sx[tf + St] with Xp[tf + St] = Xp[tf] + Xp[tf]5t, where we 
have done a first-order Taylor expansion of Xp[t] around tf for a small St. Since Xp^st[tf + St] = Xp[tf], we obtain: 



Sx[tf + St] = ~Xp[tf]St 
By following the same steps done after the Taylor expansion in Eq. (fTTl), we obtain 



p-St 



tf+5t 



Ldt 



to 



dL 
dx 



Sx[t] 



tf+St 



tf+St 



For the first term of the r.h.s. of Eq. (17), we obtain: 



tf+St 



Ldt 



to 



Sx[t] 



tf+St 



dL 
dx 



x—x[t] 
x—x[t] 



d 

di 



dL 
dx 



x—x[t] . 



dt 



Ldt- 



Ldt = Sp + LSt 



For the second term of the r.h.s. of Eq. (17): 

-I tf+St 



dL 

dx 



Sx[t] 



dL 
dx 



- x—x[tf+St] 



x—x[t f+St\ 



x{tf)St = -p[tf + St]x[tf]St 



(16) 



(17) 



(18) 



(19) 



where we use Eq. (16), fe[to] ~ and dL/dx = mv ^ p, p = mv being the particle momentum. 

Finally, since x[t] is also a physical trajectory that fulfills Lagrange equation, from to + St, we obtain for the 
third term: 



tf+St 



Sx[t] 



dt 
dx 



x—x[t] 



d_ 
dt 



dL 
dx 



x—x[t] 
x—x[t] 



dt = 



Putting together Eqs. (( 18 )) (( 19 )), and (( 20 )), we obtain: 



Sp-st = Sp + LSt — pxSt = Sp — HSt 
So, with dS/dtf = limst^o{Sp-st — Sp)/St we can conclude that: 

dS{xp[t] 



dtf 



= -H{xp[t],Xp[t]) 



(20) 



(21) 



(22) 



In summary, if Sp is the stationary value of the action function for Xp[t], then the variation of the new stationary 
value of another physical trajectory a;p_5([t], which has identical initial and final conditions but a slightly modified 
final time, is equal to the Hamiltonian (with a negative sign) evaluated at the final time of the trajectory Xp[t]. See 
problem [l] to discuss a particular example. 

Next, we will see that the evaluation of the value Sp^sx when we modify the final position Xf + Sx without modifying 
the initial and final times leads also to an interesting result (see Fig. [2]). Notice that the final position of the new 
physical trajectory Xp^sx[t\ = x[t] + Sx[t] means: 



Sx[to] — 0; — Sxf 



By following the same steps done after the Taylor expansion in Eq. (11 ), we obtain now: 



5, 



p—Sx 



Ldt 



to 



dL 
dx 



Sx[t] 



Sx[t] 



dL 

dx 



x—x[t] 
x—x[t] 



d 

di 



dL 



dx 



dt 



=x[t]^ 



(23) 



(24) 



23 



The first term of the r.h.s. of Eq. (24) is again Sp. The second term is 



- */ 


'dL' 


- to 


dx 



x—x[t j-] 



x=x[tf 



Sxf ^p[tf]Sxf 



(25) 



where we have used Eq. (23). The third term is zero. We can conclude that dS/dx — limsx^o{Sp-Sx — Sp)/Sx is 
equal to: 



dS{xp[t];xa,ta;x,t) 
dx 



mxp[t] = pp[t\ 



(26) 



In summary, the variation of the stationary value of the action function when we slightly modify t he fi nal position of a 
physical trajectory is equal to the momentum of the trajectory at the final timcj^ To check Eq. (26) for a particular 
case, see problem [2j 



IV. 2. 5. Classical Hamilton- Jacobi equation 



Now, using Eq. (22) and substituting the velocity (momentum) into the Hamiltonian by Eq. (26), we obtain the 



well-known Hamilton-Jacobi equation: 



ox 



dS{x,t) 
dt 



(27) 



We want to use the Hamilton-Jacobi equation, Eq. (27), to find all the physical trajectories. For this reason, we 



have eliminated the dependence of S on the trajectory x[t]. In addition, we have also eliminated to and Xq because 
we consider that x and t are not parameters but variables. Notice that Eq. (27) is valid for physical trajectories so 



that once we know S{x,t), we are able to compute directly all physical trajectories for all possible initial conditions. 
Therefore, Eq. (27) defines an (infinite) ensemble of trajectories rather than just a single trajectory. However, in most 
practical cases, the direct solution of the Hamilton-Jacobi equation is much more difficult than using the Newton or 
Lagrange formulation of classical mechanics. It has limited practical interest. However, it provides a direct theoretical 
connection with a classical wave equation similar to the Schrodinger equation. See problem [3] for a particular solution 
ofEq. Elh. 



IV. 2. 6. Local continuity equation for an (infinite) ensemble of classical particles 



Equation (27) can be interesting when dealing with an ensemble of trajectories. For example, it can be used when 
we have a classical (single-particle) experiment with some practical difficulty in specifying the initial position and 
velocity of the particle, such that different experimental realizations can have slightly different initial conditions. 



The ensemble of trajectories solution of Eq. (27) will provide a statistical (probabilistic) description of the classical 
(single-particle) experiment that accounts for the variability in the initial conditions. 

For such an experiment, we can reasonably assume that the variability of the initial position will be limited to a 
particular spatial region. At the initial time — 0, we can define some distribution of the initial position of the 
particles i?^(a;,0) > 0[^ Such an ensemble of trajectories will evolve in time according to Eq. (27) so that we will 
obtain a function R^{x,t) that describes the particular distribution of particles at any time. The function R^{x,t) 
is constructed by "counting" the number of trajectories present inside the interval (a;, x + dx) at time t. From Eq. 
(26), we define the particle velocity as i[t] = {l/m)dS{x,t)/dx. Finally, we know that all these classical particles will 



move, in a continuous way, from one unit of volume to another. Therefore, as we have extensively discussed in Sec. 
|III.1[ we can ensure that the ensemble of trajectories accomplishes the following local conservation law: 



dR^{x,t) 
dt 



d_ 
dx 



1 dS{x,t) 
m dx 



R^{x,t) = 



(28) 



In summary, Eq. (28) just certifies that if a classical particle goes from one point to another, it has to go through all 



the trajectory positions between these two places. 



We have recovered the subindex "p" to emphasize that Eq. | |26| l is only vaUd for physical trajectories. 

Such a distribution is positive (or zero), but it is not necessary to impose that B?(x,Ci) is normalized to unity. 
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IV. 2. 7. Classical wave equation 



Now, we have all the ingredients to develop a wave equation for classical mechanics. In previous paragraphs we 
have been dealing with two (real) functions, S{x,t) and R{x,t). The first, S{x,t), is the action function and appears 
in the Hamilton-Jacobi equation, Eq. ( p7| . In particular, its spatial derivative determines the velocity of the particles 
(see Eq. (26)). The second, R{x,t), tells us how an ensemble of trajectories is distributed at each time t 



It evolves 

In this sense, we can assume that they are some kind of "field" or "wave" 



according to the conservation law, Eq. 
that guide classical particles. 

We can construct the following classicaj^ (complex) wave function ipci{x,t) — R{x,t) exp{iS{x,t) /h). We divide 
the action function inside the exponential by h (the reduced Planck's constant) in order to provide a dimensionless 
argument. Then, it can be shown that the two previous (real) equations, Eqs. ((27)) and (([28|), for S{x,t) and 
R{x,t) are equivalent to the following (complex) classical wave equation for (x,t): 



ih 



dtp ci{x,t) 
dt 



h^d^^Mx^ 

2m dx^ 



+ V{x,t)ijci{x,t) 



d^\Mx,t)\/dx^ 



2m 



1pcl{x,t) 



(29) 



Hcl{x,t)\ 

The demonstration of this expression is left as an exercise to the reader (see problem |4|. Additionally, see problem [s] 
for a solution of the classical wave equation, Eq. (29), for a free particle. 



In conclusion, an ensemble of classical trajectories can be described with a wave function solution of a wave equation. 
Thus, a common language for classical and quantum mechanics has been obtained by using an (infinite) ensemble of 
classical trajectories with different initial positions and velocities instead of just one single classical trajectory. 



IV. 3. Trajectories for quantum systems 



Before comparing the classical and quantum wave equations, let us first discuss in this section whether trajectories 
can be also used to describe quantum systems. We will introduce such trajectories in two different ways: first, as a 
direct consequence of the local conservation of particles extracted from the Schrodinger equation and, second, directly 
following the work presented by Bohm in his original paper [8] . 



IV. 3.1. Schrodinger equation 

As discussed in Sec. |III.3[ the single-particle Schrodinger equation in a ID quantum system subjected to a scalar 
time-dependent potential, V{x,t), is: 

.^dih(x,t) h? d'^tjAxA) , , 

rn^^--^^^Vix^mx^t) (30) 



It is important to emphasize that in the orthodox interpretation of ^p{x,t), Eq. (30) does not describe a single 
experiment but an ensemble of identical (single-particle) experiments. The orthodox meaning of the square modulus 
of the wave function |'!/;(a;, t)p is the probability density of finding a particle at position x at time t when a measurement 
is performed. Such probabilities assume an infinite number of identical experiments. 



IV. 3. 2. Local conservation law for an (infinite) ensemble of quantum trajectories 



It is known that there is a local continuity equation associated with Eq. (30). Let us first mathematically derive it 
to later discuss its physical implications. 

In order to find a local continuity equation, let us work with ipix, t) and its complex conjugate 4'*{x, t). In particular, 
we can rewrite Eq. (30) as: 



Tp* (x, t)ih 



dip{x, t) 
dt 



d'^tbix t) 

-r{x,t)— " > +r{x,t)vix,mx,t) 

2m ox'^ 



(31) 



For simplicity, we assume that R{x,t) and S{x,t) are single valued. This condition is equivalent to assigning a unique velocity to each 
position of space. If this is not the case, then we will need different wave functions, that is, a mixed state, to describe the different 
velocity fields of the classical system. 
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and the complex conjugate of Eq. (31 ) as: 



ip{x, t)ih 



dt 



2m dx'^ 



+ ^P{x,t)V{x,t)^P*{x,t) 



From Eqs. (31) and (32), we obtain: 

dmx,t)\^ 

dt 



n d 



2mdxV^'' dx dx 



(32) 



(33) 



We can easily identify Eq. (33) as the local conservation of particles discussed in Eq. ([I| where p{x,t) — \^p{x,t)\'^ 
and the current density, J{x, t), is defined as: 



J{x,t) 



2m 



dx 



dx 



(34) 



Unlike other wave equations, the Schrodinger equation is compatible with a local conservation of particles due 
to the fact that V{x,t) is a real function. We have noticed above that we can interpret p{x,t) = |V'(a:,t)P as a 
spatial distribution of an ensemble of trajectories. Such an ensemble of particles can be obtained by repeating a 
single-particle experiment and measuring, at each time t, the number of particles at each x. The presence of such 
local conservation of particles is very relevant for us because it justifies our aim to look for an ensemble of continuous 
trajectories describing p{x,t) = \ip{x,t)\'^. 



IV. 3. 3. Velocity of Bohrman particles 



In Bohmian mechanics, the particle velocity is defined as: 

v{x, t) 



J(x,t) 



(35) 



where J{x,t) is given by Eq. (34). Taking into account that |V'(a:^,i)P is the distribution of the ensemble of particles, 
it is easy to show that this velocity definition is compatible with the local continuity equation, Eq. ( 33 ) , and that an 



ensemble of well-defined trajectories whose initialpositions are all selected according to the distribution |V'(x,fo)P 
will reproduce |V'(a^,t)P at all times (see problem]?). 

Notice that one could add a divergence- free term to this velocity, and it would still fulfill Eq. ( 33 ) . Even though 
this would change the path followed by each individual trajectory, the trajectory distribution and, therefore, the 
measurement outcomes would remain unchanged. However, if one develops the nonrelativistic limit of a relativistic 
treatment of Bohmian mechanics, one finds that the divergence-free term does not appear for spin-O-particles and 
that the velocity is unambiguously defined by Eq. (35). For other values of spin, the situation is a bit more complex 
[6T| [62] and far from the scope of this book. 



IV. 3. 4- Quantum Hamilton-Jacobi equation 

Following the path described by Bohm in his original paper [2] , we will now start from the quantum wave equation, 
that is, the Schrodinger equation, to arrive at a "quantum" Hamilton-Jacobi equation that describes the motion of 
quantum trajectories. 

The first step is to write the quantum (complex) wave function, ip{x, t) = ipr{x, t) -\- ii^i{x, i), in polar form: 

i?2 (x, t) = i^l {x,t)+ V? {x, t) (36) 

S{x, t) = h arctan ( "^'^"^'^M (37) 

In principle, S{x,t), the so-called quantum action, is not well defined when tpri^x^t) — ipi{x,t) ~ 0, that is at those 
points where R{x,t) = 0, meaning that no particles will reach themP^ 



We assume that the wave function is single valued so that R{x,t) is also single valued. However, the definition of S{x,t) has some 
practical difficulties. In principle, S{x, t) is a multivalued function because the function arctan(a::) itself is a multivalued function. If we 
want to use Eqs. (( |36| ) and ((jSTj) to reconstruct the wave function, then the multivalued problem can be eliminated by imposing an 
additional restriction on the definition of S |20l 1241 163| . 
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The quantum Hamilton- Jacobi equation can be found by introducing tfj{x,t) = R{x,t) exp{iS{x,t) /h) into Eq. 



(30). On the one hand, the imaginary part of the resulting equation gives the local conservation law identical to the 
one shown in Eq. (28), which we rewrite here for convenience: 

ai?2(x,t) 



ut ox \m ox 

On the other hand, the real part gives a quantum Hamilton-Jacobi equation: 



dSix,t) 
dt 



1 
2m 



dS[x,t) 
dx 



+ V{x,t) + Q{x,t) = 



(38) 



(39) 



Since the last term in the r.h.s. of the classical wave equation, Eq. (29), is not present in the Schrodinger equation, 
Eq. (30), an additional term appears in the quantum Hamilton-Jacobi equation, the so-called quantum potential, 
defined as: 



Q{x,t) 



d^R{x,t)/dx^ 



2m 



R{x,t) 



(40) 



In conclusion, identical to the classical description of a system obtained from Eq. ( 39 ) , we obtain an interpretation 
of the wave function solution of the Schrodinger equation as an ensemble of quantum trajectories with different initial 
positions and velocities. The velocity of each trajectory x[t] is defined as: 



v[t] 



1 dS{x,t) 



dx 



(41) 



x—x[t] 



Interestingly, it can be easily shown that this new expression for the quantum velocity is identical to that mentioned 
in Eq. (35). See problem |6] to show that both definitions of the velocity of Bohmian trajectories are identical: 



v{x, t) 



1 dS(x,t) J{x,t) 



dx 



(42) 



where J{x,t) is defined by Eq. (34). 



IV. 3. 5. A quantum Newton-like equation 



It is interesting to compute the time derivative of the Bohmian velocity defined in Eq. (41) to find a quantum 
Newton-like equation: 



d d 

m—-v(x\t\,t) = m—- 
dt ^ ^ ' ' dt 



' 1 OS' 






m dx 


x—x[t] 


dx"^ 



m 



x—x[t\ 



d dS 
dx dt 



x—x[t\ 



(43) 



We can rewrite Eq. (43) as: 



m-^^v{x[i\,t) 





[1- 




ds\ 


dx 


\2m 


\dx J 


" ~dt j 



(44) 



Finally, using Eq. (39) we find: 



■m—v[x[t],t) 



-^[V{x,t)+Q{x,t)) 



=x[t\ 



(45) 



We conclude here the second route for finding a common language for classical and quantum theories. The quantum 
(complex) single-particle wave function can be interpreted as an ensemble of trajectories that are all solutions of the 
same single-particle experiment but with different initial conditions. The quantum trajectories are not solutions of 
the classical Newton second law with a classical potential but solutions of the quantum Newton second law, Eq. ( 45 ) , 
where a quantum potential (that accounts for all nonclassical effects) is added to the classical potential. 
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IV. 4. Similarities and differences between classical and quantum mechanics 



In previous sections, we have provided a common language for classical and quantum theories, in terms of either 
wave functions or trajectories, to fairly compare differences and similarities between both theories. Here we emphasize 
that we are comparing either classical and quantum wave functions or classical and quantum ensembles of trajectories 
(not a single classical trajectory with a quantum wave function). The main difference between the mathematical 
description of the two theories is that the term Q{x, t) that appears in the quantum Hamilton-Jacobi equation, Eq. 



(39), is exactly the same term that appears in the classical wave equation, Eq. (29), but with a change of sign. The 



term Q(x,t) explains the exotic properties of quantum systems that are missing in their classical counterparts: 
1. Differences 

(a) Quantum superposition 



One can apply the superposition principle to find solutions of the Schrodinger equation, Eq. ( 30 ) , but it 



cannot be applied to the classical wave equation, Eq. ( 29 ) , since its last term is nonlinear. In fact, if -01 (a;, t) 
and ip2{x,t) are solutions of the Schrodinger equation, then aipi{x,t) + bip2{x,t) with arbitrary complex a 
and b is also a solution. Quantum mechanics will keep the (complex) amplitudes a and b constant at any 
time because of its linearity, which is not the case for classical mechanics. 

This is a fundamental difference between classical and quantum systems, which has deep consequences. 
Quantum mechanics can be developed in a linear vector space, while classical mechanics cannot. 

(b) Quantum wholeness 



The presence of Q{x, t) in the quantum Hamilton- Jacobi equation, Eq. (39 ), implies that Bohmian trajec- 
tories depend not only on the classical potential V{x, t) but also on the quantum potential Q{x, t), which is 
a function of the quantum ensemble of trajectories distribution, R(x, t). In fact, it is the shape and not the 
absolute value of R{x, t) that acts on each individual quantum trajectory. On the contrary, each classical 
trajectory can be computed independently of the shape of the ensemble. 

This surprising result for quantum mechanics can be illustrated with the double-slit experiment |64| . We 
assume that the initial wave functions for the classical and quantum ensembles are identical at to = 0- Their 
difference appears in the time evolution of the trajectories. For the quantum trajectories of the ensemble, 
the shape of the ensemble (whether particles are stopped or not by the double-slit screen) determines the 
shape of R{x,t), which will affect the dynamics of all trajectories, even those trajectories that are far 



from the slit (see Fig. Km). On the contrary, the classical Hamilton-Jacobi equation, Eq. (27), is totally 




FIG. 3. A particle is sent toward a thin plate with two slits cut in it. This single-particle experiment is repeated many times 
with some uncertainty on the initial conditions of the particle positions. The distribution of detected particles on the screen 
as a function of the position is different if we deal with a quantum (a) or a classical (b) system. In both cases, we use an 
identical description of the initial ensemble of particles, that is, identical R{x,to) and S{x,to). However, the classical evolution 
of the trajectories Xci[t] is d escribed by Eq. (27), which is independent of R{x,t), while quantum (Bohmian) trajectories XB[t] 
are determined by Eq. (39 1, which depends explicitly on the other particles of the ensemble through Q{x,t) that depends on 
R{x,t). These differences can also be understood by observing that the evolution of a classical wave function is described by a 
nonlinear equation, Eq. (1291), while a quantum wave function evolves with a linear one, Eq. (1301). 
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independent of R{x,t), so a single trajectory is completely independent from the rest of particles (see Fig. 
St 



2. Similarities 

(a) Uncertainty 

As far as we deal with an ensemble of (classical or quantum) trajectories, there is an uncertainty in the 
exact value of any magnitude that can be measured from the ensemble. For example, one can compute the 
mean value and the standard deviation of the position of the classical ensemble. The classical and quantum 
uncertainties can have different origins, but both ensembles have uncertainties. The classical uncertainty 
is due to the technical difficulties in exactly repeating the initial conditions of a particular experiments; 
the quantum uncertainty has a more intrinsic origin |65l I66j . 

Even assuming their different physical origin, what we want to emphasize here is that it is not licit to 
compare the uncertainty of a quantum wave function with the uncertainty of a single classical trajectory, 
saying that the classical trajectory has no uncertainty and the quantum wave function has. As we have 
discussed, we have to compare the uncertainty of classical and quantum wave functions (or classical and 
quantum ensembles of trajectories). 

(b) Initial conditions for the (classical or quantum) Hamilton Jacobi equation 

The classical Newton equation is a second-order differential equation, where both the initial position and 
the velocity have to be fixed. In the classical or quantum Hamilton-Jacobi equation, it seems that only an 
initial position is needed because the initial velocity is directly determined by the spatial derivative of the 
action. However, one can argue that in the Hamilton-Jacobi equation, two initial conditions are fixed, the 
initial position and the initial wave function (or initial action) that fixes the initial velocity. In other words, 
for one particular Hamiltonian, even if we fix one particular initial position, then it is possible to obtain 
different initial velocities if we select different initial wave functions. In conclusion, the discussion that 
classical trajectories are solved from second-order equations, while Bohmian trajectories from first-order 
equations, is somehow artificial if one recalls the quantum equilibrium hypothesis 

(c) Single-valued wave function and multivalued action function 

In principle, the action function solution of the (classical or quantum) Hamilton-Jacobi equation can be 
multivalued yielding different velocities in a particular x and t point. However, the classical wave or the 
Schrodinger equations deal with single-valued wave functions. If we want to model (classical or quantum) 
scenarios with multivalued velocities, then we have to consider several wave functions, one wave function 
for each possible velocity at x and trj On the contrary, if we only work with a single-valued (classical or 
quantum) wave function, then, the velocity itself is single-valued everywhere. Therefore, two (classical or 
quantum) trajectories that coincide in one configuration point will have identical velocities. This means 
that they will follow identical trajectories for any future time. This has the important consequence that all 
Bohmian trajectories (or classical trajectories) associated with a single-valued wave function cannot cross 
in the configuration space. 

In conclusion, the differences between quantum and classical ensembles of trajectories is not a difference between 
waves and particles, because both waves and particles can be used to study classical or quantum systems. On the 
contrary, the difference resides between a linear wave equation (for quantum mechanics), Eq. (|30[), and a nonlinear 



wave equation (for classical mechanics), Eq. (29). One of the most important consequences of such difference is that 
quantum (Bohmian) trajectories depend on the shape of the ensemble (i.e., quantum wholeness), as seen in Fig. [3^, 
while classical trajectories are independent of the shape of the ensemble, as seen in Fig. [Sja. This difference between 
classical and quantum systems has important consequences at a computational level. One can compute a unique 
classical trajectory. However, in principle, one has to compute the whole ensemble of quantum trajectories even when 
a unique quantum trajectory is required. 



IV. 5. Feynman paths 

In the introduction of this chapter we mentioned that quantum mechanics can be described either in the matrix 
formulation proposed by Heisenberg and coworkers or in the wave equation formalism developed by Schrodinger. There 



The quantum equilibrium hypothesis assumes that the initial positions and velocities of Bohmian trajectories are defined distributed 
according to the initial wave function. However, there are authors who do not want to take this assumption and then select initial 
positions and velocities not related with the initial wave function 17 65-67 . This topic will not be further developed in this chapter. 
In the quantum language, this means working with a density matrix, that is, with mixed states rather than with pure states. 
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FIG. 4. Schematic representation of physical (sohd Unes) and nonphysical (dashed hnes) trajectories in the configuration space 
{x,t) starting at {xo,to) and finishing at {xf,tf). In general, Feynman paths XFy[t] do not provide stationary values of the 
action function, while the quantum (Bohmian) trajectory XB[t] provides a stationary value of the quantum action function, Sb 
in Eq. ( |49[ ), and the classical trajectory Xci[t] a stationary value of the classical action function, Sd in Eq. ( |48| |. 



is still another alternative representation of quantum phenomena due to Feynman (see, for instance, Ref. [68] )■ The 
time evolution of a wave function can be written using the Green function (or propagator or transition amplitude) as: 



'ip{x,t) 



G{xo, to; X, t)tp{xo, to)dxo 



(46) 



Feynman provided an original technique for computing G{xo,to; x,t) from the classical Lagrangian of Eq. (10 1. One 
considers all (physical and nonphysical) paths XFy[t\ that may link the two points {x, t) and (xq, Iq). See dashed curves 
in Fig. B] Feynman associated to each path (physical or nonphysical) a complex amplitude exp(iS'(a;[i]; a;o, io! tj/fi), 
where Slx[t];xQ,tQ;x,t) is defined by Eq. ([9]). We emphasize that we are using not only the trajectory that provides 
an stationary value of S{x[t\; XQ,tQ; x,t) but all trajectories. Then, Feynman defines the Green function as: 



G{xo,to;x,t) = C 



E < 

all paths 



,iS {x[t];XQ ,to;x .t) / h 



(47) 



where C is a normalization constant. It can be demonstrated that the wave function constructed from Eq. (46) 
with Eq. (47) reproduces the Schrodinger equation [5H]. Certainly, Feynamn paths and Bohmian trajectories are 
completely different. In particular, there is only one Bohmian trajectory XB[t] solution of Eq. ( |39[ ) that goes from 
point (xo,to) to point {xf,tf), as seen in Fig. |4j while there are infinite Feynman paths a;Fy[t] that connect both 
points. In particular, one of the Feynman paths is certainly xslt] and another is the classical trajectory Xci[i]. In 
general, XB[t\ is different from a:c;[i]. The later is the trajectory that minimizes the classical action function: 



S{xci [t];xa,tQ;xf,tf) 



to 



2m 



{x[t]r~v{x[t],t) 



dt 



while the former is the trajectory that minimizes the quantum action function: 

S{xB[t];xo,to;xf,tf) 



1 



2m 



(i[ 



V{x[t],t)~Q{x[t],t) 



dt 



(48) 



(49) 



The derivation of the previous expression can be straightforwardly obtained by repeating the steps done from the 
classical action function to the classical Hamilton- Jacobi equation but using now the quantum Hamilton-Jacobi 
equation. We conclude that Bohmian trajectories are quite different from Feynman paths. 

On the other hand, one could be interested in discussing whether the Feynman approach could be a more useful 
computational tool than Bohmian mechanics for solving quantum systems. In principle, it seems that Feynman paths 
would be less practical because then we would have to compute an infinite number of trajectories from one particular 
initial position and then repeat the procedure for all initial positions. On the contrary, in Bohmian mechanics we 
only need to compute one quantum trajectory for each initial position. 
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IV. 6. Basic postulates for a single-particle 



After introducing the reader to the language of trajectories for quantum mechanics (and wave functions for classical 
mechanics), we can state the basic postulates of the Bohmian theory for a single particle. They synthesize in a few 
sentences all the discussions done up to here. In general, the postulates of any physical theory can be presented in 
different compatible ways. For example, classical mechanics can be postulated from Newton's laws or Hamilton's 
principle. We follow here the standard presentation of Bohmian mechanics that involves a small number of ingredients 

[m [501 [Ml Eg. 

FIRST POSTULATE: The dynamics of a single particle in a quantum system is defined by a trajectory x[t] that 
moves continuously under the guidance of a wave function ^{x,t). The wave function ip{x,t) is a solution of the 
Schrddinger equation: 

'^^r- = -2m^x^ + 
The trajectory x[t] is obtained by time-integrating the particle velocity v[t] defined as: 

, , J(x,t) 
\ip{x,t)\^ 



where J{x,t) is the (ensemble) current density given by: 



2m \ ' 9x ^ dx 

The initial position x[t{)\ and velocity v[to] must be specified to completely determine the trajectory x[t]. 

SECOND POSTULATE (quantum equilibrium hypothesis): The initial position and velocity of a particular trajec- 
tory cannot be known with certainty. When the experiment is repeated many times {j — 1, ■ ■ ■ , M), the initial positions 
{x^[to]} of an ensemble of trajectories {x^[t]}, associated to the same il){x,t), have to be generated so that the number 
of trajectories of the ensemble between x and x -\- dx at the initial time to is proportional to R'^{x,to) — \tp{x,to)\'^ . 
The initial velocity of each trajectory is determined by v^[tQ] — J(x^[to],to)/\ip{x^[to],to)\'^. 

The condition on the initial position can be written mathematically as: 

1 

R^{x,to)^ lim y^S{x-x^[to]) (50) 

J = l 

where j — 1, . . . , M is the number of different trajectories of the ensemble. This second postulate and the continuity 
equation, Eq. (33), ensure that an ensemble of Bohmian trajectories will reproduce R{x,t)'^ = \ip(x,t)\'^ at any time: 



M 

{x, t) = Jhn^ ]^ E - W ) (51) 

See problem [7] This is precisely the reason why it is claimed that Bohmian mechanics exactly reproduces the position 
measurement of orthodox quantum mechanics. 

These postulates represent a minimalist explanation of the Bohmian interpretation of quantum mechanics without 



mentioning either the Hamilton- Jacobi equation, Eq. (39 1 or the quantum potential, Eq. \AQ\. Certainly, we can 



formulate Bohmian mechanics without the use of the Hamilton- Jacobi formalism; however, in the authors' opinion, the 
quantum Hamilton-Jacobi and the quantum potential allows us to improve our understanding of Bohmian mechanics 
and provides clear arguments for discussing the similarities and differences between classical and quantum theories. It 
is argued by some authors that the second postulate, the quantum equilibrium hypothesis, is not a necessary postulate 
of Bohmian theory |17[ I21j . It is stated that an arbitrary initial distribution of positions and velocities can also arrive 
at quantum equilibrium quite rapidly. This topic however is far from the scope of this introductory chapter, and it 
will not be considered in this book. 

Finally, let us emphasize that no postulate about measurement is presented. In Bohmian mechanics no postulate 
about measuring is needed. For a more detailed explanation, please refer to Refs. [171 [2Ql [63l [MI and Sec. [VI 



In next section, we will explain how these postulates presented here for single-particle systems are generalized to 
many-particle systems. 
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V. BOHMIAN MECHANICS FOR MANY-PARTICLE SYSTEMS 
V.l. Preliminary discussions: The many body problem 

Up to now we have only studied a simple system composed of just one particle. However, a single-particle quantum 
system is some kind of crude idealization of a natural system formed by very large number of particles. In this section 
we will study nonrelativistic many-particle quantum systems with Bohmian trajectories. First of all, let us clarify the 
meaning of the term "many-particle" in Bohmian mechanics, as it can easily be misinterpreted. As we have explained 
in the previous section, we need an infinite ensemble of trajectories, that is, with AI — ^ oo, to describe the statistics 
of a single-particle quantum system. These different trajectories should not be confused with different particles, since 
they all refer to different dynamical realizations of the same single-particle quantum system. 

Let us consider now a quantum system with N degrees of freedom, that is, an N body quantum system. We will 
use a particular variable Xk for each degree of freedom k = 1, . . . , iV. The wave function will thus be a function of all 
xi, . . . , xn variables. Now, a many-particle Bohmian trajectory will involve N interacting particles xi[t\, . . . , XN[t\- 
The relevant point that allows us to use the adjective "many-particle" is that the N particles interact between them, 
that is, the interacting potential depends on all possible particle positions, V{xi, . . . ,xn) in a non-trivial way. Along 
this section, in order to simplify our notation, we will use x[t\ = {xi[t\, . . . ,X]^[t\) and or x = (xi, . . . ,xpf). 

Now, an ensemble of Bohmian trajectories will be composed of M — )• oo trajectories x^[t], with j — 1,...,M. 
The superscript j refers to the statistical index (counting trajectories), while the subscript k refers to each of the N 
interacting particles, that is, x^t] = {x{[t], . . . , x-'j^[t]). 

The first step to obtain the many-particle quantum trajectories is solving the following many-particle Schrodinger 
equation: 

^dMx,t) /A d"^ A , 

The solution ip{x,t) of this equation is the so-called many-particle wave function, that is defined in a iV-dimensional 
space (plus time) . The problem of N particles in a ID space is formally equivalent to the one of a single particle in an 
A'^-dimensional space. Here, we will use the term "many-particle" in a wide sense to include the single-particle N > 2 
dimensional caser3 



Equation (52) is analitically unsolvable most of the times, and its numerical integration is out of today's present 
computer capabilities, even for systems with N 5, since we need to compute the wave function Tp{x,t) in the A'^- 
dimensional configuration space. Let us roughly estimate the hard disk space that we would require to store ip{x,t). 
Considering, for instance, a system with = 10 particles confined in a ID region of 10 nm, which we discretize with 
a spatial step of Ax = 0.1 nm, we have a grid of 100 points for each dimension Xk- Then, the total number of points 
in the configuration space for the 10 particles is 100^" = 10^°. Using 4 bytes (32 bits) to store the complex value of 
the wave function at each grid point, the information contained in a 10-particle wave function would require more 
than 3 x 10* Terabytes (TB) (and more than 3 x 10^* TB for 20 particles). This practical limitation is the reason 
why our knowledge of many-particle quantum systems is quite poor. In 1929, Dirac wrote the following |70| : 

The general theory of quantum mechanics is now almost complete. The underlying physical laws necessary 
for the mathematical theory of a large part of physics and the entire chemistry are thus completely known, 
and the difficulty is only that the exact application of these laws leads to equations much too complicated 
to be soluble. 

Three decades later. Born rephrased the issue [7T| : 

It would indeed be remarkable if Nature fortified herself against further advances in knowledge behind the 
analytical difficulties of the many-body problem. 

The adjective "many-particle" can certainly be used with classical mechanics for N particles interacting through 
V{x,t). Now, it would be interesting to address the previous problem of storing the state of a 10-particle system 



The variable TV can be defined as the number of particles in a ID space, or it can be related to the number of particles in a 3D space. 
In simple words, ip{xi, X2, X3,t) can be interpreted as three particles in a ID space or just one particle in a 3D space. From a physical 
point of view, one particle in a 3D space is a "single-particle" system. However, from the computational point of view, it is equivalent 
to a three-particle system in ID. 
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but now from the classical mechanics perspective. Solving the problem means finding the many-particle trajectory 
x'^^lt]. As discussed at the end of subsection 1.2.4, a classical many-particle trajectory can be computed alone, without 
knowing the rest of the trajectories of the ensemble (with different initial positions). At a computational level, this 
feature of classical mechanics implies a dramatic simplification. Basically, at each time step of the simulation, if the 
analytical expression for the potential is known, we just need to save the position and velocity of each particle. This 
means 20 real values for the 10-particle systems mentioned above. The classical trajectory of each particle, can 
be obtained by solving a coupled system of Newton's second laws: 



dV{x, t) 



dxk 



(53) 



3=Sii [t] 



Fortunately, to compute a single classical trajectory x^'ft], we only need to evaluate the derivative of the many-particle 
potential V{x, t) along the trajectory. In the quantum case, even finding a single trajectory means solving for the 
entire wave function ^[^{x, t) that includes the information of the entire ensemble of possible trajectories. To be fair, let 
us mention that the whole ensemble of classical or quantum trajectories will require similar computational resources. 

In the scientific literature, there are many attempts to provide reasonable approximations to the many-body quan- 
tum problem. Density functional theory (72 [73] and the Hartree-Fock approximation [73H7S] are some of the most 
popular techniques among the scientific community for dealing with the many-body problem. In Sec 1.3.6 we will 
discuss whether Bohmian trajectories can help in solving the many-body problem. In what follows, we present the 
basic theory of the many-particle Bohmian trajectories. 



V.2. Many-particle quantum trajectories 



The efforts done in the single-particle section to compare classical and quantum mechanics are valid here. In par- 
ticular, we could write a classical (many-particle) wave equation or develop (many-particle) trajectories for quantum 
systems. In this section, we will only explicitly develop the latter. We start by considering (nonrelativistic) spinless 
particles. In Sec. V.4 we will introduce spin. 



V.2.1. Many-particle continuity equation 



Following a similar development as we did for the single-particle case, see Sec. |IV.3 



local continuity equation associated to the many-particle Schrodinger equation, Eq. (52) 



we can derive the foUowins 



dt ^ dxk 

fc=i 



5]— J,(f,t) = (54) 



where we have defined 



^'=(^'')=^2;^(,^(^'')^^-^ ^"^'^^^^j ^^'^ 

as the fc-th component of the current density, see problem [T4j 

Since the many-particle Schrodinger equation is also compatible with a local conservation of particles, we can 
interpret ["(/"(a^, ^)P the spatial distribution of an ensemble of trajectories. The Bohmian velocity of the k-th 
trajectory is: 

v,ix,t) = ^^0^ (56) 

In fact, the strategy followed here to develop Bohmian mechanics can be extended to any quantum theory where a 
continuity equations holds: first, look for a continuity equation for the probability density and, then, define a velocity 
for the Bohmian trajectories as the current density divided by the probability density. A particular example will be 
developed for particles with spin in Sec. |V.4|and chapter 8. 
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V.2.2. Many-particle quantum Hamilton- J acobi equation 

Alternatively, we can obtain Bohmian mechanics from a quantum Hamilton-Jacobi equation. We start by introduc- 
ing the polar form of the many-particle wave function ip{x, t) — R{x, i)e*'^*^^'*^/'' into the (nonrelativistic) many-particle 
Schrodinger equation, Eq. (52). Then, after a quite simple manipulation, one obtains from the imaginary part: 



dt 



N 



^ d f 1 dS(x,t) X , 



k=l 



dxk Xm dxk 



where we recognize the velocity of the fc-th particle as: 



Vk{x,t) 



1 dS{x,t) 
m dxk 



(57) 



(58) 



Equations ((56l) and (([58])) are identical. This is shown in problem |6] for a ID system, but it can straightforwardly 
be generalized to N dimensions. The real part of the Schrodinger equation leads to a many-particle version of the 
quantum Hamilton-Jacobi equation: 



dt 



k=l 



2m \ dxk 



(59) 



where we have defined the quantum potential as: 

Q{x,t)^Y.Qk{S,t) 

with: 

Qk{x,t) 



N 



k=l 



fi^ d^R{x,t)/dxl 



2m 



R{x,t) 



(60) 



(61) 



Again, we can obtain a system of coupled Newton-like equations (one for each component of x) for the many-particle 
Bohmian trajectories by computing the time-derivative of the Bohmian velocity of Eq. (56): 



m- 



d 
dxk 



(i/(f,t)-t-g(f,i)) 



(62) 



x—x[t] 



Both the potential V{x,t) and the quantum potential Q{x,t) introduce correlations between particles. All physical 
interactions are essentially expressed as correlations (relations) between the degrees of freedom. In the following, we 
will discuss some important differences between classical and quantum correlations. 



V.3. Factorizability, entanglement, and correlations 

There are important differences between the correlations introduced by the classical potential V{x, t) and the 
quantum potential Q{x, t): 

1. In general, the term V^(af, t) decreases with the distance between particles. A simple example is the Coulomb 



interaction. However, Q(x,t) depends only on the shape of the wave function, not on its value (see Eq. (61)) 
Thus, the quantum potential that we would obtain from Eq. (61) using either R{x,t) or aR{x,t) would be 
exactly the same, even when a — > 0. Thus, the quantum potential can produce a significative interaction 
between two particles, even if they are very far apart rj 



These highly nonclassical features of the quantum potential is what led Bohm to argue that the quantum potential interchanges 
information between systems 
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2. In general, the term V{x,t) produces classical correlations between different particles. The particular depen- 
dence of V{x, t) in all variables xi, . . . , xn imposes a restriction on the speed of the interaction. The variations 
of Xi can only affect Xj after a time larger than la;^ — Xj\/c, c being the speed of light. For example, the 
relation between the particles positions in the electromagnetic interactions ensures that there is no superluminal 
influence between particles. However, such a restriction is not present in the quantum potential. Thus, very far 
particles have an instantaneous (nonlocal) interaction between them. 

The quantum potential is at the origin of all quantum correlations, that is, entanglement, that can imply (nonlocal) 
faster-than-light interactions when two distant particles are involved. As mentioned in Sec. |III.5[ this "spooky action 
at a distance" is what bothered Einstein about Bohmian mechanics (and quantum mechanics, in general). In 1964, 
Bell elaborated his famous theorem that established clear experimentally testable mathematical inequalities that 
would be fulfilled by local theories but would be violated by nonlocal (quantum) ones |45j . All experimental results 
obtained so far confirm that Bell's inequalities are violated. Therefore, contrarily to Einstein's belief, we have to accept 
the real existence, in nature, of faster-than-light causationj^ Entanglement is an intrinsic correlation in quantum 
mechanics (whose complexity and potentialities eventually come from the fact that a TV-particle wave function lives in 
a A^-dimensional configuration space) and is at the core of quantum information science, which makes teleportation, 
quantum communication, quantum cryptography, and quantum computing possible. 

To improve our understanding of correlations, let us discuss under which conditions we cannot expect (classical or 
quantum) correlations between N particles. Let us focus our attention on a wave function 'ip{x, t) that can be written 
as a product of single-particle wave functions associated with each of the particles: 

N 



^ix,t) = Y[iPki^k.t) (63) 



k=l 



We call such a wave function "factorizable." Equation ( 63 ) expresses the physical independence of the N particles (even 
though the wave functions ipk{xk,t) may overlap). According to Born's statistical interpretation of wave functions, 
the absolute square of a wave function is the probability density of the quantum particle; thus the quantum wave 



function, Eq. (63 1, corresponds to a system without (classical or quantum) correlations between particles. It is well 
known that such a solution occurs when the potential in Eq. (52) can be written as V{x) = Vk{xk)- Next, let us 
reformulate this result with Bohmian trajectories. 



We realize that the phase and modulus of Eq. ( 63 1 are given by: 

N 



S{x,t) = J2Sk{xk,t) (64) 

and 

R{x,t) = Y[Rk{xk,t) (65) 



N 



fc=l 

by defining 

Mxk,t) ^ Rkixk,t)e'S''^-'-'^^'' (66) 
In this case, the many-particle quantum potential can be written as: 

2m Rk[Xk,t) 



Then, we can easily deduce from Eq. (591 that each ^pk{xk,t) is a solution of the following single-particle quantum 
Hamilton- Jacobi equation: 



dSk{xk,t) ^ J_ f 9Sk{xk,t) '^ ^ 
dt 2m \ dxk 



+ V{xk,t) + Qk{xk,t) = (68) 



We insist that the experimental violation of Bell's inequalities gives direct support not only to the Copenhagen interpretation but also 
to the Bohmian one, since the latter is also a nonlocal theory because of the presence of Q{x, t). 
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In addition, each {xk , t) satisfies a conservation law: 



dRl{xk,t) d (I dSk{xk,t) 



at '8x,yrn 8x, ^^(--^)j=0 (^9) 
From these equations we can easily deduce an independent guiding equation for each /c-particle: 

showing the absence of correlations between (the different components of the many-particle) Bohmian trajectory. 

V.4. Spin and identical particles 

Elementary particles, such as electrons or quarks, have spin, an internal (discrete) degree of freedom that can 
influence their quantum dynamics in a nontrivial manner. In this section we will briefly explain how to extend 
(nonrelativistic) Bohmian mechanics to include spin. 

V.4-1- Single-particle with s = 1/2 

In the orthodox formulation of quantum mechanics, the state of a single particle with spin s is described by a 2s + 1 
component vector of wave functions with: 

*i(f,i) 

^{r,t)=\ : I (71) 

*2s+i(r,t) 

where r = (a;i,X2,a;3) represents the 3D position of the particle. The time evolution of this state is no longer governed 
by the Schrodinger equation but by more involved wave equations such as the Pauli equation for s = 1/2 [571 1771 ITS] . 

The strategy to find Bohmian trajectories for spin particles will be the following: we will look for a continuity 
equation of the probability density and define the Bohmian velocity as the current density divided by the probability 
density. The idea is quite simple, but the mathematical development can be much more complicated. Therefore, we 
will focus on a particular example for a single-spin s = 1/2 charged particle, whose vectorial wave function, 

has two components and it is called a spinor. In this particular case, the spin along a particular direction takes only 
two possible values, referred to as spin-up ('f) and spin-down (\.) states. The time evolution is given by the Pauli 
equation [57] : 



1 / ^ - \2 

— a-[-ihV-qA{f,t)] +V{f,t) 
im V / 



2m 

where a = (cti, a2, cs) is a vector containing the Pauli matrices [77]: 



^ir,t) (73) 



- - ( ■ o') (^^) 

= (; _°i) m 

and V{'f,t)/q and A{r,t) = {Ai{f,t), A2{'r,t), A^ifjt)) are, respectively, the electromagnetic scalar and vector poten- 
tial. 
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From the Pauli equation we can obtain a continuity equation: 

VJ(f,t) = (77) 



dpir,t) 
dt 



with the probabihty and current densities defined as: 

p{f,t) = -^^{f,t)--^{r,t) (78) 

J(f,t) = U^r^^t) ■ V*t(f,i) _ #t(^-^t) . v#(f,t)) - ^^^p{f,t), (79) 
2m \ / mc 

where "if^ denotes the conjugate transpose of 5'. 
Defining the Bohmian velocity as: 

=^ (BO) 

one can develop Bohmian trajectories for spin particlesj^ It is out of the scope of this book to further develop 
Bohmian mechanics for the Pauli equation. In any case, practical examples of the particle trajectories for spin-1/2 
particles can be found in Refs. [BH [7H] . 
The Pauli equation can be rewritten as: 



at 



(-^W - qA{r, t)) ^-^d- B(f, t) + V{f, t) 
2m \ / 2m 



^(r,i) (81) 



B{r,t) — V X A{f,t) being the magnetic field. Note that the only term that can transfer population between spin 
states is the one with B (the so-called Stern-Gerlach term) . If this term can be neglected then each spin component 
of Eq. ( |8ip reduces to the familiar Schrodinger equation. Therefore, in the nonrelativistic regime, the spin of the 
particle only affects the motion of the trajectories in the presence of a magnetic field. 

On the other hand we can also consider the case where the Hamiltonian is separable into a part depending only on 
the particle position and a part depending only on the spin. Then if the initial state is: 

^(f,0) = ^(r,0)('"^S^n (82) 



"i(0) 

at later times one can write the spinor evolution in a simpler form as: 



^ (f, t) = ^(f, t) ) ^ ^(f, t)x{t) (83) 

where ^p{f, t) depends only on the coordinate of the particles and the function x(t) depends only on its spin. We call 
the former the coordinate or orbital wave function and the latter the spin wave function. In this case, the evolution 
of x(t) and i^{r, t) is independent. If we are not interested in the actual spin of the particles, we can just consider the 
dynamics of the coordinate function ■(/'(r, t), which will be determined by the Schrodinger equation. 

V.4-2. Many-particle system with, s = 1/2 particles 

For a nonrelativistic system of many particles, where the orbital and spin contributions of the Hamiltonian are 
initially separable]^ the many-particle wave function will be also separable at all times, that is: 

4'(ri,...fAr,t) = V(^i,---nv,t)x(t) (84) 



Alternatively, one could also consider the degree of freedom of spin as an additional three-angle variable {a, f}, 7} of the wave function 
(such that R{f,a,l3,'y,t) and S{f,a,l3,'y,t)) and look for the equations of motion of the trajectories of the positions and also of the 
trajectories of the angles (see chapter 10 in Ref. 12011. 
^'^ For example, in the absence of a magnetic field such that the initial wave function can be written as a product of an orbital part, 
4'{ri,r2, ■ ■ -TNtO), and a spin part, x(0)- 
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where Tp{fi, . . .ff^[,t) is the time evolution of the orbital state, the solution of the many-particle Schrodinger equation, 
and is the time evolution of the spin part of the system state: 



' atiT2...t]v(*) \ 
"tiT2..4«(*) 



a2{t) 
\ aw{t) J 



(85) 



W = (2s + 1) is the number of possible combinations of spin projections in one direction for all the particles in the 
system 



Notice that even if the initial many-particle wave function is not separable into orbital and spin parts, that is: 



*(fi,...fA,,0) 



■02(^*1, 



,rN,0)a2 



(86) 



1, . . . , W), is still a solution of the many-particle 



each ipj{ri^ ... ,rN,t), corresponding to a spin combination {j 
Schrodinger equation. 

In many practical situations, we are only interested in the evolution of the orbital wave function through the 
many-particle Schrodinger equation. However, even in this case, there is a peculiar dependence of the dynamics 
of the particles on the total spin. As explained in many textbooks, there is a pure quantum interaction between 
identical particles, named "exchange interaction" [771 [T^. This interaction is not classical, and we cannot find a 
term in the potential energies of the many-particle Hamiltonian of the Schrodinger equation that accounts for it. 
Alternatively, this new interaction is introduced in the "shape" of the global wave function, through the requirement 
of a particular symmetry of the global wave function. We say that a many-particle wave function is anti-symmetric 
when the interchange of the position and spin degrees of freedom associated to two identical fermions (e.g., electrons) 
results only in a change of the sign of the global wave function, 5'(ri, r2, ■ • ■ , r/v, 0- Analogously, we say that the 
many-partcile wave function is symmetric when it remains unchanged after the interchange of the degrees of freedom 
of two identical bosons. 

The crucial point why we cannot forget about spin when dealing with many-particle Hamiltonians is that for 
separable wave functions, for example, Eq. (84), the symmetry of the orbital part depends on the symmetry of the 
spin part. For example, since the total wave function must be antisymmetric, two electrons can have a symmetric 
orbital part if the spin part is antisymmetric, and vice versa. In general, note that the wave function is not separable, 
and it makes no sense to talk about the symmetry of the orbital and spin parts alone but only of the symmetry of 
the total wave function. We see that the exchange interaction induces correlations between particles by imposing 
symmetries (or shapes) to their many-particle wave function. 

A standard claim in many quantum mechanics textbooks is that identical particles, for example, two electrons with 
an antisymmetric wave function, are indistinguishable. It is affirmed that if the particles would have trajectories, 
they would automatically be distinguishable. In Bohmian mechanics, even with the symmetrization postulate, the 
adjective "indistinguishable" is inappropriate because one can label one particle's trajectory fi[t\ and the other f2[t\ 
and thus distinguish them perfectly. 

Bohmian trajectories can actually help us improve our understanding of the symmetrization postulate. Let us 
assume a two-electron system with an antisymmetric orbital wave function ^/'(ri, We assume that an electron 
labeled 1 with the initial position fi[Q] evolves into fi[t] and an another electron labeled 2 evolves from f2[0] to 
r^lt]. Then, it can be easily understood that ri'[0] — f2[0] evolves into fl'lt] — r2[i] and f2'[0] = ri[0] evolves 
into r^'lt] = fi[t]. We use primes to notice that ri[i] and fi'[t] correspond to trajectories of the same particle with 
different initial positions {fl[Q] or r^[0]). This result follows from the symmetry of the velocity, that is, the symmetry 
of the current density and modulus, when positions of the two electrons are interchanged. In summary, Bohmian 
trajectories of identical particles are clearly distinguishable, but the observable (ensemble) results obtained from them 
are indistinguishable when we interchange their initial positions. 



Additionally, one can also look for the projection of the total spin of the system. Both procedures are connected by the Clebsch— Gordan 
coefficients I77| . 
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V.5. Basic postulates for many-particle systems 



The basic postulates of Bohmian mechanics for many-particle systems are exactly the same as those we have 
developed for a single particle but adding the symmetrization postulate. As done in Sec. |IV.6[ we also present a 
minimalist formulation of the postulates. 

In order to simplify our description, we assume the simpler many-particle wave function written as: 



*(fi, . . . ,rW,0 = -0(^1, ■ • ■ ,nv,i) 



/ tlt2 ■■■tNit)\ 





V 







(87) 



The generalization of the present postulates to include the more general wave functions of Eqs. ((84|) or ((86)) follows 
straightforwardly. However, it will imply writing W different wave functions in the postulates. 



FIRST POSTULATE: The dynamics of a many-particle quantum system comprises a wave function, Eq. (81), 
whose orbital part ^(ri,r2T ■ ■ j^N^t) is defined in the configuration space {ri,f2, . . . ,risi} plus time, and a many- 
particle trajectory {r'llt], f2[t]-, ■ ■ ■ ,r/v[^]} that moves continuously under the guidance of the wave function. 

The orbital part of the wave function il^{ri, r2, ■ • • , Hv, is a solution of the many-particle Schrodinger equation: 



ih 



dip{ri, ...,rN,t) 
dt 



N 



V|. + V{ri, ...,rN,t) ■■■,rN,t) 



Each component rk[t] of the many-particle trajectory {fi[t], . . . , fjv[t]} is obtained by time-integrating the particle 
velocity Vk[t] = Ufc(ri, . . . , fjv, t) defined from: 

1^ ^ ,N Jk(fx, ...,rN,t) 

Vk{ri,...,rpf,t) = — — _^ 

|V'(ri, . . .,rN,t)\^ 

where |V'(^i, . . . , rw, i)P = il^^i^i, . . . ,rpf,t)ip*{ri, . . . ,f^N,t) is the square modulus of the wave function and Jk = 
Jk{ri, ■ ■ ■ ,rN,t) is the k-particle current density: 

h 



Jk = i 



2m 



(^ipiri, .., fN ,t)V p^ip* (ri, ..,fN,t) - ■ip*{fi, ..,rN,t)'^p^ip{fi, ..,rN,t)^ 



The initial positions {fw,f2o, . . . ,rNo} o,^d velocities {vio, V20, ■ ■ ■, "Sjvo} have to be specified in order to completely 
determine the many-particle trajectory. 



SECOND POSTULATE (quantum equilibrium hypothesis): The initial positions {r[[to], '^[^o], • • • , ^atI^o]} o-nd 
velocities {vWt^)],V2[tG]^ ■ ■ ^jvl^n]} ''Z particular many-particle j trajectory cannot be known with certainty. When 
the experiment is repeated many times, these initial positions of an ensemble of trajectories associated with the same 
'(/'(ri, . . . , r/v, t) have to be generated so that the number of trajectories of the ensemble between [ri, . . . ,riq) and 
(ri -t- (iri , . . . , rjv + c^r/v) at the initial time to is proportional to i?^(r^[to], . . . , r^[to]) = |'0('^[^o]i ■ • • j^[^o]:io)P- 
The initial velocity is determined by 'y|,[io] = JkiT^i [to], ■ ■ ■ , ^jv[^o], ^o)/|'0(^i['^o], • ■ • i ^^[^o], ^o)P ■ 

The condition on the initial position can be written mathematically as: 

■ ' " M 51 n ^^^^ ^k[to]) for t = ta 



3 = 1 fc=l 



Notice the presence of two indices, the j = 1, . . . , M for the infinite ensemble and the k — 1, 
of freedom. 



, N for the N degrees 



THIRD POSTULATE (symmetrization postulate of quantum mechanics): If the variables ri,^i andrj,^j refer to 
two identical particles of the ensemble, then the wave function, Eq. (81), is either symmetric: 



il;{.,f,,.,rj,.,t) 







V 



t) 







(89) 
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if the particles are bosons (every particle with an integer spin 0, 1, 2, . . . is a boson), or antisymmetric: 



ip{.,fi,.,fj,.,t) 







V 







-ip{.,fj,.,n,.,t) 







V 







(90) 



if the particles are fermions (every particle with a half-odd spin 1/2, 3/2, . . . is a fermion). In Eq. 1^8!^ and Eq. ( 90) 
it is understood that all other degrees of freedom of the other particles remain unchangedr^ 

Again, no postulate about measurement is needed, since in Bohmian mechanics measurement is treated as a 
particular case of the interaction between particles (see Sec. VI ) . 



V.6. The conditional wave function: 



many-particle Bohmian trajectories without the many-particle wave 
function 



In section 1.2.4, when discussing the similarities and differences between classical and quantum mechanics, we 
mentioned quantum wholeness (the dependence of each individual Bohmian trajectory on the rest of trajectories 
through the wave function) as a fundamental difference. This concept acquires an even more dramatic meaning 
when dealing with many-body systems. A particular Bohmian trajectory of the fc-particle f|,[i] depends on the rest 
of infinite j = 1,...,M possible trajectories (with different initial position) of all k — 1,...,N particles. This 
property is mathematically manifested by the fact that the many-particle wave function is defined in a iV-dimensional 
configuration space rather than in the usual real 3D space. In a 3D real space one can easily define whether or 
not two particles are far away. In general, the potential profile that determines the interaction between these two 
particles decreases with the distance. However, two distant particles that share a common region of the iV-dimensional 
configuration space where the wave function is different from zero have an interaction (entanglement) independently 
of their distance. 

As discussed in section 1.3.1, the need of computing the many-particle wave function in the A'^-dimensional config- 
uration space is the origin of the quantum many-body problem. Among the large list of approximations present in 
the literature to tackle this problem, Bohmian mechanics provides a natural, original and mainly unexplored solution 
through the use of the conditional wave function ^80, . Due to the Bohmian dual description of a quantum system as 
particles and waves, one can reduce the complexity of a many-particle wave function through the substitution of some 
of the N degrees of freedom by its corresponding Bohmian trajectories. The new wave function with a reduced number 
of degrees of freedom is named conditional wave function. For a discussion about the fundamental implications of 
working with conditional wave functions instead of full many-body wave functions, see Ref [75] . We will present now 
some preliminary work on the computational abilities of such conditional wave functions 

Hereafter, for simplicity, we return to a definition of the degrees of freedom of the many-particle quantum sys- 
tem in terms of \ xi, X2, ■ ■ ■ ,X]\[} instead of {ri,f2j • ■ • ,rw}- The direct solution of the many-particle Schrodinger 
equation, Eq. (52), is intractable numerically. As we have discussed, the main reason why such an equation is prac- 



tically unsolvable is because it looks for a solution within the configuration space associated to all particles, that is, 
{xi, . . . , Xa, . . . , xn}. On the contrary, classical equations deal with solutions in a much smaller configuration space, 
{xi [t], . . . , Xi^i[t] , Xa, Xii [t], . . . , X]^ [t]}, where all other trajectories are known parameters, except Xa- In other words, 
the Newton solution of Xa[t] just needs the spatial dependence oiV{xi[t], ... ,Xa, ■■■ ,XN[t]) on the variable Xa- In this 
section, we summarize the formalism developed by Oriols in a recent work |81j to compute many-particle Bohmian 
trajectories without knowing the many-particle wave function. It is a clear example of how Bohmian mechanics can 
be a powerful computational tool. 

The main idea behind the work developed in [8r is that any Bohmian trajectory Xa[t] that is computed from 
the many-particle wave function ^'(x, i) solution of Eq. ( |52[ ) can be alternatively computed from a much simpler 
single-particle conditional wave function (f)a{xa,t) = '^{xa,xi,[t],t). Here, we use the notation x — {xa,Xb} with 
Xb = {xi, . . ■ ,Xa~i,Xa+i, ■ . ■ ,xn} for particle positions and xi,[t] — {xi[t], . . . ,Xa-i[t],Xa+i[t], . . . ,xi\[[t]} for Bohmian 
trajectories. 

It is quite trivial to demonstrate the ability of (paixajt) in reproducing Xa[t]. By construction, when we use a 
polar form 4>a{xa,t) = ra{xa,t) q''- Sa(xa,t)/h ^ angle Sa{xa,t) is identical to the angle S{xa,Xb,t) of '^/{xa,Xb,t) 



■^^ This simple spin vector wave function is clearly symmetric so that the orbital wave function has to be either symmetric or antisymmetric. 
For general wave functions, such as the ones in Eqs. (l |84| ) and (| |86| l), this postulate implies much more complicated restrictions on the 
possible orbital and spin wave functions. 
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evaluated at {xa, xi,[t]}. Therefore, since the velocity of the trajectory Xa[t] is computed from the spatial dependence 
of S{xa-, Xb,t) on Xa whcu all other positions are fixed at Xb = Xb[t], the same velocity will be obtained from the 
spatial dependence of Sa{xajt). Interestingly, (f)a{xa,t) is solution of a single-particle (pseudo) Schrodinger equation 
because it depends only on time t and position Xa- Next, our effort will be focused on determining such single-particle 
equations. 



V.6.1. Single-particle pseudo-Schrodmger equation for many-particle systems 



First of all, we show that any (totally arbitrary) single- valued complex function (f)a{xa,t), which has a well-defined 
second-order spatial derivative and first-order temporal derivative, can be obtained from a Schrodinger-like equation 
when the following potential W{xa,t) is used: 



W{Xa,t)^ 



ih- 



dt 



2m a; 



4'a{Xa,t) 



(91) 



For an arbitrary (complex) function, the potential energy W{xa,t) can be complex, too. In fact, we are interested in 
rewriting W{xa^ t) in terms of the polar form of the wave function (pai^a, t) — fa{xa^ t) e* ^"(^i. *)/''. We obtain for the 
real part: 



Real[W{xa,t)] = - \ 



1 / dSa{Xa,t) 



2 m 



dxa 



d'^ra{Xa,t) ^ dSa{Xa,t) 

2mra(xa,t) d'^Xa dt 



From Eq. (91 ), we do also obtain for the imaginary part: 

^ ( drl{xa,t) d / rl{xa,t) dSa{Xa,t) 



Imag[W{xa,t)] 



2r 



l{Xa,t)\ dt dXa 



dXa 



(92) 



(93) 



It can be easily verified that Imag[W{xa,t)] = when the single-particle wave function preserves the norm. 

Finally, we are interested in using the above expressions when the single-particle wave function is the one mentioned 
in the introduction of this subsection, (j>a{xa,t) — '^{xa,Xb[t],t). In particular, we use ra{xa,t) = R{xa,Xb[t],t) and 
Saixa,t) = S{xa,Xb[t],t). Then, we realize that Eq. (91 1 transforms into: 



ih- 



t'aiXg.t) 

dt 



2m d'^Xa 



[t],t) -\- GaiXa,Xb[t],t) -\-iJa{Xa,Xb[t],t) ](t)aiXa,t) 



where we have defined: 



N 



Ga{x,t) = Ub{xb,t)+ 



1 fdS{x,t) 



2 m \ dxi 



+ Qk{x,t) 



dS{x,t) 
dxk 



Vk{x[t],t) 



(94) 



(95) 



with Qk{x, t) as the quantum potential energy defined in Eq. (61 ). The terms Ua{xa, Xb[t],t) and Ub{xb, t) are defined 
from the many-particle potential that appears in the original many-particle potential energy in Eq. (52) as: 



In addition, we have defined: 



N 



Ja{x,t)^ ^ 



V{x,t) = Ua{Xa,Xb,t) + Ub{Xb,t) 



h dR^(x,t) d fR^(x,t)dS(x,t) 



k—l;k^a 

In order to obtain these expressions, one has to carefully evaluate: 



dxk 



dS{Xa,Xb[t],t) _ f dS{Xa,Xb,t) 



dt 



dt 



N 



^ dS{Xa,Xb[tlt) 

+ 2^ vk{x[t],t) 



dxk 



(96) 



(97) 



(98) 



and use Eq. (59) evaluated at {a;a,Xb[i]}. Identically for dR {xa, Xb[t],t)/dt with Eq. (57). 

The single-particle pseudo-Schrodinger equation, Eq. (94), is the main result of this algorithm. Let us discuss the 
role of each potential term in Eq. (|94|: 
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1. The term Uaixa, Xb[t],t) is a real- valued potential whose explicit dependence on the positions is known. It has 
to be evaluated from the particular Bohmian position of all particles except Xa[t\. 

2. The term Ga{xa, a;b[t], t) is a real- valued potential whose explicit dependence on the positions is unknown (unless 
we know the many-particle wave function) and needs some educated guess. It takes into account, for example, 
the exchange interaction between particles. 

3. The term iJa{xa,Xh[t],t) is an imaginary- valued potential whose explicit dependence on the positions is also 
unknown and needs some educated guess. It takes into account that the norm of (j3a{xa,t) is not directly 
conserved (the norm of the many-particle wave function 'i!{xa, xi,,t) is conserved in the entire configuration 
space, but this is not true for ^!(xa,Xb[i\,t) in the Xa space). 



By using Eq. (94) for each particle, a;a[i] for a = 1, . . . , iV, we obtain a system of N coupled single-particle pseudo- 



Schrodinger equations that is able to compute many-particle Bohmian trajectories without knowing the many-particle 



wave function. The great merit of Eq. (94) is to demonstrate that such a single-particle solution of a many-particle 
problem exists, although we do not know exactly the values of the terms Ga{xa,Xb[t],t) and Ja{xaTXb[t],t). Our 
algorithms have similarities with the original work on density functional theory [721 [73]. The formidable computational 
simplification comes at the price that some terms of the potential of the corresponding single-particle Schrodinger 
equations are unknown, that is, the exchange-correlation functional in density functional theory [721 173] and, here. 



the terms in Eqs. ((95)) and ((97)) 



V.6.2. Example: Application in factorizable many-particle systems 



Let's start by discussing which will be the solution <j)a{xa,t) of Eq. (94 1 when the many-particle wave function 
^{xa,Xb,t) is factorizable (i.e., when it can we written as Eq. (63)). Then, Eq. (68) is valid for each summand of 
Ga{xa, Xb[t\,t) in Eq. (94), so it can be written in a compact way as: 



N 



Ga{Xa,Xb[t\,t) 



E 



dSk{xk[t],t) 
dt 



(99) 



Then, a (real) time-dependent term (without any spatial dependence) appears in the potential of Eq. (94). It can be 
easily demonstrated that such a term introduces the following time-dependent contribution j3a{t) into the phase of 

MXajt)- 



N 

Pa{t) = - E / 



k—l\k^a * 



dSk{xk[t'lt' 
dt' 



N 



-dt' = - Skixk[t],t) 

k—l\k^a 



Identically, Eq. (69 1 is valid for each term of Ja{xa,Xb[t],t)^ so it can be written as: 



N 



k—l:k^a 

Using ln(a h) — ln(a) + ln(6), we obtain a contribution aa{t) into the phase of (pa{xa^ t): 



(100) 



(101) 



aait) = -h f ^\nl II Rkixk[t'],t')] dt' 



= -h\nl II Rkixk[t],t)\ (102) 

yk=lM^a j 

Finally, we use that a pure (real or imaginary) time-dependent potential (without spatial dependence) in a Schrodinger- 
like equation does only introduce a pure (imaginary or real) time-dependent global phase. Thus, we obtain: 

'i fia{t) - aa{ty 



i>a{Xait) = exp 



^a(Xa-,t) = ili{xi[t],t) . ..^paiXa,t) ■ ■ ■ ^pN {x N[t], t) 



(103) 



which is, certainly, the expected result. Each term ipkixkjt) is a single-particle wave function, whose evolution is 
found from Eq. (30) after appropriately defining the initial wave packet at time t — 0. Notice the difference between 
i^aixajt) and (j>a{xa,t). The former is a single particle wave packet, while the latter has an additional time-dependent 
function multiplying ipa{xa,t). 
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FIG. 5. For 4 interacting particles without exchange interaction, the present algorithm needs A'' = 4 single-particle wave 
functions, ^pa{^a,t). The subindex a of the wave function is associated with the potential Ua{xa, xi,[t],t) and the initial wave 
packet. Each wave function %l)a(xa,t) determines the a-Bohmian trajectories Xa[t\- See also Color Insert. 



V.6.3. Example: Application in interacting many-particle systems without exchange interaction 



Up to here, we have demonstrated that a Bohmian solution to the many-particle problem exists in terms of a system 
of coupled single-particle Schrodinger equations (see Eq. (94)). The significant computational simplification comes at 
the price that the terms Ga(a;o, a;{,[<], i) and , [i] , t ) of the corresponding single-particle Schrodinger equations 

are unknown. It is in this sense that we mentioned that these algorithms have similarities with density functional 
theory [72l [73] . In this subsection we provide a simple approximation for a system of N electrons with Coulomb 
interaction but without exchange interaction. Later, we will include the exchange interaction. 

As mentioned, the solution of Eq. (94) needs educated guesses for the terms Ga{xa,Xi,[t\,t)^ Eq. ([95]), and 



Ja{xa,Xb\t\,t)^ Eq. (97). Since no exchange interaction is considered, we assume that the origin of the correlations 
between the a electron and the rest is mainly contained in the term Ua{xa,Xf,[t],t). We develop a Taylor expansion 
of the other two terms, Eqs. (([95])) and ((97)), in the variable Xa around the point Xa[t]: 



G t)=Ga{Xa[tlt) + 



dGa{Xa,t) 
dXa 



Xa—3^a[i] 



(Xa - Xa[t]) + 



(104) 



and 



Ja{Xa,t) = Ja{Xa[t\,t) + 



dJa{Xa,t) 



dx„ 



{Xa - Xa[t\) + . . 



(105) 



The simplest approximation is just a zero-order Taylor term Ga{xa,t) ~ Ga{xa[t\,t) and Ja{xa,t) ~ Ja{xa[t],t). 
The conditional wave function (paixajt) solution of Eq. (94) can be constructed in two steps. First is solving Eq. 



(94) without considering the purely time-dependent potential terms, Ga{xa[i\,t) and Ja{xa[t\,t), to find ipaixajt): 



ih 



dtjja{Xa,t) 

dt 



2m dx'i 



[t\,t) ) ^pa(Xa, t) 



(106) 



where the term Ua{xa,Xb[t\,t) is defined in Eq. (96). Second is multiplying the wave function Tpai^Xa^t) by time- 



dependent (real or imaginary) values (without any spatial dependence) for the final solution: 



(t)aiXa,t) W i}a{Xa,t) exp(zQ(t)) 



(107) 



with Za{t) = i/3a{t)/h — aa{t)h defined according to Eq. (103). Again, we have used the well-known result that a 



pure (real or imaginary) time-dependent potential term (without spatial dependence) added into the Hamiltonian of 
the Schrodinger-like equation does only introduce a pure (imaginary or real) time-dependent global phase, that is, 
exp(z£i(t)). Since the velocity of Bohmian trajectories does not depend on these pure time- dependent terms exjp{za(t)), 
we do not have to compute /3a (i) and aa{t) explicitly. 
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FIG. 6. Two-particle interacting Bohmian trajectories in a triple barrier tunneling scenario computed from our ID approach 
(symbols) and from exact 2D results (solid line) for three different lateral areas that modifies the Coulomb interaction. Reprinted 
with permission from [81]. Copyright 2007 American Physical Society. 



For example, for a system of iV = 4 interacti ng el ectrons, as depicted in the scheme of Fig. [5] we need to solve 
N — A Schrodinger-like equations defined by Eq. ( 106 1 to find = 4 wave functions ipa{xa, t). Many-particle Bohmian 
trajectories are computed from the Bohmian velocity Eq. (35) using ipai^a, t) to compute the current and the square 
modulus. The term Ua{xa,X}j[t],t) can be solved from a Poisson equation (if we deal with Coulomb interaction) and 
it introduces correlations between particles. The initial wave function -tjjaixajt ~ 0) has to be specified. For example, 
if we assume that the initial wave function is defined in a region of the space where the many-particle wave function 
is factorizable, we can use ipa{xa, 0) = V'a(a^a, 0). 

We will show now an example of the goodness of this simple approximation. We control the strength of the Coulomb 
interaction between two electrons by changing the lateral area of a 3D electron device, as described in 'ST . Figure 
[6] shows the excellent agreement between exact Bohmian trajectories and those computed within our algorithm. In 
Fig. [6^, the lateral area is so large that it makes the Coulomb interaction quite negligible. The first electron is 
transmitted, while the second is reflected. However, as shown in Fig. [6j;, the smaller lateral area provides strong 
Coulomb interaction between the electrons, and the second is finally transmitted because of the presence of the first 
one in the barrier region. 

Certainly, we have used the simplest approximation in Eqs. ((104)) and ((105)). Any possible improvement of this 
approximation will imply an even better agreement between this algorithm and the exact computation. 



V.6.4- Example: Application in interacting many-particle systems with exchange interaction 



Now, we generalize the previous result to an arbitrary system with Coulomb and exchange interactions. For 
simplicity, we will consider only the wave function defined in Eq. (87), that is, only the symmetry of the orbital wave 
function is considered. A simple Taylor approximation for the terms Ga{xa,t) and Ja{xa,t), as we assume in Eq. 
(104 1 and Eq. (1051, is not valid here. The reason can be easily understood, for example, for a system of (identical) 
electrons. Due to the Pauli exclusion principle, the modulus of the wave function tends to zero, R{xa, Xb[t]^t) — 0, 
every position where Xa — > Xk[t\. Thus, the term Ga{xa-, Xb[t],t) has asymptotes at Xa — Xk[t] that repel other 
electrons, that is, Ga{xa, xi,[t],t) — > ±oo at the same particular positions x^, invalidating the approximation Eq. 
(104 1. For the same reason, Ja{xa,Xb[t],t) — >■ ±00. 

Here, the strategy is assuming that an antisymmetric wave function, '^{x, t), can be constructed from permutations 
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FIG. 7. When Coulomb plus exchange interactions are considered, a 4 - 4 matrix of the single-particle wave functions, ^a,h{xa,t) 
(with a = 1, ... ,4 for the potential energies and j — 1, ... ,4 for the initial wave packets) is needed to compute (through a 
Slater determinant) the wave function associated with each Bohmian trajectory. See also Color Insert. 



of a many-particle wave function without any symmetry \E'no-sym(a?, i): 

= C'^^no-sym{Xp(l),,Xp(^i)^,. . . ,Xpi^i)^,t) S {p{l)) (108) 
= 1 



The subindex "no-sym" reminds that such wave functions have no spatial symmetric or antisymmetric property. The 
constant C is a normalization constant that will become irrelevant for the computation of the Bohmian velocity. The 
sum is over all A^! permutations p{l) = {p(l)^,p{l)^, . . . ,p{l)j^}, and s {p{l)) = ±1 is the sign of the permutations. 
Then, each wave function ^'no-sym(a;p(;)j , a;p(;)2 , . . . , X p(i)^ , t) eval uate d at £" = {a;Q,Xb[i]}, can be computed following 
the previous (no-exchange) algorithm, that is Eqs. (( |106| ) and ((107)). In particular: 

*no-sym(a;p(/)Jt], . . .,Xp(i)-,. . . , Xp(^i) ^^[t], t) = 1pa,p(l) , {Xa, t) exp {Za,p(l){t)) (109) 

where the permutation p(l) gives a;p(z)j = Xa, that is, p{l)j = a. Now, we have to use two labels in the subindex 

of tpa,h{xa,t) to Specify the solution of Eq. (106). The first label a accounts for the degree of freedom, that is, the 
particular trajectory, that we are computing and it also fixes the potential energy Ua{xa,Xb[t],t) in Eq. (106). The 
second label h fixes the initial wave function that we will consider. If the initial many-particle wave function can be 
defined, in a region without interactions, as: 

*no-sym(z^l, • • • ,Xn,0) = ll^i{xi,0) . . . V'JV (a;Ar , 0) (110) 

then ipa,p{i)j [xa, 0) = ''l'p(i)j {xa^ 0). In other words, identical initial wave functions ^'^^ (xi, <o) = '^k.j{xk^ to) can evolve 
differently when usin g Ui{ xi, Xi[t], t) or Uk{xk,Xk[t\,t)- 

Finally, using Eq. (107), the many-particle wave function (p(xa,t) = '^{xa,Xii[t\,t) can be written as: 

N\ 

(l){Xa,t) = ■i>{XaXb[t],t) = C^-iija,p(l)^(Xa,t) exp {Za^p(l){t)) S {p{l)) (111) 

1 = 1 
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FIG. 8. Two-particle Bohmian trajectories with Coulomb and exchange interactions computed from our ID approach (symbols) 
and from 2D exact results (solid line) for three lateral areas. Reprinted with permission from |81] . Copyright 2007 American 
Physical Society. 



Here, the angles Za^p(i){t) are relevant and cannot be ignored (only a global phase can be ignored). We emphasize 



that j is selected according to the condition p{l)j 



which depends on the index I of the permutations. 



In fact, we can guarantee the expected symmetry of the quantum wave function by defining appropriately the 
phases. It is quite simple to realize that the following definition of the angles: 



exp(za,p(;)(0) 



N 

n 



4- 



a,p{l), 



(112) 



accomplishes the symmetry requirements. Finally, putting Eq. ( 112 1 into Eq. ( 111 ), we obtain the final wave function: 

NI. 

<P{Xa,t) = C ^ ^l,p(l), ■ ■ ■ *a,p(/)„ {Xa. t) . . . ^ N ,p(l) ^ (xjv [t] , i)s(p(0) (US) 



1 = 1 



,4 



This Eq. ( 113 ) can be simply computed from a (Slater) determinant of the N ■ N matrix, ipa.hixa, t) with a = 1 
and j — 1, ... ,4 for the initial wave packets, which are defined in Fig. [7| 

The reader can wonder why we need N'^ different wave functions and not iV!. The reason is because the term 



Ua{Xa,Xb[t],t) in Eq. (106) does not change when we interchange the positions of two of the trajectories in Xb[t\ (only 
the interchange between Xa[t] and a trajectory of Xb[t\ becomes relevant). Figure |8] shows the excellent agreement 
between the exact Bohmian trajectories and those computed from this second algorithm. See details in |81j . 

In order to clarify the meaning of indistinguishable particles in the Bohmian language, in Fig. |9j we simulate 
Bohmian trajectories with the initial positions and with interchanged initial positions. By our own construction of 
the wave function of Eq. (113), the velocity has symmetry when we interchange particles, as discussed in Sec. V.4 



Therefore, although we can distinguish each electron by its own trajectories, one cannot discern between them in the 
final observable results. For example, the number of events where the first trajectory is transmitted and the other 
reflected is exactly identical to the events where the first is reflected, while the second transmitted. On the contrary, 
as seen in Fig. |9j this symmetry is broken without the exchange interaction algorithm. 



VI. BOHMIAN EXPLANATION OF THE MEASUREMENT PROCESS 



The Bohmian explanation of a quantum measurement is, perhaps, the most attractive (and also ignored) feature 
of the Bohmian explanation of the quantum nature |57l [531 [HlHSlj- Let us start by noting that Bell disliked the 
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FIG. 9. Two-particle Bohmian trajectories with Coulomb and exchange interactions computed from our ID approach (symbols) 
and from 2D exact results (solid line) when initial positions are interchanged. Reprinted with permission from ^81) . Copyright 
2007 American Physical Society. 



word "measurement" [U]- He preferred the word "experiment" because "when it is said that something is measured 
it is difficult not to think of the result as referring to some pre-existing property of the object in question." On 
the contrary, in an experiment, it is natural to think that everything can change, that is, time-evolve, because of 
interactions. 



VI. 1. Differences between orthodox and Bohmian measurements 



According to Bohmian postulates, the evolution of the wave function is uniquely determined by the Schrodinger 
equation with the appropriate Hamiltonian. Surprisingly, this is not true for the orthodox formulation of quantum 
theory. Let us discuss this point in detail in the next two subsections. 



VI. 1.1. The orthodox measurement process 



In the measurement process, one usually separates the quantum system and the measuring apparatus (see Fig. 

The orthodox prediction of some experimental property of the quantum system is described through the use of a 
proper operator G whose eigenvalues give the possible outcomes of the measurement. When we measure a particular 
eigenvalue, the initial wave function is transformed into an eigenfunction of the operator. This is the so-called von 
Neumman (or projective) measurement. Thus, the time evolution of the wave function of a quantum system is 
governed by two (quite) different laws: 

1. The first dynamical evolution is given by the Schrodinger equation. This dynamical law is deterministic in the 
sense that the final wave function of the quantum system is perfectly determined when we know the initial wave 
function and the Hamiltonian of the quantum system. 

2. The second dynamical law is called the collapse of the wave function. The collapse is a process that occurs when 
the wave function interacts with a measuring apparatus. The initial wave function before the measurement is 
substituted by one of the eigenstates of the particular operator G. Contrarily to the dynamical law given by 
the Schrodinger equation, the collapse is not deterministic, since the final wave function is randomly selected 
among the operator's eigenstates. 
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FIG. 10. (a) The Bohmian measurement assumes that the quantum system and the measuring apparatus are exphcitly 
simulated, (b) The orthodox measurement assumes that only the quantum system is explicitly simulated, but the measuring 
apparatus is substituted by a proper operator acting on the wave function of the system. 



The duality in the time evolution of a quantum system in the orthodox interpretation is certainly a persistent 
controversial issue, the so-called measurement problem or collapse of the wave function [47\ I85j . As repeatedly 
stressed by Bell, the orthodox theory is unprofessional because it does not explain with total accuracy which parts 
of the entire quantum system form the the measuring apparatus and which the quantum system itself. It is even not 
clear if the measuring apparatus needs some kind of human activity {with a PhD? [47 ) to be accepted as a proper 
measuring apparatus. The problematical way in which measurement is treated in the orthodox formulation has been 
stressed by Bell [H]: 

The concept of "measurement" becomes so fuzzy on reflection that it is quite surprising to have it appearing 
in physical theory at the most fundamental level. Less surprising perhaps is that mathematicians, who 
need only simple axioms about otherwise undefined objects, have been able to write extensive works on 
quantum measurement theory- which experimental physicists do not find it necessary to read. [. . . ] Does 
not any analysis of measurement require concepts more fundamental than measurement? And should not 
the fundamental theory be about these more fundamental concepts? 

In short. Bell argued that the separation between the quantum system and the measuring apparatus in Fig. |10| 3 
is arbitrary. The encapsulation of the rest of the world (except the quantum system) into a mathematical entity 
called an operator, G, is a very clever trick that allows for straightforward calculations of the results of quantum 
measurements [SSI HI] without considering the rest of the world. 



VI. 1.2. The Bohmian measurement process 



In Bohmian theory, the measurement process is treated just as any other quantum process of interacting particles 
and the previous measurement difficulties of the orthodox interpretation simply disappear. There is no need to 
introduce operators [63 [ I86| , l87]. Here, the entire quantum system is described by a trajectory plus a wave function 
(rather than a wave function alone). The wave function and the trajectory are both associated with the entire system, 
that is, the quantum system plus the measuring apparatus. Then, there is one dynamical law for the evolution of the 
wave function and another for the evolution of the trajectory: 

1. The Schrodinger equation (with the appropriate Hamiltonian of the quantum system plus the measuring ap- 
paratus) determines the time evolution of the wave function independently of whether a measurement process 
takes place or not. 

2. The time evolution of the particle is determined by the time integration of the Bohmian velocity independently 
of whether a measurement process takes place or not. 

For example, imagine that some kind of pointer indicates the measured quantity; then, the particles (degrees of 
freedom) that form the pointer must be present in the Hamiltonianpj In this sense in Bohmian mechanics, a physical 



In modern electronic measuring devices, the pointer could be represented by a seven-segment array of light-emitting diode (LED) 
displays, each one with two possible states, ON and OFF. When electrons are present inside the PN interface of one of the LEDs, a 
radiative transition of the electrons from the conduction to the valence band produces light corresponding to the ON state. The absence 
of electrons is associated with an OFF state. 
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(a) Bohmian measurement explanation in [xs,xa] space 
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(b) Orthodox measurement explanation in [xs] space 



▲ Xs 



4 Xs 



Wave function collapse 




II 




a ga 



FIG. 11. (a) Bohmian measurement in the [a::s,a^A] space: From the nonoverlapping many-particle (system + apparatus) wave 
function, only the Qa part of the wave function where the trajectory is present is needed to compute the evolution of the 
Bohmian system, (b) Orthodox measurement in the [xs] space: The (system) wave function collapses into the Qa part when 
the measurement takes place. See also Color Insert. 

quantum system must be described by a many-particle Hamiltonian. Once the Bohmian trajectories associated with 
the positions of the pointer are known, the value of the measurement is already predicted. We just need knowledge 
of the positions of the pointer particles. Again, Bell provided us with one of his didactic sentences [Hj : 

In physics the only observations we must consider are position observations, if only the positions of 
instrument pointers. It is a great merit of the de Broglie-Bohm picture to force us to consider this fact. If 
you make axioms, rather than definitions and theorems, about the "measurement" of anything else, then 
you commit redundancy and risk inconsistency. 

Therefore, a proper modeling of a Bohmian measurement just needs the explicit consideration of the degrees of 
freedom of the pointer in the many-particle wave function and many-particle Bohmian trajectories that define the 
entire system. See Fig. 11 In particular, we have to introduce into the Hamiltonian the interaction of the particles 
of the pointer with the rest of the particles of the system. The back reaction of the measurement process on the 
wave function is trivially considered. Certainly, a Hamiltonian with or without the measuring apparatus will provide 
a different evolution of the quantum system wave function. For this reason. Bell preferred the word "experiment" 
instead of "measurement." 

The Bohmian and the orthodox explanations of a measurement produce the same probabilistic predictions. However, 
the mathematical implementation of the equations of motion in each case is quite different. The orthodox quantum 
theory requires an operator to describe the effect of the measuring apparatus, but this operator is not needed in the 
Bohmian explanation. 



VI. 2. Theory of the Bohmian measurement process 



Now, we provide a mathematical demonstration of the ideas mentioned earlier. In the standard interpretation of 
quantum theory, the von Neumann (projective) measurement process is defined in a particular quantum region, which 
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is defined through the degrees of freedom xs- See Fig. [TT] 3. The state of the quantum system in this particular region is 
determined through the wave function ipsi^s, t). The process of measuring a particular magnitude is mathematically 
defined through an operator, for example, G, acting on the wave function. The possible outcomes of the measurement 
process correspond to one of the possible eigenvalues g of this operator that satisfy the equation Gipg{xs) = gi'gi^s) 
with ipg{xs), an eigenvector of this operator. The set ipgi^s) forms an orthonormal basis of the Hilbert space of the 
quantum system so that the arbitrary initial wave function can be written as: 

V's(^s,i) (114) 
a 

with Cg{t) a complex value with the only restriction |cg(t)P = 1 in order to ensure that ipsi^s, t) is well normalized. 
When measuring the eigenvalue ga the total wave function rpsixs^t) collapses into Tpg^(xs)- Then, the probability 
of obtaining ga in the measuring apparatus is just Pg^ = \cg^{t)\'^. In order to avoid unnecessary complications, 
hereafter, we have assumed that this basis has no degeneracy. 

In order to understand how this (von Neumann) measurement can be described by means of Bohmian mechanics, 
we repeat that the measurement process in Bohmian mechanics is treated just as any other quantum process, which is 
described by the Schrodinger equation plus a trajectory alone, defined both in the enlarged configuration space that 
includes the quantum system and in the measuring apparatus. There is no privileged status for the observer, that is, 
the measuring apparatus, over the quantum system. This is the reason why some authors describe Bohmian theory 
under the title of "quantum mechanics without observers" |83j . 

Surprisingly, this idea was, somehow, anticipated in the book Quantum Theory, published in 1951 [34 , by Bohm. 
Although in this book Bohm followed an orthodox presentation of quantum mechanics, he also provided a detailed and 
conceptually intricate treatment of the measurement process in chapter 22. First, Bohm demanded in unequivocal 
terms a quantum theoretical treatment of the entire process of measurement within the Schrodinger equation. In 
particular, he argued that a good measuring apparatus will force the wave function to decouple into several nonover- 
lapping wave packets. Each one of these wave packets can be associated to one particular value of the measurement 
process. Second, he supplemented his measurement process by adding some obscure discussion about decoherence 
where only one of the several wave packets survives. After 1951, when Bohm presented his "hidden variables" theory, 
O [9], he kept the first part of his previous explanation of the measurement process and then provided a simple 
explanation about the final selection of the available wave packets: the selected wave packet is the one that contains 
the trajectoryp^ see Fig. 1.10(a). 

Hereafter, we provide a more detailed mathematical demonstration on how Bohmian mechanics explains the (von 
Neumann) measurement. Our explanation follows Refs. [501 [331 [Ml US] ■ In order to be able to define the measurement 
process in the Bohmian formalism, apart from the degrees of freedom xs, we need the degrees of freedom of the 
positions of the pointer belonging to a measuring apparatus Thus, we define a total wave function ^t{xs, xp,t) 
in a larger configuration space that includes the quantum region plus the measuring region, {xs, xp}. According to the 
Bohmian postulates, we do also select a trajectory {^^[t], a;p[t]} in this larger configuration space. The (Bohmian) 
time evolution of the total wave function and the total trajectory is all we need to explain any quantum process, 
including the (von Neumann) measurement process described earlier. 

The explicit consideration of the measuring apparatus is not trivial. In fact, there are some necessary conditions 
that the entire system has to satisfy to be able to say that a measuring apparatus is able to correctly determine the 
eigenvalues g. First, the pointer positions xp\t] of such an apparatus have to be restricted to a particular region, 
xp[t] G Sg^ , every time that the quantum system is in the eigenstate V'gi {xs)- Let us define ^g^ {xs, XA,t) as the total 
wave function that fits with the property that any experiment whose quantum system is described by ipg-^ {xs) implies 
that the pointer points in the particular region, xp [t] £ Sg-^ . A second necessary condition for a good measuring 
apparatus of the eigenvalue g is that Sgi n Sg2 — when we measure. We have defined the restricted region allowed 
by the pointer positions associated with a second eigenstate, i}}g^{xs), as xp[t] S Sg^. This implies that the states 
$gj^ (x5, x^, t) and (^g^{xs,XATt) do not overlap in the larger configuration space during the measurement. 

Since the eigenstates 'ipg{xs) form a complete basis, we can decomposed any function, and in particular each 
^g{xs,XA,t), into the following sum: 

^g{xs,XA,t) = ^fg'{xA,t) 1pg'{xs), (115) 

g' 

with fg{xA,t) ~ J ^g{xs,XA,t)ipg{xs)dxs- However, from our previous discussion about the properties of a "g ood" 
von Neumann measurement apparatus, ^g{xs,XA,t) cannot be a sum over different eigenstates 'ipg{xs) in Eq. (115) 



Very roughly speaking, Bohmian mechanics avoids the need to select randomly the final wave packet because the initial position of the 
trajectory is already selected randomly. 

In fact, we would have to include all others degrees of freedom of the "environment," xe, needed to be able to assume that the considered 
entire system is a closed one described by a pure (not mixed) state. 
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because then the measuring apparatus would provide the same pointer position xp[t\ G Sg for different eigenstates. 
Thus, during the time of the measurement, the only good decomposition for ^g{xsT XA,t) is: 

^g{xs,XA,t) = fg{xA, t) ^g{xs) (116) 

We emphasize that fg{xATt) is a normalized function because ^g{xs,XA,t) and tpg{xs) are also normalized functions 
in their respective configurations spaces. By its own construction, fg^{xA,t) H fg^(xA,t) = during the measuring 
time. Thus, even if ipgii^s) and '032 (^s) overlap, the states ^g-^{xs,XA,t) and ^g-^{xs,XAjt) do not overlap in the 
larger configuration space. See Fig. [ll]a. Thus, we can ensure that a general function in the quantum system, Eq. 



(114 1, can be rewritten in the whole, that is, quantum plus measuring, configuration space associated to a good 



measuring apparatus as: 

'^T{xs,XA,t) = '^Cg{t) fg{xA,t) i^gixs) (117) 
g 

In summary, during the time of measurement, the only total wave functions that can live in the entire quantum 



system that includes a good measuring apparatus of the eigenvalues g are the ones written as indicated in Eq. ( 117 1. 
An example of such wave functions is depicted in Fig. |11^ . It is important to notice that Eq. (117) implies no 



restriction on the wave function t/js^xs, t) but only on the total wave function $t(x5, XA,t)- If these restrictions are 
not respected, we can find other types of total wave functions in the configuration space {x5,xp}, but they would be 
incompatible with stating that we have an apparatus that is able to measure the eigenavalue g with certainty at time 
t. 

After this discussion, we can show quite trivially how the von Neumann measurement is exactly reproduced with 



Bohmian mechanics. As we have mentioned, apart from the wave function Eq. (117), we have to select an initial 



trajectory {2?s[0], a?^[0]}. Such a trajectory will evolve with the total wave function, and during the measurement. 



the particle trajectory {a?5[i], a;^[f]} will be situated in only one of the nonoverlapping wave packets of Eq. (117), 
for example fg^{xA,t)'ipg^{xs) as depicted in Fig. Illk Thus, the pointer positions will be situated in xp[t] £ Sg^^ 
and we will conclude with certainty that the eigenvalue of the quantum system is ga- In addition, the subsequent 
evolution of this trajectory can be computed from fg^(xA,t)^g^{xs) alone. In other words, we do not need the entire 



wave function Eq. (117) because the particle velocity can be computed from fg^{xA,t)^g^{xs)- The rest of circles of 
Fig. Ill are empty waves that do not overlap with fg^{xA,t)tpg^(xs) so that they have no effect on the velocity of 
the Bohmian particle. This is the simple explanation of how the complicated orthodox collapse is interpreted within 
Bohmian mechanics. This discussion is depicted in Fig. [TTJ 

We have mentioned before that the von Neumann measurement predicts that the probability of finding ga is 
Pg^ = \cg^{t)\'^. We can easily recover this result from our Bohmian picture. This probability is obtained after the 
repetition of the same experiment, with the same ^t{xs, XA,t)^ many times. The second (quantum equilibrium) 



postulate of Bohmian mechanics (see Sec. V.5 ) ensures that the modulus of the wave function will be reproduced by 
the trajectories when we repeat the same experiment many times. Then, the probability of finding the eigenvalue ga, 
that is, the eigenstate ipg^{xs), is just: 



p ^ I I\CgAt)fgaixA,t)^pgAxs)?dxsdXA ^ |Cg„(^)l^ .,x|2 /i.o^, 

nEJCg{t)fgixA,t)UxsWdxsdXA EJc,(OP ' ^ ^ ' ^ ' 



where we have explicitly used the condition that the each <I>gj (^5, xa, t) does not overlap with the others for the 
computation of the modulus of the total wave function, |$T(a?s, a?A, i)P, in the denominator. Additionally, we use 
/ \fpg{xs)\'^dxs = 1, / \ fg{xA,t)\'^dxA = 1, and |cg(*)P = 1- Thus, the same probabilistic predictions are obtained 
from Bohmian and orthodox quantum formulations when a von Neumann measurement is performed. 

Finally, we want to enlarge the explanation of the role played by the empty waves belonging to the wave packet of 
Eq. ( |117 ), which do not contain the particle. In principle, one can argue that such empty waves can evolve and, in 



later times, overlap with the original wave that contains the particle. If we are interested in doing subsequent (i.e., 
two times) quantum measurements, a good measuring apparatus has to avoid these "spurious" overlaps. This can be 
understood as an additional condition for qualifying our measuring apparatus as a good apparatus. In the examples 
of the subsequent subsections Sec. VI.2.1| and Sec. |VI.2.2 we will realize that the assumption of nonoverlapping 



between empty and full waves can be very reasonable in many practical implementations of the measuring apparatus. 
In fact, in a real system with a very large number of particles associated to the apparatus (or the environment), the 
probability that such overlap occurs in the configuration space of all the particles is almost null. See a discussion on 
page 79 in Ref. [82] . 

Finally, let us mention that we have only considered von Neumann (i.e., projective) measurements. Other types of 
measurements are also possible, which do not collapse into an eigenstate. Such measurements can also be explained 
within Bohmian mechanics with an extension of the ideas discussed here [531 [521 [551 ■ 
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VI. 2.1. Example: Bohmian measurement of the momentum 



Let us explain the momentum measurement of the stationary state of an electron initially in an energy eigenstate 
of a square well of size L. For example, ipn{x) — C sin(n7ra;/L) within the well and zero elsewhere, C being the 
normalization constant and n an integer (denoting the vibrational state). Since the wave function is real, we obtain 
p = dS{x)/dx = 0. Thus, the Bohmian particle is at rest, meaning that its velocity and kinetic energy are zero. 
However, we know that for a high enough value of n, the previous wave function ^{x) can be roughly approximated 
by a sum of two momentum eigenstates with eigenvalues ±nh/L. Therefore, when we do a (Copenhagen) momentum 
measurement, we will obtain an outcome ±nh/L and the previous wave function will collapse into one of the two 
momentum eigenstates. 

We have shown that the Bohmian measurement has to provide the same eigenvalues and eigenstates. How? The 
answer is that according to the Bohmian measurement, we have to specify in detail how the measurement process 
develops. In principle, for a formal Bohmian discussion of measurement, we would have to discuss the wave function 
and particles of the whole quantum plus measuring apparatus system. However, in this and the next didactic examples, 
we will avoid the wave function and trajectory related to the measuring apparatus. 

A measurement of the momentum, for example, can be undertaken by removing the walls that form the well and 
detecting the electron somewhere in a screen far from the initial walls. The time interval between removing the walls 
and detecting the particles will allow us a computation of the "electron velocity." The time-dependent process of 
removing the walls will imply that the stationary wave function will evolve into two time-dependent wave packets 
moving on opposite directions, which will become completely separated in space. The particle will end in one wave 
packet or the other with a momentum very close to ±n?i/L, the sign depending on which wave packet the initial 
position of the particle enters Finally, the electron will be detected far from the walls after the appropriate 
interval of time. For this particular experiment, we are sure that the empty wave will not affect the wave that carries 
the particle. As mentioned earlier, the relevant point is that we have been able to measure the same eigenstates and 
eigenvalues predicted by the orthodox interpretation without collapsing the wave function. 

In order to understand the Bohmian explanation of this experiment, it is important to discuss the roles played by 
the kinetic and quantum potentials. The initial Bohmian kinetic energy is zero, while the initial quantum potential 
is Q{x,t) = {nfi / / {2m) for the initial state il){x) = C sin{mTx/ L). On the contrary, the final energy of the system 
appears in the form of kinetic energy of electrons {nh / L)'^ / {2m) , while the final quantum potential from Eq. is 



negligible Q{x,t) 0. See problem 10 



VI. 2. 2. Example: Sequential Bohmian measurement of the transmitted and reflected particles 



Finally, we provide another didactic explanation of how Bohmian mechanics explains sequential measurement. The 
wave packets in Fig. 12 represent the solution of the (unitary) Schrodinger equation for a wave packet incident upon 
a tunneling barrier, at three different times. The initial wave packet (with norm equal to one) is divided into a 
transmitted plus a reflected wave packet. According to the Copenhagen interpretation, when the system is measured 
at time ti, a nonunitary evolution appears in the wave function and, randomly, the reflected wave packet disappears. 
Only the transmitted wave packet describes the electron at time ti. Then, when the system is measured again at t2, 
the electron is still represented by the transmitted wave packet. 

Alternatively, the same unitary and nonunitary evolution can be explained with Bohmian mechanics. The initial 
position of the Bohmian trajectory is selected randomly at the initial time Iq. Then, at times ti and ^2, the evolution 
of the trajectory is only determined by the transmitted wave packet. The reflected wave packet is an "empty wave" 
that has no effect on the evolution of the trajectory. Here, we have implicitly assumed that the eigenstates of the 
measuring apparatus are transmitted wave packets for positive currents, and reflected ones for negative currents. 

As expected, the probability of measuring, first, the particle as being transmitted at time ti and measuring, after, 
the same particle as being reflected at time ^2 is zero either with orthodox or Bohmian mechanics. 



VI. 3. The evaluation of a mean value in terms of Hermitian operators 



VI. 3.1. Why Hermitian operators in Bohmian mechanics? 



In the previous section, we have emphasized the ability of Bohmian mechanics to address quantum phenomena 
just from the Schrodinger equation and the equation of motion of the trajectories. The operators, which are an 
indispensable tool in the orthodox formulation of quantum mechanics to deal with the measurement process, become 
unnecessary in the Bohmian formulation. Once one realizes that operators are unnecessary, the conceptual difficulties 
associated with the wave function collapse of the orthodox formulation disappear. However, to be fair, the Bohmian 
measurement process explained in the previous subsections has quite limited practical utility. Including the measuring 
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FIG. 12. Schematic explanation of the abihty of Bohmian mechanics to discuss the unitary and nonunitary evolution of a wave 
packet incident upon a tunnehng barrier for three different times to < ti < t2. See also Color Insert. 

apparatus into the entire Hamiltonian is not an easy task, because it implies increasing the number of degrees of 
freedom that one has to simulate from {xs} to {xs,xa}- How many degrees of freedom are needed to describe a 
measuring apparatus? The answer is too many degrees of freedom. The exact Bohmian simulation of the measurement 
process is, most of the time, a quite difficult computational task. 

In these circumstances, the use of a Hcrmitian operator acting only on the wave function of the quantum system with 
the ability of providing the outcomes of the measurement process without the explicit simulation of the measuring 
apparatus is very welcomed. Let us emphasize that we are talking only at the computational level. In simple 
words, operators are not needed in Bohmian mechanics, but they are a very helpful mathematical trick in practical 
computational issues. These ideas are emphasized by Goldstein, Diirr, Teufel, and coworkers when they refer to the 
"naive realism about operators" [HSl [Ml EZ] . 



We now want to reproduce the mean value, obtained from orthodox Hcrmitian operators, within Bohmian mechan- 
ics. Someone who dislikes what we have explained up to here about the Bohmian measurement can just focus on the 
mathematical demonstration developed here on how Bohmian trajectories reproduce exactly the mean value of ortho- 
dox Hcrmitian operators. For simplicity, we will consider single-particle systems. The generalization to many-particle 
systems can be done straightforwardly. 

First, according to the orthodox measurement, we assume that a particular experiment is described by the quantum 
system wave function ip{x, t) alone. After a large number, M, of repetitions of the same experiment, we can elaborate 
a large list i = 1 . . . , / of the possible outcomes {a^} with their occurrence probabilities {Pi\- From this data we can 
compute the mean value as: 



VI. 3. 2. Mean value from the list of outcomes and their probabilities 



I 




(119) 



This is the standard (classical or quantum) definition of the mean value. 



VI. 3. 3. Mean value from the wave function and the operators 



According to the orthodox quantum theory, an observable is associated with a Hcrmitian operator A that describes 
the measurement process. Such an operator fulfills v4|iii) = ai\ui), {\ui)} being its eigenstates. Therefore, Pi — 
{^\ui) {ui\'\\)) so that we can write: 




(120) 



53 



where we use ai\ui) = A\ui) and we have identified the expression X]i=i 1'"*) ('"il ^ the identity because of the 
orthonormahty of the eigenstates. Therefore, we can compute the mean value of an ensemble of experiments from Eq. 
(1201 by only knowing the wave function and the operator (without knowing either the eigenstates or the outcome 
values and their probabilities). Therefore, we have been able to compute mean values within the orthodox quantum 
interpretation without directly mentioning the wave function "collapse." 



VI. 3. 4- Mean value from Bohmian mechanics in the position representation 



We can always write the Hermitian operator A and the mean value (A) in the position representation. Then, the 



(121) 



mean value of this operator over the wave function 'ip{x,t) is given by: 



A) = 

■4' 



ip*{x,t)A ( x,—ih— ) '4>{x,t)dx 



Alternatively, the same mean value can be computed from Bohmian mechanics by defining a spatial average of a 
"local" magnitude Ab{x) weighted by R^{x,t): 



A) = 



R^{x,t)AB{x)dx 



(122) 



In order to obtain the same value with Eqs. ( 121 ) and ( 122 ), we can easily identify the local mean value Ab{x) as: 

'r{x.t)A{x,~in-§^),l}{x,t) 



Ab{x) = Real 



rix,t)^ix,t) 



S(x.t) 

J 'ip{x,t)=R.{x,t)e R 



(123) 



We take only the real part, RealQ, because we know that the mean value must be real, and Eq. ( 123 ) without RealQ 
could take complex values. 



VI. 3. 5. Mean value from Bohmian trajectories 



For practical purposes, we will compute the mean value using Eq. (122 1 with a large j = 1,...,M number 

of Bohmian trajectories with different initial positions. We will select the initial position x^ [to] of the Bohmian 

trajectories according to the second postulate. Therefore, we can use Eq. (50l to write R^{x,t) in Eq. (122). Finally, 
we obtain: 



1 *^ 



(124) 



By construction, in the limit M oo, the value of Eq. (124 1 is identical to the value of Eq. (122 1. 

Now, we provide several examples of how some common mean values are computed from the orthodox quantum 
formalism and from Bohmian trajectories. First, we compute the mean value of the position: 



'(l]*(x, t)xip{x, t)dx 



(125) 



with xb{x) = X so that: 



R^{x, t)xdx 



(126) 



Identically, the mean value of the momentum: 



{P)^^ I V{x,t) [~ih-^]i}{x,t)dx 



(127) 
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with pb{x) = dS{x,t)/dx: 



iP)^- I R'M^-^^dx (128) 



For the classical potential, we have: 

/oo 
r{x,t)V{x,t)i^{x,t)dx (129) 
-OO 

with Vb{x) — V{x,t) so that: 

/oo 
R^{x,t)V{x,t)dx (130) 
-oo 

Now, we compute the mean value of the kinetic energy: 

(K)^, = J J*{x,t) (^-^^j ^{x.t)dx (131) 

It is important to notice that the local mean value of the kinetic energy takes into account the Bohmian kinetic energy 
plus the quantum potential. In particular, Kb{x) can be obtained from the expression: 

/ i?(x,t)e-^^(-'*)/''|^ (#)'i?(x,t)e*^(-'*)/''\ , , 



The real parQof Kb is: 

2m \ dx 

so that, finally, we obtain the Bohmian expression of the mean kinetic energy of the ensemble of trajectories 



^^^1/9^^ +Q(:^,t) (133) 



OO 



W„ = J- (^J +0(.,,)j<i. (134) 

See problem |10| for a detailed calculation and the explanation of the meaning of the two terms. 

Finally, we compute the mean value of the current density operator. First, let us remark that the probability 
density operator can be written as |a;)(a;| and its expected mean value is (iplx) {x\ip) — \ip{x,t)\^ or, in the Bohmian 
language, (tplx) {x\ip) — R^{x,t). The (Hermitian) current operator can be written as J = 1/ {2m){\x){x\p + p\x){x\). 
It can be easily demonstrated that: 



2 



1 

{J)^^J{x,t)^v{x,t)R\x,t)^\nTi —y^v{x^[t])5{x~x^[t]). (135) 



The average value of the current density depends on the position, and it is equal to the average Bohmian velocity 
multiplied by the square modulus of R{x,t). At a particular position "x," this current is just the sum of all particles 
that reside arround this position x = x^ [t] at time t. 



VI. 3. 6. On the meaning of local Bohmtan operators Ab(x) 

It is important to emphasize that the local Bohmian operators Ab{x) are not the eigenvalues of the operator A. 
In general, the eigenvalues are not positi on de pendent, while Ab{x) are. The expression Ab{x) is what is needed to 



compute the mean values of A with Eq. (122). This is its exact meaning. 



It can be demonstrated quite easily that the imaginary part of Eq. | |l32[ | is equal to the spatial derivative of the current that becomes 
zero when integrated over all space. We use J{x = itoo,t) = 0, which is always valid for wave functions that are normalized to unity, 
but it is not true for other types of wave functions such as plane waves. 
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There is a well-known example, mentioned in the previous Sec. VI.2.1[ that emphasizes the differences between 



eigenvalues and local Bohmian operators. We consider the particle between two walls, separated by a distance L, 
whose wave function is ipi^) — C* sin(n7ra;/i) within the walls and zero elsewhere. Since the wave function is real, we 
obtain p — dS{x)/dx = 0. Thus, the local Bohmian momentum is pb{x) = dS{x,t) / dx = 0, and the mean value in 



Eq. ( 128 1 is also zero. On the contrary, the wave function can be written as: 

^{x) w C (e""^/-^ - e''^^/-^) (136) 

We know that e^^-^^/L g^j^jj ^-m-nxjh eigenstates of the momentum operator with eigenvalues a\ = nhn/L and 
a2 — —nhn/L, respectively. Thus, the probabilities for each of these eigenvalues are identical, and the mean value of 



the momentum from Eq. (119) is again zero. 

In conclusion, in general, Ab{x) cannot be identified with a^. However, by construction, the mean values computed 
from ai and Ab{x) are identical. A similar discussion about the evaluation of the mean value of the quantum power 
can be found in Ref. lOTl. 



In this Sec. VI. 3 we have only discussed how the mean values (i.e., the first moment) obtained from orthodox 
quantum mechanics can be exactly reproduced within Bohmian mechanics. More complicated higher moments (e.g., 
the second moment called the variance, whose positive square root is the standard deviation) computed from orthodox 



quantum mechanics can be equivalently reproduced within Bohmian mechanics, as we anticipated in Sec. VI. 2. 2 
ESIIMIIHS. See also Ref. |92l. 



VII. CONCLUDING REMARKS 



In this chapter we have provided a fully comprehensive and didactic overview of the Bohmian formulation of 
quantum mechanics that can be useful for understanding the following chapters and also for any newcomers interested 
in using Bohmian mechanics in their daily research activity. 

To be fair, most of the research on quantum phenomena is developed within the orthodox machinery, and few 
significative contributions have been obtained with Bohmian mechanics. Many researchers ignore Bohmian mechanics, 
or they decide not to use it because they believe that someone demonstrated that there is something strange about 
Bohmian mechanics. In the introduction of this book, we illustrated this situation using Hans Christian Andersen's 
tale "The Emperor's New Clothes." The tale explains how a naked emperor strolling in procession on the streets 
believes he is well dressed, until a small child cries out, "The emperor is naked!" Children tales end happily with a 
moral that pretends to improve their behavior. Almost a century has passed, and we are still somewhere in the middle 
of the Bohmian tale. Bohm and Bell were the first to show to the scientific community that "Bohmian mechanics is 
a correct formulation of quantum phenomena." As in the tale and despite the whispering among the townsfolk, it 
seems that the idea that the orthodox formulation is the only licit tool to analyze quantum phenomena has gone on 
and the "orthodox emperor" continues to walk more proudly than ever. 

It is still not clear what ending awaits the Bohmian tale and what moral it will teach us. In our opinion, Bohmian 
mechanics will only be able to have a happy ending once it can be explicitly demonstrated that a particular meaningful 
quantum problem is better calculated (in terms of computational abilities) or formulated within Bohmian mechanics 
than within its orthodox counterpart. 

We hope that this happy ending is not very far and that this chapter (and the entire book) can contribute to 
shorten the time needed to arrive to it. If so, we might learn a moral similar to the one deduced from Andersen's 
tale: "It is useful to look for answers outside of the main stream." Sometimes the scientific revolutions are hidden there. 
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VIII. PROBLEMS 



Problem 1 Check that Eq. [22) is fulfilled for a free particle with V{x^t) = 0. We define the initial time as to = 



and the initial position as Xq ~ 0- 



Solution 1 For a free particle, the Lagrangian is just the kinetic energy with a constant velocity, v[t] — vq. We have 
Sp — K(vo)tf so that Sp = mxf'^/{2tf) and Sp-st — rnxf'^/{2{tf +6t)). Notice that modifying the final time, without 
modifying the final position, means changing the initial velocity. If we use the Taylor expansion 1/(1 + a;) ~ 1 — a;, 
then we find Sp-st = rnxf"^ (1 — St/tf)/{2tf). Finally, Sp-st ^ Sp = —mxf'^ 6t/(2tf'^). This is the kinetic energy 
with a negative sign that, in our simple system, is also the Hamiltonian with a negative sign. Therefore, Eq. \2^ is 
exactly recovered. 



Problem 2 Following the example of problem^ test Eq. (26). 



Solution 2 The new value will he SpSx — + (5a;)^/(2i/). If we neglect we obtain Sp-Sx = '7ia:/^/(2iy) + 



mxfSx/tf . So that Sp-Sx — Sp = m vS x. Finally, Eq. {26) is recovered. 



Problem 3 Solve Eq. (21) for the free particle case developed in problems^ and[£| to find S{xq, t^; x, t). 



Solution 3 Equation Eq. \21/^ can he rewritten as dS{x,t)/dt = —H{x,t) = —E{vq). The system is conservative, 
and the value o f th e Hamiltonian at any time is equal to the initial kinetic energy E{vq) = mi'o/2. Therefore, we can 



transform Eq. (21) into {dS{x,t) /dx)"^ / {2m) = E{vq) so that dS/dx = ±y^2mE{vQ) — mvQ. If we assume that the 
initial velocity (energy) is independent of the initial (final) position, then we easily integrate the partial derivatives of 
S to obtain finally, S{x,t) = zLmv^x — E{vo)t. 

As expected, this action function corresponds to trajectories with constant velocity ±vo, that is, x[t] — xq ± Vq t. 
The important point of our calculation here is that we have found all trajectories with an arbitrary initial position xo 
but a fixed initial velocity dS/dx — imwo- Let us emphasize that, although the classical function S is multivalued, if 
the initial velocity was -\-vo, then we can use this knowledge and keep always the positive velocity for such trajectories. 



Problem 4 Show that the two real Eqs. ({21^) and (\2^) are equivalent to the following complex classical wave Eq. 
(29 



.^dii)ci{x,t) 



d^^ci{x,t) 

2m dx^ 



+ V{x,t)'4jci{x,t) + 



d^\Mx,t)\/dx^ 



2m 



\'ipcl{x,t)\ 



1pclix,t) 



(137) 



Solution 4 14^6 start the formal demonstration by introducing the polar form of the wave function, ipci { x, t) = 
R{x,t)exp{iS{x,t)/h), into Eq. (29). We have explicitly included the time in the scalar potential of Eq. ^2S^ in 



order to compare it with the Schrddinger equation. Then, we evaluate the temporal and spatial derivatives: 



d_ 
dt 



(Re'^'f^) 




i_dRdS_ 
h dx dx 



^R^ 
h dx"^ 




js/h 



(138) 



(139) 



Now, if we introduce Eqs. ((138)) and ((139)) into Eq. (29) and eliminate the factor e'"^/^, which is present in all 
terms, we obtain: 



ih 



dR 
~dt 



R 



as; _ _^d^ „ A^^^ 

dt 2m dx"^ m dx dx 



^±R^ 
2m dx"^ 



Zm \ ox 



V R 



r d^R 

2m dx^ 



(140) 



The real part of Eq. (I4O) is 



R 



as 



2m V ox 



V R 



(141) 



which is exactly the Hamilton-Jacobi Eq. [21). Note that the last term in the r.h.s. of Eq. (131) is eliminated by the 
first term in the r.h.s. of Eq. (139). 
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On the other hand, the imaginary part of equation Eq. { I4O) gives: 

1 



OR 
~dt 



1 dRdS 



R 



m dx dx 2m dx 



(142) 



Now, if we multiply both terms of Eq. (14-2) by 2R, we directly recover the conservation law of Eq. (28). Notice how 
the constant h does not appear either in Eq. (I4I) or in Eq. (142). 



Problem 5 Show that ^ci{x,t) = exp{i{±px — Et) /h) is a solution of the classical wave equation of Eq. (29) for the 
free-particle system with V{x,t) — 0. 



Solution 5 We compute the first two terms of the classical wave Eq. (29). 

d 



dt\ J h 



2m dx^ 



2m 



(143) 
(144) 



Since the other two terms are zero, we obtain E = P^ / {2 m). Thus, ipci{x,t) = exp{i(zLpx — Et)/h) is a solution of 
the classical wave Eq. (29). 

Surprisingly, a plane wave, apart from a free-particle quantum solution of the Schrodinger equation, is also a licit 
solution of a wave equation for classical mechanics. The fact that both classical and quantum solutions are identical is 
due to the fact that for the particular solution discussed here, R — 1, the last term of Eq. (29^ , that is, the quantum 
potential, is zero. Whether we want to model a free particle as a plane wave or as an ensemble of trajectories is just 
a matter of convenience. 



Problem 6 Show the following identity for the computation of the Bohmian velocity. 

1 dS{x,t) J{x,t) 



v{x, t) 



m dx \ili{x,t)\'^ 



(145) 



where J(x,t) is defined by Eq. (34)- 



Solution 6 We start by performing the time derivative of S{x,t): 



dS{x,t) 
dx 

dS{x,t) 
dx 

dS{x,t) 



ft— — tan 
ax 



^ijXyt) 

A{x,t) 



= h- 



dx 



d 4'i(x,t) 
dx ^r(x,t) 

( iHix^Y 
^ \ i^Ax,t)) 

^ ^.(x,t)g^-V,.(x,t)g^ 
i;^{x,t)+^f{x,t) 



Finally, using the following identity: 



ipix,t) 



dip*{x,t) 



rix,t) 



dilj{x,t) 



- A{x,t) 



dip^jx^t) 
dx 



dx dx 
we obtain that both definitions of the velocity of quantum trajectories are identical: 

1 dS{x,t) _ J{x,t) 



ipiix,t) 



dlpr{x, t) 

dx 



v{x, t) = 



m dx 



\i^{xM 



(146) 



where J{x,t) is defined by Eq. (34). 



Problem 7 Show that the Bohmian velocity definition 



v{x, t) 



J{x,t) 

mx,tw 
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given in Eq. (35) fulfills the continuity equation, Eq. l3!^\j : 

dJ{x,t) 



dt 



dx 



J{x,t) being the current defined in Eq. [34): 



As a conclusion, we learn that a proper (i.e., quantum equilibrium) ensemble of Bohmian trajectories will reproduce 
R{x,t)^ = \ip{x,t)\'^ at any time. 

Solution 7 We rewrite the continuity equation as: 



dt 



R^{x',t)dx' + J{x + Aa;,i) - J{x,t) = 



By using the Bohmian definition of the velocity, we rewrite Eq. (l^l) as: 

i*x+Ax 



d_ 
di 



R\x',t)dx' = + Aa;,t) v{x + l\x,t) + R\x,t) v{x,t) 



(147) 



(148) 



Now, we make use of the ensemble of Bohmian trajectories by using Eq. (50) to write i?^(a: + Aa;, i) and R^(x,t) as 
a sum of delta functions at t = to- Finally, we obtain: 



_9 
dt 



/ R^(x' ,t)dx' ^ + lim -^^^ Six - xAtOvAt) - lim — V (5(a; + Ax - a;Ji])w,(t) 

M-)-oo M ^-^ M^oo M ^-^ 



We can rewrite the 8 as time derivatives of the Heaviside (or unit step) functions Q: 

rx+Ax 1 M 



d_ 
di 



/ R^{x',t)dx' ^ + lim —y" — (e{x,[t]-x)e{x + Ax-x,[t])) 

J ^ A/— i-oo M ^ — ' Ot 



(149) 



(150) 



Thus, counting the number of Bohmian i particles inside the region [x, x + Ax] does exactly reproduce the square 
modulus of the wave function and its time- dependent variations. In conclusion, if we know that the ensemble of 
Bohmian trajectories reproduces R{x,t()) at the particular time to, then it will reproduce R{x,t) at any other time: 



1 

R\x,t)= S^^i^-^dt]) 



(151) 



This is precisely the reason why it is claimed that Bohmian mechanics exactly reproduces the position measurement 
of orthodox quantum mechanics. Since an ensemble of Bohmian trajectories reproduces 'ip{x,t) at any time, the 
ensemble of trajectories will also reproduce the mean value of other observables when they are written in the position 
representation. 



Problem 8 Demonstrate using Bohmian trajectories that the (average) dwell time (t) that a particle described by the 
wave function ip{f,t) spends in a particular region r G can be written in an orthodox language as: 

(r) ^ r dt [ dv\^{f,t)\' (152) 



Solution 8 For a particular Bohmian trajectory rk[t\, the time spent in a particular region for this k-particle can 
be unambiguously written as: 

Tfe = / dt^ dte{t - 4)0(ifc - 1) (153) 
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where t^, and are the times when the electron enters and leaves the volume Q, respectively. The function Q{t) is a 
Heaviside function (or unit step function). We can relate the Heaviside function to the delta function S as: 



e{t^ti)e{ti-t)= / dvs{f-r,[t]) 



(154) 



In order to obtain the average value of the dwell time, (r), we have to make a sum over the ensemble k = 1, . . . , M 
of Bohmian trajectories: 



^ M ^ M poo n 

1™ TrT^'^fc^ 1™ 77^1 / '^M dv6{r-fk[t]) 
M->oo M ^ M^oo M J ^ Jo 



(155) 



k=l k=l' 

Finally, by noting that the ensemble of Bohmian trajectories reproduces the modulus of the wave function: 

M 



|V'(r,i)|'= lim i^E'5(f~rfc[t 

A/— ^oo ivl ^ — ^ 



(156) 



k = l 



obtain the expected results, Eq. (152). See more information on tunneling times in f88l 1 90\/ . 



Problem 9 Let us consider a single Bohmian trajectory Xk [t] accounting for the quantum dynamics of a ID particle 
in a potential V{x). Show that the left probability density: 



Pl 



Xklt] 



\'ilj{x,t)\ dx 



(157) 



is constant, which, in turn, means that ID Bohmian trajectories do not cross each other. 
Solution 9 We are interested in .showing that: 

dPL d /-^"W 



We know that: 



d /■^'=[*1 



dt dt 



dx\i;{x,t)\^ ^\i>{x,t)\^ 



dx\i}{x,t)\^ = 



dxk[t] , T'^W, d\4'{x,t)\ 
' ' dx- 



dt 



+ 



dt 



= 



By the continuity equation, Eq. [33), with the current density in Eq. (34-), we obtain: 

dxk[t] dj{x,t) 



dx ^X^'"' = 
ax 



We finally obtain -^Pl — 0, that is, Pl is a constant when the Bohmian trajectory moves according to: 

dXk[t] J{Xk[t],t) 



dt \lb{xk[t],t)(^ 



(158) 



(159) 



(160) 



(161) 



This is exactly the Bohmian velocity of Eq. (4-2). 



Problem 10 Show that the expected value of the kinetic energy, where the kinetic operator is given by T = 
-{h^/2m)V'^, yields: 



T) = (Q) + 



2m 



(162) 



where Q is the quantum potential. Equation (162) indicates that the quantum potential Q could be also interpreted as 
an additional quantum kinetic energy. 
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Solution 10 By writing the wave function in polar form, ip = Re^^/^ , we calculate the expression of the expected 
value of the kinetic energy by averaging over all the possible positions of the particle, that is, in a (infinite) volume V: 



T 



'4}*T%l]dV 







j r 




2m 



2m 



dV 
dV 



2m 



V • ^Re:^'^- + -R^Se'^'^ 



R 



2m 



2i. 



-VR-VS 



h 



:RV'S- 



1 



R' - 



2m R 



dV 



2m 



V • ( R^VS] dV 



dV 



R{VS)' 



R' 



dV 



2m 



{vsy dv 



(163) 
(164) 
(165) 
(166) 
(167) 



The first and third terms in the r.h.s. of Eq. 
theorem on the integral in the second term, leads us to: 



V • ( R'^VS] dV 



167) are, respectively, {Q) and V^- The application of the divergence 



ndS — 0, 



(168) 



S being the surface boundary enclosing the volume of integration and ft the outward pointing unit vector perpendicular 
to it. Equation (168) is zero since the population must be conserved, and there can be no flux of probability outside 
the (infinite) volume V . 

Thus, we see that the expected value of the quantum kinetic energy operator corresponds to the sum of the expected 

values of the classical kinetic energy and the quantum potential. This discussion is related to what we have 
explained in Sec. \ VI.2.1\ 

Problem 11 Show that for a stationary state, ijjix,!) = ip{x)e~^^^^^ , that is, for an eigenstate of the Hamiltonian, 
it follows that: 



P'^ix) = 2m[E- V{x) - Q{x)] 



(169) 



where P{x) = dS/dx is the Bohmian momentum and E — —dS/dt the energy eigenvalue of the corresponding energy 
eigenstate. Note that Eq. 1 16£lj) su ggests that the semiclassical (i.e., WKB) approximation for the particle momentum 
given by p(x) — ±y^2m[E ~ V{x)] is accurate at those positions where the quantum potential is relatively small. 



Solution 11 For an eigenstate of the Hamiltonian, it follows that: 

dR 



dt 



= 



dS 

Ik 



= -E 



Then, by using the quantum Hamilton- J acobi Eq. \3(^ : 

dS 



-V 



with 



one obtains 



(dS\ 
2m \dx) 



1 d'^R 
2m R dx"^ 



Q 



P^{x) = 2m [E - V{x) - Q{x)] 



(170) 
(171) 

(172) 

(173) 
(174) 



where we have used the Bohmian definition of the momentum P{x) = dS/dx. 
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Problem 12 Recently, there have been several investigations where a quantum trajectory extension of Bohmian me- 
chanics has been performed taking the quantum action S to he complex U06\. \108\j . Assuming the following expression 
for the wave function: 



■)p{x,t) 



,- S(x,t) 



(175) 



derive the complex quantum Hamilton- J acobi equation associated with it. 

Then, assuming that v{x,t) = (S/ S {x , t)) / m and taking the time derivate of this Hamilton- Jacobi equation, obtain 
the Newtonian-like equation of motion for the velocity: 



dv [xq , t] 
dt 



m 



ih 
2m 



(176) 



where the quantum features, such as nonlocality, are manifested in the second term of the r.h.s. of Eq. [176). 



Solution 12 Casting il;{x,t) = e* i 



where 



into the time- dependent Schrddinger equation yields: 

ih 



V 



2m 



2m 



(177) 



(178) 



is the "new" quantum potential. The first step to obtain the "new" Newtonian-like equation of motion is to take the 
spatial derivative of the complex quantum Hamilton- Jacobi equation: 



dv - \ ih-±^ 



Then, by identifying the total derivative of the velocity: 

dv dv 
dt ~ dt^ 

one finally arrives at 

dv[xo,t] VV 



vVv 



ih 



-V^v 



(179) 



(180) 



(181) 



dt m 2m 

The terms at the r.h.s. can be identified as the classical and the quantum force acting on the particle, respectively. 

Problem 13 The Gross-Pitaevskii equation governs the dynamics of a Bose-Einstein condensate (BEC) in the so- 
called mean field approximation. In the zero temperature limit, it reads: 



^^ojppt) 

at 



_ v2 + y(f,i) + 5|^(f,t)p 

2m 



(182) 



where g is the nonlinear interaction between the atoms of the BEC, V{f, t) is the trapping potential, and the BEC 
wave function is normalized to the number of atoms as ^ d^r\ip{f,t)\ — N. Note that for g = 0, one formally 
recovers the Schrddinger equation. Apply Madelung 's formulation to the BEC by writing down its wa ve fu nction as 
'ip(f,t) = yjn(f, t)e*^'- and obtain the continuity and Hamilton-Jacobi equations associated with Eq. (182). 



Solution 13 First of all, notice that following the usual nomenclature of the literature on this Eq. (182), tp{r,t) = 
y/n{r,t)e''^^^^^\ here we use different functions for the role of R{r, t) and S{r, t) that are related, hut not identical, to 
= n{r,t) and 9{r,t), respectively. Following the same steps as in Sec. IV. 3. 4 it is straightforward to obtain: 

dv 



dn 
di 



V + ng 
V(ntT) 



1 



2m y/n 



(183) 
(184) 



where we have defined v{f, t) = ^^V6{r, t). Equation ( 183) is the classical Euler equation for an inviscid compressible 
fluid with an additional pressure term whose origin is purely quantum (the last term in this expression) that is very 
often called quantum pressure. These equations will be used later in chapter 3. 
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Problem 14 Use the many-particle Schrddinger Eq. (52) to deduce a many-particle local continuity equation. 



Solution 14 In order to find a local continuity equation, let us work with ijj{x,t) and its complex conjugate ip*{x,t). 
In particular, we can rewrite Eq. (30) as: 



. dip{x, t) 



dt 



^fc=i -^k 



r{x,t)in^^^^^ = r{s,t) [Y,-^^^Ks^t) 1 + r{x,t)v{x,t)i;ix,t) 



and the complex conjugate as: 



Tpix, t)ih 



d^*ix, t) 
dt 



^r^{x,t)[Y,- 



\k=l 



2^^^2^*(2^'*) 1 + ^{x,t)V{x,t)r{x,t) 



When the first equation is subtracted from the second, we obtain 

d\i,{x,t)\ 



dt 



Y,Jk{x,t) = 



k=l 



where we have used: 



d d 

^(^'*)^ — - V{x,t)- — ■0(^,i, 
oxk oxk 



and defined: 



j/^,. f ,.9}P*(At) ^^ dip{x,t) 



(185) 



(186) 



(187) 



(188) 



(189) 



as the k-th component of the current density. 



A.l. APPENDIX: NUMERICAL ALGORITHMS FOR THE COMPUTATION OF BOHMIAN 

MECHANICS 



In this section, we will describe the basic numerical algorithms available in the literature for the computation of 
Bohmian mechanics. An excellent source of information on this issue can be found in Ref. |16j . There are many possible 
algorithms that can be classified following different schemes. One can distinguish between analytical or synthetic 
algorithms]^ An analytical algorithm, first, solves the Schrodinger equation and, then, evaluates the Bohmian velocity 
directly from the numerical wave function. In such algorithms, the computational difficulties are directly related to 
the numerical solution of the Schrodinger equation. Alternatively, one can choose computational algorithms where 
the Schrodinger equation is substituted by the quantum Hamilton- J acobi equation (synthetic approach). Then, the 
difficulties of the algorithms are intrinsic to the numerical resolution of the (nonlinear) Hamilton-Jacobi scheme. In 
addition, the quantum Hamilton-Jacobi equation can be solved with or without trajectories. Another classification 
divides algorithms between time-dependent and time-independent ones. From a computational point of view, the 
former is basically related to an eigenvalue problem, while the latter is related to an initi al va lue problem. Finally, a 



classification between ID spatial grids and 2D (or even 3D) is also possible In Table (A.l), we have summarized 
the characteristics of the numerical algorithms that we have discussed. 

Much of our attention in this section will be devoted to time-dependent scenarios. In such scenarios, either the 
Schrodinger equation or the quantum Hamilton-Jacobi equation can be written as a first-order differential equation 
in time of the form: 



The names "analytical" and "synthetic" are related to philosophical concepts. The original name of analytic propositions referred to 
those propositions that were true simply by virtue of their meaning, while synthetic propositions were not. 

Certainly, the direct solution of more than three degrees of freedom is very computationally demanding. To the best of our knowledge, 
the only Bohmian algorithm that takes into account correlations for more than three degrees of freedom is the one mentioned in Sec. 
1.3.6. This many-particle Bohmian algorithm does not fit with the previous classification because it is a mixed analytical and synthetic 
algorithm. 
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TABLE A.l. ID and 2D Bohmian computational algorithms (TD — time dependent; TI = time independent; SE = Schrodinger 
equation; QHJE — quantum Hamilton-Jacobi equation; BT — Bohmian trajectries). 







Analytical 


Synthetic 






no-trajectory trajectory 


no-trajectory trajectory 


Spatial ID 


Time-dependent 
Time-independent 


— TDSEiB + BT 

— TISEiD + BT 


TDQHJEiD TDQHJEio + BT 
TIQHJEiD TIQHJEiD + BT 


Spatial 2D 


Time-dependent 
Time-independent 


— TDSE213 + BT 

— TISE2D + BT 


TDQHJE2D TDQHJE2D + BT 
TIQHJE2D TIQHJE2D + BT 



TABLE A. 2. Bohmian time-dependent computational algorithms. 







Type of solution 


SE-> linear f{W[t),t) 


QHJE-> nonlinear f{W{t),t) 


Implicit 


dW{t) 
dt 


Unconditional stable 


Matrix inversion 


Newton-like algorithm 


Explicit 


dW(t) 
dt 


Conditional stable 


Explicit solution 


Explicit solution 






Unstable 


No solution 


No solution 



where W{t) can take the role of ip{x,t), R{x,t) and S{x,t)). f{W{t),t) is a generic function that depends on W{t) 
and time t. We define a mesh for the time variable tj = jAt for j — 1, . . . , Af and assume that the initial value of 
VF(0) is known. Then, one can use an exphcit finite difference of the time derivative: 



dW{t) 



dt 



W{t,+,)-W{t,) 
At 



(A.2) 



giving a simple and explicit numerical solution of the time-dependent Eq. ( A.l ) as: 



Wit,+,)^Wit,) + f{Wit,),t,)At 



(A.3) 



where the solution W{tj^i) at tj+i can be found directly from the knowledge of W{tj) and f{W(tj),tj) at the 
previous time tj. These kinds of solutions are named explicit solutions (also known as forward Euler method) and 
are, in general, unstable because the error during successive time integrations grows exponentially. Only for a very 
dense time mesh, that is, a very small At, the solution becomes conditionally stable. 
On the contrary, we can look for an implicit approximation of the derivative: 



dW{t) 



dt 



At 



(A.4) 



that provides an implicit solution (also known as backward Euler method) of the time-dependent Eq. (A.l) as: 



Wit,)^W{t,^,) + f{W{t,),t,)At 



(A.5) 



The practical numerical solution of Eq. (A.5) depends on the type of equation. For instance, for the Schrodinger 
equation, the function f{W{tj),tj) is linear with respect to W{tj), and then, the relationships between W{tj) and 
W{tj-i) can be written as a set of coupled linear equations, one for each grid point. Then, the entire system of linear 
equations can be written as the product of a (known) matrix and an (unknown) vector. The (unknown) vector is 
formed by the unknown value of W{tj) at each grid point. The solution follows from inverting the matrix. On the 
contrary, for the quantum Hamilton-Jacobi equation, where f{W{tj),tj) is nonlinear, one has to look for the numerical 
solution of W{tj) — W{tj-i) — f{W{tj),tj)At = through a Newton- like procedure. Apart from the important increase 
of the computational burden, both types of implicit solutions provide a nongrowing error independently of the value 
of At. For this reason, they are called unconditionally stable algorithms. See a summary of the classification in Table 



(A.2) 



Hereafter, we will explain in detail some of these numerical algorithms, discussing their main advantages and 
drawbacks. We will mention some very simple algorithms that can be useful for those newcomers who are just 
interested in computing Bohmian trajectories. 
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A. 1.1. Analytical computation of Bohmian trajectories 

According to the first postulate of Bohmian mechanics, the computation of trajectories is quite simple once the 
wave function is known. As we mentioned in the introduction, here, the utihty of using Bohmian trajectories is that 
allows an intuitive and different explanation of many quantum phenomena. In particular, in this section, we will 
discuss three simple numerical methods for solving the ID Schrodingcr equation: 



A. 1.1.1. Time- 



ndent Schrddinger equation for a ID space (TDSEid-BT) with an explicit method 



In this subsection, we will explain a very simple explicit algorithm to compute the time-dependent Schrodinger 
equation. Apart from a mesh in the time variable, an additional mesh is imposed on the ID spatial degree of freedom 



with a constant step, Ax, that is, — k Ax for k = 1, . . . , N . Then, we define ip(x,t)\x= 



:Xk;t = tj 



^pj{xk), and the 



temporal and spatial derivatives present in the time-dependent Schrodinger equation can be approximated by: 



d^/j {x, t) 



dt 

d^^l^{x,t) 



dx^ 



2At 

tpj{xk+i) - 2tljj{xk) + •0j(xfc_i) 
Ax^ 



Inserting Eqs. ((A.7|) and ((A. 8)) into Eq. (A. 6), we obtain the following simple recursive expression: 



Ax m ii 



(A.7) 
(A.8) 

(A.9) 



Once we know the wave function at the particular times tj a nd tj-i for all spatial positions in the mesh, we can 
compute the wave function at next time ij+i, using Eq. (A.9 1. This algorithm solves an initial- value problem. The 
recursive application of Eq. (A.9) provides the entire time evolution of the wave packet. Equation (A.9) is not valid 



for the first, Xi, and last, cc^v, points. To avoid discussions on the boundary conditions, we can use a very large spatial 
simulating box so that the entire wave packet is contained in it at any time. Then the wave function at the borders 
is negligible. 

As indicated in Tables A.l and A.2[ this explicit solution can be unstable, and its error grows in each recursive 
application of Eq. (A.9 1. To provide a (conditional) stable solution we have to deal with small values of At and 
Aa;. For example, to study electron transport in nanoscale structures as depicted in Figs. |A.l| and |A.2[ this recursive 
procedure provides accurate results (the norm of the wave packets is conserved with high precision) when Ax is on 
the order of 1 or 2 A and the temporal step. At, is around lO^^^s. 

In order to define the initial value of the wave function, we can consider that the wave packet at times t = {ti,ti+At} 
evolves in a flat potential region contained in a much larger simulation box. Then, for example, we can define the 
initial wave function as a time-dependent Gaussian wave packet |42| : 



^{x,t) = 



2a' 



1/4 



1/4 



exp 



xa 



{t-h)Y 



2ih{t-ti) 



(A.IO) 



where a is the spatial dispersion of the wave packet, m the particle effective mass, xq the central position of the wave 

packet at the initial time ti, kc — \J '^j^ the central wave vector in the x direction related with the central energy 

E, and (f) = —9 — Tikc'^t/ {2m) with tan(26') ~ 2ht/ {ma?) (see [H]). In particular, at the initial time t = fi, we obtain 
the simplified expression: 



il;{x,t = ti) = 



1/4 



e^ik^^-^o)) exp 



{x - Xo)^ 



(A.ll) 



We define the wave packet spatial dispersion a.^ — a/^/2 and the wave packet width in the reciprocal space as 
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Distance (nm) 

FIG. A.l. Representative Bohmian trajectories associated with double-packet scattering in a double-barrier potential. The 
position of the barriers is indicated by vertical lines. Reprinted with permission from [9U] . Copyright 1996 American Physical 
Society. 



To finish this practical computation of Bohmian trajectories from the direct solution of the time-dependent wave 
function, we can discuss some features of Bohmian trajectories related to the fact that they do not cross in the 
configuration space. We consider a quite "exotic" initial wave packet. We use a sum of two Gaussian wave packets 
defined by the expression Eq. (A. 10) with different central positions and central (momenta) wave vectors [50] • The 
wave packet is certainly quite exotic because it describes just one particle. In Figs. |A.l| and |A.2| we see that Bohmian 
trajectories can be reflected for two different reasons: first, because of their interaction with the classical potential 
(the particles collide with the barrier) and, second, because of the collision with other trajectories traveling in the 
opposite direction (Bohmian trajectories do not cross). The second process is responsible for the reflection of those 
particles of the first packet, which never reach the barrier, and for the reflection of the entire second packet. These 
collisions between Bohmian particles are related to the quantum potential in regions where the classical potential is 
zero, but for them to occur, there should be particles coming from right to left. In this regard, if the initial wave 
packet is prepared as a superposition of eigenstates incident from left to right (as is always assumed in scattering 
thought experiments), and the classical potential is zero on the left-hand side (l.h.s.) of the barrier, then finding 
particles coming from the r.h.s. in this region will be at least very uncommon. 



A. 1.1. 2. Time-independent Schrddmger equation for a ID space (TISEm) with an implicit (matrix inversion) method 



Next, we consider another very simple computational algorithm similar to the one discussed in the previous section 
(ID, analytic and explicit) but for time-independent systems. First, let us discuss under which circumstances we can 
expect that a time-independent Schrodinger equation can provide a reasonable physical description. Let us assume 
that the wave function ip(x,t) can be written as a product of a spatial function ip{x) and a temporal function C(^)i 
that is, ip{x,t) = ip{x)C,{t). Then, the Schrodinger equation, Eq. (A.6), with a time-independent potential V{x) can 
be written as: 



m 



2m dx^ 



V{x)ip{x) 



ip{x) 



(A.12) 
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FIG. A. 2. Picture motion of the same double-packet wave function considered in Fig. |A.1| calculated by numerical integration 
of the time-dependent Schrodinger equation. Five representative "snapshots" obtained at different times are shown with the 
vertical scale arbitrarily changed in each case for clarity. Although the norm of the wave function is always unity, it does not 
seem so because of the changes in the scale. The marks are visual aids that indicate the position of some related Bohmian 
trajectories shown in Fig. |A.1[ The double-barrier position is indicated by the vertical dashed lines, and the arrows indicate 
the sense of motion of the two packets. Reprinted with permission from |93. Copyright 1996 American Physical Society. 



Since the l.h.s. of Eq. (A. 12) only depends on t and the r.h.s. of Eq. (A. 12) only depends on x, both sides are equal 



to a constant that we refer as E. Therefore, we obtain two separate equations: 



2m dx^ 



+ V{x)(p{x) = Elp{x) 



(A.13) 
(A.14) 



The first of these equations has a trivial analytical solution ({t) — e'^ n . However, the time-independent Schrodinger 
equation needs a numerical method to be solved for arbitrary shapes of the potential V{x). 

The second-order derivative can be approximated as: 



dx^ 



ip{xk+i) - 2ip{xk) + ^pjxk-i) 
Ax^ 



(A.15) 



Then, Eq. (A.14 1 can be written in a matrix representation as: 



2m 



V' + iV{xk)~E)I\^^0 



(A.16) 



where (p = {(p{x2), ■ ■ ■ ,(p{xN^i))'^ is a — 2 vector, / is the {N — 2) x {N — 2) identity matrix and is the 
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(A'' — 2) X (iV — 2) matrix constructed from Eq. (A. 15 1 as: 



\ 



2 10 

-2 10 



1-21 

1 -2 y 



(A.17) 



The solution of can be found by inverting the matrix equation, Eq. (A. 16). The (nonzero) solution of Eq. (A. 16) 



requires that the determinant of the matrix be equal to zero. The determinant requirement is called the characteristic 
equation (less often, the secular equation). In our case, the determinant can only be zero for some particular values 
of the energy E (i.e., quantized energies). Thus, this algorithm finds the eigenvalues E and eigenvectors ip, typical of 
a boundary value problem]^ 

In the previous vector definition, only N — 2 unknowns are considered. We have assumed that f{xi) = and 
f{xN) — 0. This implies that is a, real vector. Then, if we select an eigenstate, we realize that the current density 
is zero, J{x,t) — 0, because ip{xk) = ip{xk)*- AH Bohmian velocities are zero, dS/dx = 0, meaning that Bohmian 
trajectories do not "move." This result is obviously consistent with the fact that the amplitude R{x, t) of the global 
wave function ip{x,t) = ip(x)e^^^*^^ is time independent. One is tempted to affirm that there is no kinetic energy 
for these Bohmian trajectories. However, as we see in Eq. (134), the evaluation of the mean value of the kinetic 
energy takes into account not only the Bohmian velocity but also the quantum potential. We can explain this point, 
differently, from the quantum Hamilton-Jacobi equation, Eq. (39). In our particular time-independent case, we have 
an energy conservation law: 



V{x) 



d^R{x)/dx^ 
2m R{x) 



E 



where we have used dS/dt = —E in Eq. (39). Then, we can obtain the following (time-independent Schrodinger-like) 
equation for the real wave function R{x): 



d^Rjx) 
2m dx^ 



+ V{x)R{x) = E R{x) 



(A.18) 



Comparing this expression with the previous Eq. ( A.14 ), we realize that in this particular problem (when all Bohmian 
trajectories have zero velocity), the quantum potential takes the role of "kinetic" energy. 



A. 1.1. 3. Time-independent Schrddinger equation for a ID space (TISEid) with an explicit method 

We can also look for a very simple explicit algorithm to find the solution (f of the time-independent Schrodinger 
equation. Equation (A.14) is a homogeneous second-order differential equation that can be written as: 



-^V'(x) = f{x)ip{x) 



where f{x) is the constant of proportionality between the solution and the second derivative at each point: 

ym 

f{x) = -j{V{x)-E) 



From the Taylor expansion for (p{xk+i) around x^, we get: 



(p{xk+i) = (p{xk) + Ax 



dip{x) 



dx 



Ax-^ d^(p{x) 



21 dx^ 



Ax^ d^(p{x) 



3! dx^ 



Ax* d-^ifix) 



4! dx^ 



Ax^ d^(p{x) 
~J\ dx^ 



0{h') 



(A.19) 



(A.20) 



(A.21) 



The practical implementation of this algorithm can be done, for example, quite easily with the software MATLAB. The instruction 
[V,D] = eig(A) produces matrices of eigenvalues (D) and eigenvectors (V) of matrix A, so that A * V = V * D. Equation | |A.16[ | can 
be rewritten a.s A*V = V*D. 
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Identically, for ip{xk-i) around Xk, we get: 

dip(x) 



ip{xk-i) = ip{xk) - Ax 



dx 



Ax'^ d'^(f{x) 



dx'^ 



3! dx^ 



Ax4 d^^{x) 



4! 



dx^ 



Ax^ d^Lp{x) 



5! dx^ 



O(Ax^) 



The sum of these two equations gives: 

(^(xfc+i) + ip{xk-i) = 2ip{xk) + Ax- 



2 d'^vix) 



dx'^ 



Ax* d'^ip(x) 



12 dx^ 



0{Ax^ 



(A.22) 



(A.23) 



We use Eq. (A. 19) to define, directly, the second-order derivative. We compute the fourth-order derivative from a 



second-order derivative written as a difference quotient. Finally, we find a recursive expression for the wave function 
ip{xk-i): 



ip{Xk- 



[2 + f^Ax^f{xk)\ v{xu) 



1 



12 



f{xk + l) tp{Xk+l) 



1 



12 



(A.24) 



We have explicitly written Lp{xk-i) in Eq. (A.24) because later, in a numerical example, we will be interested 
in specifying the boundary conditions of ip{x) at x — )■ ±(X). The Numerov method |931 194j provides a numerical 
algorithm to solve ordinary differential equations of the second order in which the first-order term does not appear. 
It is a fourth-order linear multistep method. The above expression, Eq. (A.24), allows a simple solution of Lp{xk-i) 
once we know ip[xk), Lp{xk+i)^ f{xk), f{xk+i), and /(xfe_i). Finally, repeating the process, we can determine the 
entire solution. 

We will apply this last explicit computational algorithm for tunneling scenarios. For the numerical simulations. 



we will use Ax = 1 A. A gain, the amplitude of the global wave function ■0(x, t) = (p{x) 



-iEt/h 



is time independent. 



Thus, the current density J{x,t), from Eq. (34), is also time independent. Therefore, from the continuity equation, 
Eq. (28), we deduce that J{x,t) is also uniform (position-independent). Then, the velocity v(x) — constant/i?-^ (x) 
does not depend on time but depends on position. We have S{x,t) ^ 0. The velocity can be computed from: 



dx 



^2m{E-V{x)-Q{x)) 



(A.25) 



See problem 11 for a detailed calculation. This quantum expression, Eq. ( A.25 ), is different from the similar expression 
extracted from the classical Hamilton-Jacobi equation dS{x,t)/dx — ^2m(E — ^(x)). When E < V{x), we will not 
be able to find a (real) velocity in a classical system (classical particles are unable to pass through spatial regions 
where the potential barrier is higher than its energy). Interestingly, when E < V{x) in Eq. ( A.25[ ), we will still find 
(real) velocities because of the term Q{x) in Eq. ( |A.25[ ) (Bohmian particles are able to pass through spatial regions 
where the potential barrier is higher than its energy). This is the well-known tunneling effect. 

As a numerical example, let us concentrate on (time- independent) scattering states traveling from left to right in 
a double-barrier structure. For a simulation box < x < 40 nm, these states can be written as: 



(p(x) = (e*'^-^ + r(fc)e-*'=^) ; x < : 
V 27r 



(x) = {t{k)e'^^) ; X > 40 : 



(A.26) 
(A.27) 



where r{k) and t{k) are (unknown) complex coefficients accounting for reflection and transmission, respectively. 
However, up to an irrelevant (complex) constant that can be fixed by "normalization," we know the results <y3(x) at 
X — ±cx). In the central region, we can use the Numerov numerical method using the fact that the wave function and 
its derivative have to be continuous in the matching points between the regions where we have analytical or numerical 
knowledge of the wave functions. 

In Fig. |A.3[ we have represented one Bohmian trajectory for a state incident from left to right upon a double- 
barrier potential profile. Its energy E = 0.05 eV corresponds to the first resonant level of the double barrier [55]. 
In this particular case, where the transmission coefficient is very close to unity, the results obtained within the 
Bohmian approach are quite compatible with our intuitive understanding of the tunneling phenomenon: the particle 
is transmitted, and its velocity decreases in the well. We do also plot the total potential-that is, the sum V{x)-\-Q{x)- 
which is lower than the electron energy E = 0.05 eV. 
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FIG. A. 3. Bohmian trajectories associated with stationary scattering states impinging on a typical GaAs/AlGaAs RTD with 
2 nm barriers of 0.3 eV separated by a 7 nm well. The classical potential and the total potential (the sum of classical 
and quantum terms) are also depicted in dashed and solid lines, respectively, (a) Resonant eigenstate, Ek = 0.05 eV; (b) 
nonresonant eigenstate, Ek = 0.06 eV. Reprinted with permission from i95 . Copyright 1996 Elsevier. Abbreviation: RTD, 
resonant tunneling diode. 



However, the situation is far from intuitive for nonresonant states. In Fig. |A.3| 3, we have represented one of these 
nonresonant trajectories for the same potential profile. Now, although the transmission coefficient is much smaller 
than unity for the nonresonant states, they present the same features as the resonant ones: all Bohmian trajectories 
are transmitted through the barrier. This fact can be easily understood from a mathematical point of view. As we 
mentioned, the current and particle densities are time independent and positive everywhere. So, the Bohmian velocity 
is always positive. In this regard, although Bohmian trajectories associated with scattering states perfectly reproduce 
the presence probability and the current density, they do not reproduce our particle-intuitive understanding of the 
tunneling phenomenon, since in principle, we would expect reflected as well as transmitted trajectories. The origin 
of this lack of intuitive behavior of Bohmian trajectories is only due to the initial wave function that we select. The 
velocity of Bohmian particles is neither related to the incident electron nor related to the reflected one but to the 
average of both because the wave function itself is deflned as a weighted sum of waves with positive and negative 
momentum. In particular, if the barrier is inflnitely high, then the reflection coefficient is equal to one and Bohmian's 
velocity is zero (one would expect positive and negative electron velocities simultaneously). This argumentation is 
consistent with the first example: If the reflected plane wave is not present (the transmission coefficient equal to 
unity), then Bohmian trajectories reproduce our intuitive picture for tunneling. 

All the ID numerical algorithms previously discussed can be easily generalized to 2D and 3D systems (see, for 
example, the solution of the Schrodinger equation for computing Bohmian trajectories in atomic scenarios in Refs. 
[96L I97j). However, the computational burden of such techniques with finite-difference meshes is enormous, and 
the solution using adaptative meshes seems very interesting. For example, the triangular mesh generated using an 
algorithm inspired by Ref. [98 , with Delaunay refinement algorithms, yields a final mesh that is surprisingly well 
shaped, resulting in an excellent support to solve partial differential equations, as it was remarked in Ref. |99) . 



A. 1.2. Synthetic computation of Bohmian trajectories 



Up to here, we have only discussed analytical algorithms to find Bohmian trajectories. However, it is possi- 
ble to compute Bohmian trajectories without knowing the wave function. The key point is using the polar form 
ip{f,t) = R{f,t) e'}cp{iS{r,t)/h) to obtain the quantum Hamilton- Jacobi equations. As we anticipated, these new 
equations are nonlinear in R(f,t) and S{f,t). We will now discuss the advantages and disadvantages of numerically 
solving the quantum Hamilton- Jacobi equation, Eq. (39 1. Hereafter, we will analyze two different algorithms to solve 
numerically Eq. (39). First, directly solving these equations (the Eulerian formalism) or some equivalent variable 
change . Second rewriting the Hamilton-Jacobi equations within the Langrangian formalism and solving them. 
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The first point that we want to clarify is that, although the function S'(r, t) is in principle multivalued, the numerical 
solution of S{r, t) has a unique value at each point r and t. The reason is because of an important theorem on the 
existence and uniqueness of solutions to first-order equations with given initial- value problems. See page 70 in [lOOj . 



A. 1.2.1. Time- dependent quantum Hamilton- J acobi equations (TDQHJE\d) with an implicit (Newton-like fixed Eulerian 

mesh) method 



In the Eulerian formalism, we use a fixed grid. We are interested in providing a discussion valid for ID, 2D, and 
3D spatial scenarios. Therefore, we define a grid in real space through a set of discrete vectors for i = 1, . . . , iV-in 
particular, pi = (xi,yi,Zi) for a 3D space, pi = {xi,yi) for a 2D space, and pi = (xi) for ID. In order to provide a 
compact notation, we define: 

i?, = iR,{pi),...,R,ipN)f 



where Rj{pi) = R{r,t) at t 
a mesh in time tj = jAt for j = 1, 



tj and r 



,M 



(A.28) 

■,Sj{pN)f (A.29) 
V,^{Vj{p,),...,V,ipN)f (A.30) 
identically for Sj{pi) — S{f,t) and Vj{pi) = V{r,t)). We have also used 



In the definition of the temporal derivatives of the quantum Hamilton-Jacobi, we can use an implicit expression. 



Eq. ( A.5[), or an explicit one, Eq. (A.3|. The explicit expression will lead to an algorithm similar to the one described 
in Sec. A. 1.1.1 with the same advantages (simplicity) and disadvantages (i.e., conditional stable). In this section. 



we will explore the implicit route. Then, the discrete version of the quantum Hamilton-Jacobi equations, Eq. (39 1, 
reads: 



At 



At 



= 



(A.31) 



2m 



2to i?, 



The operations such as power, multiplication, and division between the vectors Xj and Yj are understood over each 
point in the spatial mesh, that is: 



[X,ip,)Yj{p,), Xj{pn)Y,{pn)]' 



(A.32) 



Here the symbol V is the matrix representation of the differential operator V in the ID, 2D, or 3D spatial mesh. 
The representation of V depends in general on the choice of the spatial mesh and on how the spatial derivatives are 
approximated. For example, in a ID system with a uniform mesh step size Ax, the discrete version of the spatial 
derivative in one dimension is: 



dRjx, t) __ Rj{pi+i) - R-jjpi^i) 
dx ~ 2Ax 

with Pi = xq + iAx. Then, the nabla operator in one dimension Vid is represented by the following matrix: 



(A.33) 



/ 



Viz 



2Ax 




-1 



V 



-1 
-2 2/ 



(A.34) 



Notice that we have changed the definition of the spatial derivative from Eq. (A.33) in xi and xj^. In two or three 



dimensions the matrices still have many zeros with the entries 1 and —1 properly located. 

Hereafter, as we mentioned, we only discuss the computational problems of the implicit solution of the Hamilton- 
Jacobi equations. Such an equation is nonlinear with respect to Rj and Sj, and the matrix inversion explained in 
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Sec. A. 1.1. 2 cannot be done. Then, in order to find the numerical solution of the discrete version of Eq. (A. 1.2.1 1, 



we have to employ some Newton-like method. If we rewrite in compact form the Eq. (A. 1.2.1) as: 

F{R„S,) = 



(A.35) 



then wc have to find the roots (zeros) of the nonlinear expression in Eq. (A.35 1. It is proved that the sequence: 



ih+i) 



R 



{h) 
ih) 



(A.36) 



where JpiR'j'^K is the Jacobian of F {r'j'^\ Sj'^^ ) , for suitable initial vectors R.'j^^ and Sj^\ converges quadratically 



However evaluating explicitly the Jacobian: 



to the solution of Eq. (|A.35|) for h — > cxi. Thus the key of this approach is to solve the linear system in Eq. (A.36 1 

diag{Rj) + V (^diag ^ ( \7 \ 

JFiRj,Sj) ' 



' 2m 



d{ag (4-) 



diag 



where: 



diag{Xj) 



\ 








Vdzag (^^ j V 
^^diag{Sj) + diag V 



(A.37) 



Xj{pN-i) 







(A.38) 



X.ipN) / 



we can readily recognize that for small components of Rj, the Jacobian is ill conditioned and the system of Eq. 
(A.36 1 is impossible to be solved numerically. Unfortunately this is the case of localized particles where Rj become 
exponentially small far from the central position of the particles. This problem does not appear with the explicit 
algorithm for solving the quantum Hamilton- Jacobi equations. The problem can also be avoided by choosing the 
variable change ip{r,t) = exp(C(r, t) + iS{f,t)/h) allowing to the equations: 



dt 



d_ 
Ft 



Sir,t) 



C(f,t) 



h^S{r, t) + 2V5(f, t) ■ VC(r, t)) 
Ira \ J 



(A.39) 



(V5(r,t))^ 
2to 



+ 1^ (v^C(r, t) + (VC(r, t)f) - ^ (r, t) 



However, the advantage of this method carries out also its disadvantage. The limitation of this approach is that 
the functions S(f^ t) and C'{f, t) strongly vary in the entire simulation space. Thus when we account for a realistic 
situation for which a large simulation space is used, the computational problem can become very hard because the 
number of space mesh points pk will be quite large to obtain an accurate solution. 



A. 1.2. 2. Time- dependent quantum Hamilton-Jacobi equations (TDQHJEid) with an explicit (Lagrangian mesh) method 

A different approach to solve the quantum Hamilton-Jacobi equation is based on the Lagrangian picture. The 
basic idea is to use as a grid the Bohmian trajectories themselves. In this picture, the point pk is no longer a fixed 
spatial point but a time-dependent point pk — Pk(t) with k = 1, . . . , N moving with instantaneous velocity Vk{t)- For 
a detailed analysis see Refs. |101fll04] . 

To pass from the Eulerian to the Lagrangian hydrodynamic picture, we must be careful with the time derivative of 
the functions R{pkit),t), S{pk{t),t) and V{pk{t),t). We have to use the chain rule: 

i^i+v{r,t)-y (A.40) 
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which inchides the convective term v{f, t) ■ V with the velocity defined: 

g{r^^t) = — VS'(f,t) 



(A.41) 



The new equations derived from the quantum Hamihon-Jacobi equations (or equivalently from Eq. (A. 1.2.1 )) can be 
written as: 



dt 



^Pfe(i) = v{pk{t),t) 



2m 



2m R{pk{t),t) 



ViPk{t),t) 



(A.42) 



(A.43) 



v{pk{t),t) = -^S{pk{t),t) 
m 



(A.44) 



dt 



RiPk{t),t) 



-R{pk{t),t)V -vipkit)) 



(A.45) 



The differential system in Eqs. 



( A.42 A.45[ ) is written for the trajectory pk{t) with a particular initial position p/c(0). 
However to (numerically) reproduce the wave function ipi^y t) = ^(^i t) exjp(iS{f, t)/^) defined in the entire space, we 
have to deal with fc = 1, . . . , iV trajectories with N oo. The distribution of their initial positions Pk{0) has to be 
obtained according to the sq uare m odulus of the initial wave function \ipipkiO), 0)^, that is, the quantum equilibrium 
hypothesis described in Sec. IV. 6 See also Ref. [20]. Thus, by construction we are using a mesh that is continuously 



adapting to the wave function time evolution because the Bohmian trajectories themselves follow the wave function. 

In this picture the strategy to solve Eqs. ( A.42| - A.45[ ) can be very complicated |101H104] . However, some simple 
consideration can be done in order to understand the limits of this approach. The method widely used [10111104] 
to integrate Eqs 



(A.42 A.45) is to solve them in cascade with the following approximations: 



at time tj the new 



trajectories positions are computed through: 



Pj,k = Pj-i.k + Atu,-_i, 



(A.46) 



with pj^k = Pk{ij)- In the updated grid defined by Pj,k, the variables Sj{pj_k), ^j{Pj,k), and Rj{pj,k) are evaluated by: 



SjiPj.k) = Sj-i{pj^k) + Ai 

RjiPhk) = Rj-liPj.k) - At 
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(A.47) 

(A.48) 
(A.49) 



where Vj is the gradient operator evaluated in pj^k through efficient interpolation methods [101H104j . It can be realized 
that we have used an (unstable) explicit (forward Euler) time integration method [105[ to avoid the computational 
burden associated with implicit methods with nonlinear equations discussed before. This limitation (small temporal 
step to ensure stable solutions) is made worse since the time step must be further reduced when some trajectory 
approaches points for which the potential V and/or the quantum potential Q strongly varies provoking large variations 
of some Vj. 

The variables 5* and v vary in the entire simulation space; thus when we account for a realistic situation such that 
a large simulation space is used, the computational problem can become very hard due to the fact that the number of 
space mesh points pk can diverge to obtain an accurate solution. On the contrary, this method has some fundamental 
advantage over methods using fix grid points. In the latter, the extension of ID algorithms to 2D and 3D algorithms 
means increasing exponentially the number of grip points. However, in the former methods (i.e., when the grid points 
are Bohmian trajectories), the number of points can be chosen by the user. Thus, for exemple, one can simulate a ID 
system with 1,000 Bohmian trajectories (grid points) and a 2D system with 3,000 Bohmian trajectories (grid points). 



A. 1.3. More elaborated algorithms 



We have presented very simple and easily implementable codes for computing Bohmian trajectories in simple 
scenarios. These algorithms can be used for newcomers to be introduced to the computation of quantum dynamics 
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with Bohmian mechanics. Certainly, there are many more elaborated and "professional" computational algorithms 
in the literature. An excellent source of information on this issue can be found in the book written by R. E. Wyatt 
entitled Quantum Dynamics with Trajectories: Introduction to Quantum Hydrodynamics [T^. 

To end this section, we want to mention briefly some original and powerful techniques, inspired by Bohmian 
mechanics, available in the literature, whose explicit description is far from the scope of this introductory appendix. 
Recently, B. Poirier showed how Bohmian trajectories can be computed for ID eigenstates without the wave function 
|107j . The work of N. Makri shows the possibility of using imaginary time algorithms to estimate eigenstates and 
eigenvalues |109j . With a similar goal, I. P. Christov uses quantum trajectories to develop time-dependent quantum 
Monte Carlo algorithms [1101 llllj . We also mention the work of E. L. Bittner |112| developed in chapter 5 in this 
book. The work presented in Sec. |V.6| can also be interpreted as an approximating technique for computing many- 
particle systems |81j . This last work has been used by the group of Oriols to develop a powerful quantum electron 
transport simulator. ( http:/7europe.uab.es/bitlles). See Refs. [11411115] or chapter 6 for practical results on electron 
correlation in nanoelectronics. Another line of work is also devoted to the approximation of the quantum potential. 
See, for example, the work of S. Garashchuk and V. A. Rassolov |113j . The correlation between classical and quantum 
systems has also been studied extensively. See, for example, Ref. jll6j . We do also highlight the work by D. J. Tannor 
and coworkers on the use of imaginary quantum trajectories (complex action) for quantum dynamics [1061 1108] . 

All the works mentioned have been elaborated during the past 10 years after the pioneering work of Wyatt [1011 
I117j . Much work is still needed to evaluate the real capability of synthetic applications of Bohmian mechanics as a 
computational tool. 
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